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CHAPTER  I 

INTRODUCTION 

A.  General  Statement  of  the  Problem 

Carbon,  in  its  elemental  or  allotropic  forms,  plays  a 
very  important  part  in  our  everyday  existence.  In  the 
form  of  coal,  it  is  still  one  of  the  major  energy  sources 
in  the  world.  Its  uses  range  from  being  a  moderator  (in 
the  form  of  graphite)  in  a  nuclear  reactor  to  that  of 
printer's  ink  (in  the  form  of  lampblack).'*'  Its  study  as 
a  chemical  element  forms  the  basis  of  organic  chemistry 
and  its  ever  widespread  use  in  industrial  applications 
continually  demands  a  better  understanding  of  its  properties. 

One  related  study  of  considerable  interest  is  that  of 
graphite  oxidation  and  erosion.  Graphite's  remarkable  re¬ 
sistance  to  chemical  attack  under  severe  conditions  has 
always  made  it  an  indispensable  material  in  the  electrolytic, 
smelting  and  foundry  industries.  In  the  aerospace  industry, 
the  search  for  highly  refractory  materials  for  aerospace 
and  interplanetary  travel  applications,  such  as  re-entry 
heat  shields  on  nose  cones,  nozzle  liners  for  rocket  engines 
and  thermal  protection  space  suit  clothing,  has  also  cast 
graphite  and  its  various  forms  as  an  attractive  material 

due  to  its  high  temperature  capabilities — large  heat  capacity 

2 

in  sublimation,  inertness  and  high  strength-to-mass  ratio. 
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In  most  cases,  of  interest  is  a  prediction  of  the  mass  loss 

rate  of  the  surface  material  for  a  given  external  atmosphere 

in  order  to  estimate,  for  example,  electrode  wear  or  a  heat 

shield  requirement.  A  precise  knowledge  of  gas-surface 

reaction  kinetics  would  then  be  ideal. 

Typically  the  atmosphere  of  interest  is  that  of  air, 

and  the  primary  reactions  involve  oxygen  and  nitrogen.  For 

a  re-entry  spacecraft,  the  dominant  mass  loss  of  a  graphite 

heat  shield  is  caused  by  chemical  reactions  with  the  hot 

shocked  gas  flowing  around  the  vehicle  for  surface  temperatures 

-2  2 

between  600°C  and  2500°C  and  pressures  in  the  10  to  10 

3 

atmospheres  range.  Sublimation  proves  to  be  the  primary 
mass  loss  mechanism  for  surface  temperatures  above  2500°C. 

For  the  chemical  reaction  dominated  temperature  range,  which 
does  cover  a  fairly  wide  range  of  interest,  reaction  with 

4 

oxygen  is  dominant  as  compared  with  nitrogen.  Furthermore, 
such  reactions  may  involve  dissociated  oxygen  generated 
downstream  of  the  shock  wave  preceding  the  re-entry  vehicle, 
even  though  dissociation  levels  in  the  undistrubed  air  are 
negligible  below  altitudes  of,  say,  200,000  ft.  Alternatively, 
a  vehicle  may  pass  through  atmospheric  regions  in  which 

5 

dissociation  has  been  induced  by  solar  radiation.  Thus, 
a  study  of  graphite  oxidation  should  desirably  include  the 
combined  but  separable  effects  of  atomic  and  molecular 
oxygen,  i.e.,  in  a  dissociated  atmosphere,  in  view  of  its 
prevalence  in  many  situations. 
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Examination  of  the  present  available  experimental  data 
on  graphite  oxidation  in  the  reaction  rate  controlled 
regime  shows  that  the  oxidation  process  always  is  assumed 
to  follow  a  rate  law  in  the  Arrhenius  form:** 

"eA/RTT 

m  =  ke(pT)  e  A  (1.1) 

(Explanation  of  all  symbols,  unless  otherwise  noted,  is 
given  in  the  List  of  Sumbols  and  Abbreviations  Section.) 
Depending  upon  the  investigator,  the  type  of  graphite 
studied,  and  the  test  conditions,  the  order  of  the  reaction, 
n,  has  been  reported  in  the  range  zero  to  unity;  the 
activation  energy ,EA , falls  within  the  limits  9  to  60  Real/ 
mole;  and  the  effective  collision  frequency,  ke,  varies 
over  several  orders  of  magnitude. 

Several  explanations  can  be  offered  for  the  huge 
scatter  in  the  experimental  data.  One  concerns  the  validity 
of  the  Arrhenius  rate  law  itself.  The  rate  expression  is  a 
carryover  from  homogeneous  reactions  and  some  questions  can 
be  raised  about  its  applicability  to  heterogeneous  reactions. 
Though  the  exponential  component  of  the  equation  may  be 
justified  by  concepts  from  kinetic  theory  as  applied  to  two 
reacting  gases**  (treating  the  surface  as  an  array  of  atoms, 
for  example) ,  no  distinction  is  then  made  between  gas  and 
surface  temperatures,  which  often  is  not  the  case.  Constant 
activation  energy  was  usually  assumed  for  macroscopic 


7 

applications  in  the  temperature  range  of  interest  but,  in 

fact,  experimental  values  do  show  a  dependence  on  tempera- 

8  9 

ture.  For  example,  both  Blyholder  and  Eyring  and  Horton, 
in  their  respective  experiments,  show  that  logarithmic  plots 
of  m  versus  1/T^,  (so-called  Arrhenius  plots)  are  nonlinear, 
imp lying  a  "varying"  activation  energy.  Furthermore,  the 
reaction  order  was  observed  to  be  temperature  dependent. 

Rosner  and  Allendorf,1®  in  their  experiments  on  graphite 
oxidation,  have  also  shown  that  oxidation  probability  (a 
measure  of  m)  in  an  Arrhenius  plot  is  non-linear.  The 
effective  collision  frequency,  ke,  is  essentially  an 
empirical  steric  factor  and,  of  course,  does  not  offer 
insight  into  the  interaction  details. 

A  second  major  reason  for  extensive  scatter  of 
collected  data  is  related  to  the  experimental  techniques. 

A  survey  of  the  literature  on  graphite  oxidation  rate 
studies,  a  brief  summary  of  which  is  presented  in  Section  I-B, 
reveals  that  a  majority  of  the  experiments  were  carried  out 
in  the  continuum  regime  (small  Knudsen  number  based  on  the 
mean-free-path  of  the  reacting  gas  and  a  characteristic 
dimension  of  the  macroscopic  surface)  in  which  diffusion 
effects  may  be  important.  Multiple  gas-surface  collisions 
as  well  as  gas-gas  collisions  are  then  common,  the  implica¬ 
tion  being  that  the  species  detected  were  not  necessarily 
the  primary  products  of  the  surface  reaction  but  rather  a 
near  equilibrium  or  possibly  equilibrium  distribution  of 


xw*‘  r.»f  ^V*'*1*^ 
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the  products.  Very  reactive  free  radical  intermediates  of 
the  reaction  were  always  masked  and  never  detected.  Since 
the  primary  goal  of  all  kinetic  studies  is  the  delineation 
of  a  reaction  mechanism  consistent  with  kinetic  delta,  a 
somewhat  more  accurate  and  controllable  experimental 
technique  is  in  order.  It  was  the  purpose  of  the  present 
work  to  devise  such  an  experiment  so  as  to  permit  study  of 
the  graphite  oxidation  process  in  a  dissociated  atmosphere. 

B.  Summary  of  Past  Carbon  and  Graphite  Oxidation 

Inves  tigations _ 

The  continuous  interest  in  carbon  and  graphite  oxida¬ 
tion  over  many  decades  had  generated  a  voluminous  amount  of 
literature  on  the  subject.  A  literature  survey  reveals 
that  almost  every  common  form  of  carbon  has  been  the  subject 
of  oxidation  studies.  Such  work  naturally  included  both 
scientifically  oriented  studies  in  which  the  reaction 
mechanisms  were  of  primary  interest  and  engineering  design 
type  studies  in  which  measurements  of  erosion  of  devices 
fabricated  from  particular  forms  of  carbon  were  aimed  at. 

Only  the  former,  more  fundamental,  work  will  be  reviewed 
here  and  even  so,  an  exhaustive  review  will  not  be  attempted. 
Instead,  what  follows  will  be  a  presentation  of  representa¬ 
tive  work  and  of  major  contributions.  Readers  interested 

in  more  discussions  are  referred  to  articles  by  Blackwood^ 

12 

and  Walker  et  al  for  work  done  prior  to  1959  and  by 
Thomas^  for  more  recent  investigations. 
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The  crystallography  of  graphite  is  referred  to  any 

standard  handbook,  such  as  Ref.  1. 

14 

Smith,  in  1863,  was  probably  the  first  to  suggest 
that  when  o'.vgen  is  adsorbed  on  a  carbon  surface  it  under¬ 
goes  a  chemical  change.  He  found  that  adsorption  of  oxygen 
on  charcoal  did  not  cease  even  after  a  month;  and  that 
while  nitrogen  and  other  gases  adsorbed  on  charcoal  could 
be  easily  removed  as  such,  oxygen  could  be  reiaoved  only 
on  strong  heating,  and  then  only  as  carbon  dioxide.  Rhead 
and  Wheeler15  concluded  that  oxygen  combines  with  a  mass 

of  carbon  directly  to  form  a  physicochemical  complex,  CO, 

x  y 

of  variable  composition ,  which  decomposes  over  a  wide 

range  of  temperature,  giving  a  mixture  of  carbon  monoxide 

CO  and  carbon  dicxide  C02.  These  observations  were  con- 

16 

firmed  by  Lowry  and  Hulett,  who,  while  attempting  to 
determine  the  content  of  the  moisture  and  gas  in  charcoal, 
noticed  that  the  gases  evolved  up  to  900°C  consisted  almost 
entirely  of  C02  and  CO,  the  former  predominating  at  lower 
temperatures  and  the  latter  at  higher  temperatures. 

Further  evidence  for  chemisorption  or  fixation  of 

17 

oxygen  on  carbons  were  furnished  by  the  work  of  Langmuir 
and  others.18'15  Langmuir,  in  his  work  on  the  interaction 
of  oxygen  with  carbon  filaments,  noticed  that  the  gas  was 
held  so  firmly  on  the  carbon  surface  that  it  came  off  when 
heated  only  as  oxides  of  carbon.  He  also  found  that  the 
adsorption  rate  of  oxygen  on  the  carbon  filament  was  quite 
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irreproducible  and  depended  on  the  previous  history  of  the 

filament.  Nevertheless,  he  proposed  that  on  the  surface 

of  the  filament  one  adsorhable  atom  or  molecule  of  the  gas 

would  be  bound  for  each  carbon  atom.  Combination  at  the 

surface  then  produced  CO  and  CO2 ,  which  were  desorbed  when 
20 

heated.  Shah,  using  charcoal  as  a  source  of  carbon, 
proposed  the  chemisorption  of  oxygen  on  the  surface  and 
suggested  that  the  mechanism  was  similar  to  that  put 
forward  by  Langmuir. 

By  then,  it  was  quite  well  established  that  when 
carbon  was  exposed  to  oxygen  three  products  could  form: 
gaseous  carbon  monoxide,  gaseous  carbon  dioxide,  and  a  carbon- 
oxygen  surface  complex.  Also,  the  amount  of  gaseous  products 
formed  relative  to  the  amount  of  surface  complex  formed  was 
shown  to  depend  strongly  upon  temperature.  Though  the 
exact  nature  or  composition  of  the  surface  complex  was  not 
known,  the  concept  as  such  was  an  important  step  in  under¬ 
standing  the  reaction  mechanisms. 

Studying  the  oxidation  of  thin  layers  of  activated 
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charcoal,  Chukhanov  and  Karzhavina  concluded  that  oxida¬ 
tion  of  the  carbon  could  take  place  at  a  temperature  as  low 
as  300 °C  and  that  it  proceeded  through  the  formation  of  a 
surface  oxide,  C-jO^,  which  formed  very  rapidly  and  then 
decomposed  in  the  presence  of  carbon  and  oxygen,  giving  CO 
and  CO2  in  equal  amounts.  At  temperatures  above  1100°C 
the  reaction  favored  the  formation  of  CO.  Increasing  the 
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velocity  of  the  gases  flowing  over  the  carbon  suppressed  the 

formation  cf  C02‘  The  fact  that  the  ratio  cl '  CO  to  C02  was 

not  constant  and  could  be  varied  by  a  change  in  the  condi - 
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tions  was  confirmed  by  Tsukhanova. 

Realizing  that  the  oxidation  of  a  carbon  surface  had 

been  shown  to  be  a  complicated  process  and  the  possibility 

that  primary  steps  could  be  masked  by  secondary  px/<ces3es 

23 

in  previous  experiments,  Meyer  studied  the  oxidation  of 
graphite  in  the  form  of  electrically  heated  elements  at 
low  pressures.  He  found  that  the  pressure  in  the  reaction 
cell  had  * o  be  low  enough  for  gas  molecules  leaving  the 
hot  surface  to  reach  the  cold  walls  of  the  cell  without 
collision  with  other  molecules  and  the  subsequent  formation 
of  secondary  reaction  products.  He  also  found  that  the  gas 
velocity  through  the  apparatus  had  to  be  great  enough  to 
prevent  the  reaction  products  returning  to  the  filament  and 
that  the  graphite  surface  had  to  be  smooth  to  avoid  such 
occurrences.  CO  and  C02  were  observed  in  equal  proportions 
independent  of  temperature  and  pressure  at  temperatures 
below  1200°C.  The  rate  of  reaction,  however,  was  propor¬ 
tional  to  the  oxygen  pressure.  The  reaction  was  explained 
by  supposing  that  the  oxygen  penetrated  deeply  into  the 
graphite  lattice  by  solution  and  at  the  same  time  passed 
into  a  reactive  state  and  disintegrated  the  lattice, 
preferentially  at  the  basal  planes.  As  the  temperature  was 
raised  above  1200 °C,  the  reaction  rapidly  slowed  down  and 
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the  CC^/CO  ratio  fell  to  approximately  0.5  at  about  1500°C, 
the  reaction  then  being  of  the  zeroth  order.  This  change 
was  attributed  to  the  diminution  of  the  solubility  of  the 
oxygen  in  the  carbon  until  ultimately  the  surface  oxidation 

mechanism  became  solely  responsible  for  the  reaction  rate. 

24 

The  suggestion  by  Mayers  that  the  rate  of  oxidation 

of  carbon  under  many  conditions  war.  determined  not  by  the 

speed  of  chemical  reaction  but  by  the  rate  at  which  oxygen 

could  diffuse  onto  the  surface  helped  to  explain  dis- 
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crepancies  in  several  obsrrvad  rates.  Kuchta  et  al, 
using  high  velocity  air  streams  passing  over  spectroscopic- 
grade  carbon  between  900°C  and  1200°C,  were  able  to 
demonstrate  the  change  of  the  reaction  rate  from  diffusion- 
limited  to  chemically-controlled  as  the  air  velocity  was 
varied.  The  reaction  of  carbon  and  oxygen  was  diffusion- 
limited  at  low  air  velocities,  but  at  high  air  velocities 
the  rate  of  reaction  approached  that  predicted  on  the 
basis  of  the  Arrhenius  equation. 

26 

Strickland-Cons table1 s  earlier  work  in  1944  showed 
substantial  agreement  with  recent  results  obtained  in  the 
study  of  graphite  oxidation.  Using  filaments  of  fully  out- 
gassei  graphite  resistively  heated  to  900-2000°C  and  burned 
in  an  oxygen  background  of  20-500  millitorr,  he  found  that 
CO  was  the  principal  reaction  product.  There  was  always  a 
trace  of  C02  but  the  CO/CO^  ratio  was  invariably  greater 
than  7  and  varied  inversely  with  temperature.  He  could  not 
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determine  an  origin  for  the  C02  as  the  snail  amounts  could 
have  resulted  from  secondary  oxidation  of  CO.  The  reaction 
was  observed  to  obey  fir'-t  order  pressure  dependence  with 
a  maximun  rate  occurring  at  about  1200°C.  The  rate  of 
the  reaction  was  also  found  to  exhibit  a  hysteresis 
behavior,  the  extent  of  which  depended  on  the  past  tempera¬ 
ture  history  of  the  filament.  On  proceeding  from  a  lower 
to  a  higher  temperature  a  new  rate  lower  than  the  steady 
state  value  was  obtained.  Conversely,  on  reducing  tempera¬ 
ture,  the  initial  rate  was  found  to  be  higher  than  the 
steady  state  rate.  Such  hysteresis  was  explained  by  the 
existence  of  an  inhibiting  surface  layer  which  was  removable 

at  high  temperatures.  Similar  results  were  obtained  by 
27 

Duval,  again  using  a  hot  filament  technique,  but  a 

reverse  hysteresis  was  noted.  This  reverse  hysteresis  from 

that  of  Strickland-Cons table  was  confirmed  by  recent  re- 
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suits  of  dander  et  al  and  Shih  using  completely 

different  experimental  techniques.  Strickland-Constable , 

in  addition,  noticed  that  the  amount  of  oxygen  adsorbed 

was  somewhat  less  than  previously  suggested. 

A  study  by  van  Loon  and  Smeets^0  of  burning  carbon 

particles  in  a  vortex  chamber  had  tended  to  confirm  the  view 

that  CO  was  the  primary  product  but  the  controversy  was  by 
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no  means  over.  Meyer  showed  that  the  results  of 
Strickland-Constable  should  not  be  taken  without  reservation 


mainly  because  that  the  latter  was  unable  to  prevent 
multiple  collisions  with  the  hot  filament  due  to  che  lack 
of  a  flow  technique.  In  contrast  to  Strickland-Constable, 
he  found  the  carbon-oxygen  reaction  to  be  of  first  order 
from  700-1200°C  while  from  1500-22G0°C  a  zeroth  order  was 
observed.  No  maximum  in  the  reaction  rate  was  observed. 
Instead/  the  two  regions  of  first  and  zerot'.  order  pressure 
dependence  were  characterized  by  constant  activation 
energies  of  25  and  90  kcal/mole  respectively.  The  CO/CC^ 
ratio  was  nearly  unity  and  was  independent  of  temperature 
and  pressure.  He  attributed  the  zeroth  order  dependence 
at  high  temperatures  to  evaporation  of  reaction  product 
as  being  the  rate  limiting  step. 
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Oxidation  of  pyrolytic  graphite  was  studied  by  Horton 
in  a  flowing  system  at  atmospheric  pressure.  By  measuring 
dimensional  changes,  the  reaction  rates  were  obtained  for 
the  basal  and  edge  planes.  Neither  gas  velocity  nor 
diffusion  greatly  affected  the  observed  rates,  which  were 
of  one-half  order  with  an  activation  energy  of  37  kcal/rr.ole 
at  temperatures  between  600-1550 °C.  A  mechanism  consisting 
of  rapid  adsorption  of  oxygen  with  rate  limiting  desorption 
of  CO  or  CO2  was  suggested.  Also,  the  reaction  rate  of 
spectroscopic  graphite  was  found  to  be  two  orders  of 
magnitude  faster  than  that  of  the  basal  plane  of  pyrolytic 
graphite. 
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Many  attempts  '  '  had  been  made  to  study  the  kinetics 
of  the  reaction  by  determining  the  order  of  the  reaction 
and  the  energy  of  activation.  The  reaction  appeared  to  be 
of  the  first  order  up  to  temperatures  of  approximately 
1300 °C,  but  above  this  was  zeroth  order.  Energies  of 
activation  varied  between  14  and  36  kcal/mole.  The 
difficulties  involved  in  preventing  secondary  reactions 
and  the  dependence  of  the  reaction  order  on  the  surface 
coverage  by  the  carbon-oxygen  complex,  which  in  turn  was 
a  function  of  the  absolute  pressure  level,  c  illed  for  some 
reservation  in  accepting  the  results  of  those  investigations. 
The  study  of  the  reaction  order,  however,  did  provide  a 
means  for  estimating  the  percentage  of  surface  area  covered 
by  stable  surface  complex. 

34 

Loebenstein  and  Deitz,  in  their  study  of  the  chem¬ 
isorption  of  oxygen  on  amorphous  carbon  at  2')0°K,  indicated 

that  less  then  6%  of  the  total  carbon  surface  proved  active 
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for  adsorption.  Bonnetain,  Duval  and  Letort,  by  sub¬ 
jecting  carbon  to  oxidation  and  succes  ive  degassing  at 
temperatures  between  450  and  700°C,  found  that  the  relative 
production  rates  o*  CO  and  C02  depended  on  pressure  levels 
as  well  as  temperatures.  They  concluded  that  the  simplest 
model  accounting  for  these  effects  would  be  based  on  a 
surface  complex  state  and  estimated  that  the  surface  com¬ 
plexes  covered  about  5%  of  the  total  surface  area  at  a 

-4 

pressure  level  of  10  torr. 
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Laine  et  al,  monitoring  the  course  of  the  oxygen- 
graphitized  carbon  (Spheron  C)  reaction  mass  spectro- 
metrically,  determined  that  CO  was  the  principal  product 
of  the  surface  complex.  Their  outgassing  experiments  also 
yielded  the  total  active  surface  area  which  at  300°C  was 
about  2.3%  of  the  total  area.  Due  to  the  change  in  the 
available  area  with  coverage,  the  reaction  was  found  to 
deviate  from  first  order.  The  surface  complex  was  found 
to  be  stable  at  the  reaction  temperature  but  decomposed 
and  desorbed  as  the  surface  temperature  was  increased 
with  nearly  all  of  it  desorbed  by  S50°C. 

Strickland-Constable  et  al  r~  studied  the  oxidation 
rate  of  pyrolytic  graphite  between  1000-2400°C  by  means  of 
a  directed  oxygen  jet  at  a  heated  carbon  bar  and  measured 
the  rates  from  photographed  erosion  profiles.  They  were 
able  to  correlate  their  results  empirically  with  modifica¬ 
tions  of  the  theory  proposed  earlier  by  Blyholder,  Binford 
39 

and  Eyring.  The  modified  theory  assumes  that  surface 
reaction  sites  can  be  classified  into  two  groups:  a  more 
reactive  A-site  and  a  less  reactive  B-site.  At  low 
temperature  the  oxidation  rate  is  dominated  by  the  more 
reactive  A-sites,  which  represent  a  fraction  of  the  surface. 
At  higher  temperatures  the  A-sites  rearrange  thermally  to 
form  less  reactive  B-sites.  A-sites  are  also  created  when 
B-sites  are  oxidized  and  desorb  CO.  It  is  the  competition 
between  the  generation  of  active  A-sites  by  oxidation  of 
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B-sites  and  their  deactivation  by  thermal  annealing  that 

resulted  in  the  observed  maximum  in  the  reaction  rate  with 

varying  graphite  temperatures.  At  still  higher  temperatures 

the  oxidation  rate  increases  with  an  increase  in  the 

B-site  reactivity.  Though  the  explicit  nature  of  the  two 

different  sites  was  not  stated,  it  was  encouraging  to  see 

that  experimental  results  over  a  wide  range  of  pressures 

and  temperatures  fitted  quite  well  to  this  two-site  model. 

In  order  to  determine  the  nature  of  the  active  sites 

on  which  surface  complexes  could  form,  topographical  studies 

by  means  of  optical  and  electron  microscopy  have  been  made 

for  the  carbon-oxygen  reaction.  Perhaps  best  known  in  this 
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area  is  the  work  of  Hennig.  '  '  Using  naturally 

occurring  graphite  single  crystals  as  well  characterized 

specimens,  he  found  that  the  edge  planes  oxidized  at  least 

twenty  times  faster  than  the  basal  planes.  The  etch  pits 

on  the  basal  planes  were  hexagonal  in  shape  bounded  by 

vertical  (1120,  1210,  2ll0)  faces,  suggesting  that  the  more 

active  sites  were  lattice  defect  points,  such  as  vacancies 

43  44  45 

and  dislocations,  and  grain  boundaries.  Thomas  et  al  '  ’ 

drew  similar  conclusions  using  optical  microscopy.  In 

45 

addition,  Thomas  found  that  atomic  oxygen,  generated  by 
microwave  discharge  of  oxygen,  formed  uniform- y  spaced  pits 
over  the  entire  basal  plane  surface  with  no  preferred  sites. 
No  quantitative  measure  of  the  reaction  rate  was  made. 


however. 
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The  suggestion  of  different  types  of  adsorption  sites 

led  to  the  adsorption  experiments  by  Hart,  Vastola  and 
50 

Walker.  They  found  two  levels  of  adsorption  in  the 

temperature  range  from  25-400°C.  Dissociative  immobile 

adsorption  was  assumed  and  suggestions  were  made  that 

different  types  of  adsorption  sites  might  be  associated 

with  the  various  carbon-carbon  distances  in  graphite. 

51 

Lussow  et  al  studied  the  chemisorption  of  oxygen  on 

graphon  in  the  temperature  range  from  300-6 80 °C.  Oxygen 

chemisorption  was  found  to  increase  rapidly  above  400 °C, 

indicating  two  types  of  active  sites.  A  suggestion  was 

made  that  an  oxygen  molecule  first  formed  a  mobile 

adsorption  complex  which  later  dissociated  to  form  two 

mobile  atomic  surface  complexes,  seemingly  in  contrast  to 
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the  assumption  of  Hart  el  al.  The  conclusion,  however, 
was  that  different  type  of  active  sites ,  dominating  in 
various  temperature  ranges,  might  have  different  fundamental 
behaviors . 
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Bansal  et  al,  using  preoxidized  graphon,  studied  the 
desorption  rate  of  surface  oxides  from  500-1000°C.  Graphs 
of  cumulative  CO  desorption  vs .  temperature  revealed  several 
distinctive  slopes  of  the  curve,  implying  that  CO  desorbed 
from  different  types  of  sites.  Also,  oxides  formed  at 
higher  temperatures  were  found  to  be  less  stable  than  those 
formed  at  lower  temperatures. 
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Walker  et  al  investigated  the  effects  of  previously 

adsorbed  surface  complexes  on  later  adsorption  patterns 

and  subsequent  desorption  products.  Their  approach  made 
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use  of  heavy  oxygen  (02  )  as  a  tracer.  A  graphon 

sample  at  280°C  was  first  exposed  to  ordinary  oxygen  and 
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then  subsequent  adsorption  was  made  using  02  .  Upon 
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desorbing  the  surface  complex  formed,  the  fraction  of  O 
in  each  product  was  found  to  be  invariant  with  desorption 
temperature ,  indicating  that  the  time  of  formation  of 
surface  complexes  had  little  effect  on  the  reaction 
process. 

The  thermal  annealing  or  deactivation  of  active 

sites  as  proposed  in  the  two-site  model  of  Strickland- 

Constable  et  al  may  be  attributed  to  the  reduction  of 
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crystal  imperfections  by  heat  treatment.  Bowman,  in 
his  study  of  the  imperfections  in  graphite  structure 
using  X-ray  diffraction,  noticed  that  most  strains  in  the 
crystal  structure  were  eliminated  by  heating  in  the 
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temperature  range  of  2000-3000 °C.  Nightingale  and  Snyder, 
in  addition  to  X-ray  diffraction  technique,  monitored  the 
dimensional  changes  of  the  bulk  graphite  alter  exposing  it 
to  neutron  radiation,  which  resulted  in  Jattice  defects 
such  as  vacancies  and  interstices.  It  was  found  that  higher 
temperatures  or  longer  time  of  annealing  had  a  parallel 
effect  in  the  rate  of  restoring  the  bulk  graphite  to  its 
original,  preradiation  dimension  and  structure.  Similar 
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conclusions  were  drawn  by  Price  and  Bokros5^  for  neutron 

irradiated  pyrocarbons,  primarily  based  on  measurements 

of  physical  properties  such  as  density ,  resistivity  and 

microhardness.  Specifically,  restoration  took  place 

after  one  hour  at  about  1400° C  with  a  rapid  increase  in 

recovery  rate  at  over  200 0°C. 

Oxidation  on  carbon  surface  was  known  to  cause  lattice 

inperfections  which  tend  to  accelerate  the  reaction  rate. 

In  fact,  this  process  was  responsible  for  the  generation 

of  A-sites  from  B-sites  as  proposed  in  the  two-site  model 
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of  Strickland-Cons table  et  al.  Antonowicz  discovered 

electron  spin  centers  on  carbon  surface  after  oxidation. 

This  was  attributed  to  the  removal  of  carbon  atoms  leaving 

neighboring  atoms  with  unpaired  spins.  As  expected,  the 

spin  centers  were  eliminated  upon  thermal  annealing. 

The  effects  ot  impurities  in  the  carbon  on  the 

reactivity  with  oxygen  had  also  been  investigated. 
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Harker  noticed  an  increase  in  reactivity  of  coconut 
charcoal  and  oxygen  with  the  presence  of  alkali  metals  but 
attributed  the  result  to  the  attack  of  the  surface  by  free 
metal.  Amariglio  and  Duval*  found  that  metal  oxides  of 
high  vapor  pressures  were  generally  catalysts  in  acceler¬ 
ating  the  graphite  oxidation  reaction  and  proposed  that  the 
metal  oxide  furnished  a  "reactive'’  oxygen  to  the  graphite. 
Roscoe  studied  the  catalytic  effects  of  boron  and 
molybdenum  on  oxidation  of  natural  graphite  crystals  and 
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concluded  that  the  catalyst  acted  on  existing  active  sites. 

To  prevent  the  impurities,  which  existed  in  most 

forms  of  carbon,  from  masking  the  true  kinetics,  Evans  and 
49 

Phaal  employed  diamond  as  a  clean,  well-defined,  non- 

porous  surface  in  studying  the  carbon-oxygen  reaction. 

Unfortunately,  at  temperatures  over  850°C  a  stable,  black 

carbon  film  formed  on  the  surface,  drastically  changing 

the  surface  characteristics  and  destroying  the  original 

idea  of  a  well-characterized  surface. 

Recently,  a  totally  different  approach  from  previous 

experimental  techniques  has  been  made  in  studying  graphite 
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oxidation.  Olander  et  al,  using  modulated  molecular 

beam-mass  spectrometry,  investigated  the  reaction  kinetics 

of  oxidation  of  pyrolytic  graphite  in  both  the  basal  and 

pricm  planes.  Suggestion  of  reaction  mechanisms  was 

based  on  comparison  of  model  calculations  and  experimental 

results  of  a  reaction  product  vector.  The  experiment 

covered  an  effective  oxygen  beam  pressure  range  at  the 
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target  between  5  x  10  and  2  x  10  torr  together  with  a 
surface  ten^erature  range  of  1000  to  1800°K.  CO  was  found 
to  be  the  dominant  reactant  product  and  a  first  order  de¬ 
pendence  of  pressure  was  observed  for  surface  temperatures 
above  1200°K.  The  reaction  probability  exhibited  a 

maximum  around  1450 °K  and  a  hysteresis  behavior  similar  to 
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that  observed  by  Duval  was  confirmed.  A  two-branch, 
two-site  mechanism,  similar  to  the  model  used  by  Strickland 
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Constable,  '  was  proposed  for  the  basal  plane  reaction 

in  which  the  molecular  oxygen  first  dissociated  on  the 

surface  forming  one  bound  CO  and  one  mobile  adsorbed  0  atom. 

The  adsorbed  atomic  oxygen  was  said  to  then  migrate  about 

the  surface  until  it  reacted  with  a  suitable  carbon  atom 

and  subsequently  both  bound  CO  were  desorbed.  Oxidation 

of  the  prism  plane  was  found  to  be  affected  by  the  diffusion 

of  the  surface  species  into  the  bulk  graphite.  While  the 

work  represented  an  elegant,  self-contained  theory  and 

experiment,  the  interpretation  of  data  depended  heavily 

on  the  modelling  phase.  Nevertheless,  it  signifies  an 

important  step  in  unraveling  the  carbon-oxygen  reaction. 

29 

Shih  carried  out  a  similar  molecular  beam  experiment 
in  which  the  re-emitted  particles  from  the  target  surface 
were  modulated  instead  of  the  incident  beam.  The  phase 
shift  measurements  were  avoided  but  so  were  their  subjec¬ 
tive  interpretations.  General  agreement  of  results  existed 
between  Shih  and  Olander  et  al  and  in  particular,  Shih 
showed  that,  depending  on  the  temperature  range,  several 
minutes  to  many  hours  were  required  to  reach  steady  state 
reaction  rates  with  regard  to  surface  annealing  and  site 

generation.  While  the  empirical  parameters  chosen  by 
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Nagle  and  Strickland-Cons table  in  their  two-site  model 
appeared  to  correlate  quite  well  for  results  of  incident 

CO 

C>2  pressure  greater  than  10  atm.,  the  same  parameters, 
extended  to  the  lower  pressure  range  of  Shih  (order  of 
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10  atm. ) ,  seemed  to  predict  a  lower  reaction  probability 

than  was  measured.  It  was  suggested  that  the  discrepancy 

arose  due  to  the  relative  insensitivity  of  the  chosen 

parameters  to  the  detailed  dynamics  on  the  surface  at  high 

-5 

pressures  (above  10  atm.)  as  implied  by  the  short  adjust¬ 
ment  times  while  this  was  not  true  for  the  case  of  low 
—  8 

pressures  (10  atm.  or  lower)  where  a  somewhat  more 
accurate  description  of  the  surface  processes  was  necessary. 
CO  was  also  observed  to  desorb  in  a  lobular  fashion  towards 
the  surface  normal  and  the  flux  could  be  described  quite 
well  by  a  cosine  square  distribution. 

While  numerous  investigations  have  been  carried  out 
for  graphite  oxidation  by  molecular  oxygen,  investigations 
of  the  effects  of  atomic  oxygen  on  the  rate  of  oxidation 
have  been  relatively  rare,  primarily  due  to  the  difficulty 
in  generating  atomic  oxygen. 

59 

In  1959  Streznewski  and  Turkevich  studied  the 
reaction  of  carbon  with  oxygen  atoms.  Evaporated  carbon 
films  and  carbon  particles  obtained  by  burning  benzene  in 
air  were  reacted  with  atomic  oxygen  produced  by  D.C.  dis¬ 
charge  at  5000  volts.  Optical  transparency  of  the  carbon 
films  was  monitored  for  the  rate  of  erosion  and  0  concen¬ 
tration  was  determined  by  the  rise  of  temperature  of  a 
platinum  filament  introduced  in  the  gas  stream.  The 
oxidation  rate  of  the  carbon  black  samples  was  inferred  by 
noting  the  change  in  size  using  electron  microscopy.  The 
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rate  of  erosion  of  carbon  was  found  to  be  independent  of 
the  specimen  temperature  in  the  range  from  20  to  100°C 
and  to  depend  directly  on  0  concentration. 

Vastola  et  al^^  reacted  spectroscopically  pure  graphite 
rods  at  room  temperature  with  0  from  a  microwave  discharge 
and  analyzed  the  products  of  reaction  mass  spectrometrically. 
CO  was  found  to  be  the  major  product.  However,  no  quantita¬ 
tive  results  were  attempted  at  all. 

Rosner  and  Allendorf'*''^  have  made  comparative  studies 
of  the  attack  of  pyrolytic  and  isotropic  graphite  by  atomic 
and  molecular  oxygen.  Atomic  oxygen  was  generated  by  a 
microwave  discharge  in  a  gas  flow  system  and  the  ablation 
rate  of  the  electrically  heated  graphite  filament  was 

-2 

monitored  optically.  Typical  oxygen  pressure  level  was  10 
torr  and  the  target  temperature  varied  from  1100  to  2000 °K. 
Their  principal  findings  were:  (1)  0  was  an  order  of 
magnitude  more  reactive  than  (>2  for  both  pyrolytic  and 
isotropic  graphite  at  equivalent  surface  temperatures; 

(2)  the  oxidation  rates  of  pyrolytic  and  isotropic  graphite 
by  O  were  comparable  whi le  pyrolytic  graphite  was  more 
oxidation-resistant  than  isotropic  graphite  in  O^;  (3)  all 
oxidation  probabilities  exhibited  a  maximum  in  the  tempera¬ 
ture  range  covered;  (4)  the  surface  temperature  dependence 
of  the  oxidation  probabilities  by  0  were  much  less  pronounced 
than  that  by  02.  These  results  implied  that  0  and  02  oxida¬ 
tion  may  very  well  proceed  along  different  mechanisms.  The 
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surface  structure  of  graphite,  which  proved  to  be  important 
for  the  02  oxidation  rate,  had  a  minor  effect  on  the  0 
oxidation  rate.  It  was  suggested  that  while  the  0 2  oxida¬ 
tion  rate  was  limited  by  dissociative,  activated  adsorption 
of  oxygen  molecules  on  the  surface,  the  0  oxidation  followed 
a  Rideal-type  mechanism  in  which  the  impact  of  a  gaseous  0 
upon  an  adsorbed  0  led  to  the  formation  and  desorption  of 
C02. 

Marsh  et  al,61'^2'®3  a  series  0f  investigations,  had 
made  careful  studies  on  the  surface  conditions  of  graphite 
after  exposing  it  to  a  stream  of  atomic  oxygen.  Again 
microwave  discharge  was  used  to  generate  O  and  its  concentra¬ 
tion  measured  by  calorimetric  method  using  a  catalytic  probe 
of  freshly  prepared  silver-coated  platinum  wire  which 
limited  the  maximum  gas  temperature  to  be  measured  at  about 
100°C.  The  total  weight  loss  of  the  specimen  was  measured 
by  a  vacuum  microbalance  and  the  surface,  structure  was 
examined  by  an  electron  microscope.  For  a  surface  tempera¬ 
ture  range  of  0  to  400°C,  it  was  found  that,  unlike  molecular 
oxygen  which  attacked  preferred  regions  of  the  cleavage 
surface  such  as  defect  points,  atomic  oxygen  produced  etch 
pits  over  almost  the  entire  basil  surface.  The  pits  were 
conically  shaped,  smaller  in  size  but  greater  in  number 
than  the  well-defined,  though  somewhat  shallow,  hexagonal 
pits  produced  by  molecular  oxygen.  Additionally,  the 
atomic  oxygen  tended  to  display  preferential  attack  of  the 
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basal  planes  along  twin  bands.  A  zero  or  very  small 
activation  energy  was  observed  for  the  C-0  reaction  and 
the  formation  of  surface  oxide  complexes  was  established 
as  both  CO  and  CO2  were  desorbed  when  the  target  was  heated 
up  to  high  temperatures  subsequent  to  oxidation.  Judging 
from  the  shape  of  the  etch  pits  observed  after  oxidation, 
it  was  concluded  that  0  and  oxidation  must  involve 
different  mechanisms  and  the  0  oxidation  rate  should  be 
relatively  insensitive  to  the  forms  of  graphite.  These 
conclusions  indirectly  substantiated  the  results  obtained 
by  Rosner  and  Ailendorf  though  the  two  groups  worked  in 

distinctively  different  surface  temperature  regimes. 

64 

Recently,  McCarroll  and  McKee  studied  qualitatively 
the  effects  of  dissociated  oxygen  on  natural  graphite 
crystals  at  temperatures  below  1000°C.  Topographical 
changes  of  the  surface  were  monitored  by  electron  microscopy 
after  the  specimen  was  exposed  to  a  stream  of  dissociated 
oxygen  from  microwave  discharge.  Specifically,  it  was 
found  that  for  surface  temperatures  below  700°C,  no 
reaction  took  place  for  02  whereas  0  produced  unoriented 
pits  over  the  whole  basal  surface.  It  was  concluded  that 
C>2  oxidation  took  place  at  sites  of  lattice  defect  points 
such  as  vacancies,  non-basal  screw  and  edge  dislocations. 

No  such  sites  were  necessary  for  initiation  of  0  oxidation. 

A  common  criticism  for  all  0  oxidation  studies  is  that 


there  is  in  general  no  precise  knowledge  of  the  composition 
and  state  of  the  gas  generated.  Electrode  contamination. 
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carrier  gas  in  discharge  cavity,  ionization  and  electronic 
excitation  of  oxygen  may  all  have  significant  effects  on 
the  surface  phenomena.  Perhaps,  this  is  why  very  few 
quantitative  investigations  have  been  attempted.  There  is 
no  doubt  that  a  better  way  of  generating  oxygen  atoms  would 
aid  the  study  of  the  C-0  reaction  in  a  more  definitive  way. 

C.  Molecular  Beam-Mass  Spectrometry  Technique 

Molecular  beam  techniques  have  been  used  in  different 
areas  of  research  for  a  number  of  years.  Much  of  our 
present  knowledge  of  interatomic  and  intermolecular  inter¬ 
actions  comes  from  atomic  and  molecular  beam  scattering 
studies.  Measurements  of  the  elastic  differential  cross 
section  now  yield  very  accurate  potential  curves  for  systems 

characterized  by  a  single  potential;  Pauly's  work  on 

65 

alkali-rare  gas  systems  is  a  notable  example.  For  non¬ 
reactive  gas-solid  scattering  studies,  sophisticated 
techniques  such  as  time-of-flight  measurements  were  used 

to  characterize  accurately  the  state  of  the  gas  before  and 

6  6 

after  surface  scattering.  Moran's  work  involving  argon 

on  platinum  is  an  example  of  results  that  can  be  used  to 

assess  the  assumptions  made  in  a  scattering  theory  such  as 

6  7 

that  of  Logan  and  Stickney.  In  chemical  research, 
molecular  beam  techniques  have  been  used  for  a  number  of 
years,  particularly  in  the  study  of  gas  phase  reactions, 
surface  ionization  and  chemical  resonance  effects.  However, 
heterogeneous  chemical  kinetic  studies  have  been  made  possible 
only  within  the  past  decade  by  the  development  of  several 
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important  experimental  tools:  high  sensitivity,  compact 

mass  spectrometers;  clean,  high  speed  vacuum  pumps;  and 

phase-sensitive  detectors  or  lock-in  amplifiers.  Recently, 

6  8 

Merrill  reviewed  the  applications  of  molecular  beam 
techniques  to  catalysis  and  reaction  studies  and  predicted 
considerable  sophistication  in  future  experimental  systems. 

Typically,  a  molecular  beam-mass  spectrometry  system 
consists  of  two  vacuum  chambers  separated  by  an  orifice  or 
slit.  A  beam  generating  source  is  housed  in  one  chamber 
while  the  target  and  detector  assembly  is  located  in  the 
other.  The  reactant  gas  from  the  source  interacts  with 
the  surface  as  a  collimated  beam  of  particles  in  parallel 
flight  and  the  products  of  reaction  are  monitored  by  the 
mass  spectrometer.  The  pressure  level  in  the  reaction 
chamber  is  such  that  there  is  virtual  absence  of  gas-gas 

-7 

collisions  (typically  p  ^  10  torr) . 

A  major  advantage  of  this  method  in  studying 
heterogeneous  reaction  is  that  very  reactive  intermediates 
of  the  reaction  can  be  analyzed  as  readily  as  stable  prod¬ 
ucts.  Detection  in  high  vacuum  insures  that  the  observed 
species  represent  primary  products  of  reaction  rather  than 
those  resulting  from  wall  or  gas  phase  collisions  prior  to 
analysis.  Important  experimental  parameters,  such  as 
intensity  of  the  molecular  beam  (the  pressure  of  the  gaseous 
reactant),  temperatures  of  the  beam  and  target  surface, 
incident  beam  and  detector  angles  with  respect  to  the 
surface,  and  composition  of  the  beam  are  all  reasonably 
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controllable.  Again,  due  to  low  pressure  levels,  gas 
diffusion  limited  reactions  are  eliminated  and  the  test 
surface  is  less  prone  to  contamination. 

The  molecular  beam-mass  spectrometry  technique  has, 
however,  one  great  disadvantage:  the  reaction  product 
signal  is  usually  small  and  buried  in  a  relatively  large 
noise  background,  thereby  limiting  the  resulting  signal -to- 
noise  S/N  ratio.  Thus,  any  study  of  chemical  reactions 
using  a  molecular  beam  is  predicated  upon  the  attainment 
of  a  satisfactory  S/N  ratio.  This  problem  can  be  overcome 
in  part  by  use  of  a  "lock-in"  detection  technique.  By 
modulating  or  mechanically  chopping  the  signal  flux  and 
applying  some  form  of  phase-sensitive  (lock-in)  detection, 
the  weak  signal  may  be  extracted  from  the  noise.  An 
analysis  of  this  lock-in  amplification  technique  as  applied 

to  molecular  beams  scattered  from  solid  surfaces  was  given 
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by  Yamamoto  and  Stickney. 

Successful  work  using  this  technique  in  this  same 

29 

laboratory  was  first  done  by  Shih,  who  studied  graphite 
oxidation  by  molecular  oxygen.  Among  efforts  by  other 

groups  was  the  study  of  dissociation  of  hydrogen  on  tantalum 

70  71  72 

by  Krakowski  and  Olander.  Smith  and  Fite  and  McKinley 

investigated  the  nickel-chlorine  reaction.  Germanium 
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oxidation  was  studied  by  Anderson  and  Bondar t  and  the 

dissociation  of  nitrous  oxide  on  tungsten  by  Coltharp  et  al. 

Recently,  an  investigation  pertinent  to  the  present  work 

2  8 

was  completed  by  Olander  et  al,  who  made  a  comparative 


study  of  the  oxidation  mechanisms  of  the  basal  and  prism  | 


planes  of  pyrolytic  graphite. 


In  view  of  the  success  of  the  molecular  beam-mass 
spectrometry  technique  in  studying  heterogeneous  chemical 
reactions/  it  was  decided  that  the  present  investigation 
of  graphite  oxidation  by  dissociated  oxygen  was  feasible. 

D.  Thesis  Objective 

The  object  of  the  thesis  is  to  study  experimentally  the 

reaction  mechanisms  of  dissociated  oxygen  on  graphite  by 
/ 

means  of  molecular  beam-mass  spectrometry  technique. 


Separate  effects  of  atomic  and  molecular  oxygen  on  the 
oxidation  rate  were  to  be  determined  and  isotropic  graphite 
was  to  be  used  as  the  target  material.  In  order  to 
achieve  the  thesis  objective  with  some  insight  into  tne 
dependence  of  the  reaction  on  relevant  parameters/  the 
following  experimental  quantities  as  functions  of  beam  and 
target  temperatures/  equivalent  beam  pressure  and  incident 
angle  were  to  be  measured: 

1.  The  reaction  probabilities  of  atomic  and  molecular 
oxygen . 


2.  The  products  of  reaction  and  their  relative  amounts. 

3.  The  spatial  flux  distribution  of  the  products  and 
reactants  of  reaction  emitted  from  the  target. 

4.  Hysteresis  behavior  patterns. 


A  major  portion  of  the  study  involved  the  design  and 
development  of  a  vacuum  system,  a  high  temperature  0  and 
0^  source,  together  with  a  sensitive  and  stable  detection 
scheme.  To  minimize  design  complications,  all  experimental 
variables  were  restricted  to  the  principal  plane  only,  i.e. 
the  plane  defined  by  the  beam  axis  and  the  normal  to  the 
target. 
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CHAPTER  II 

SYSTEM  DESIGN  AND  SOME  PERTINENT  CONSIDERATIONS 


A.  Generation  of  an  Atomic  Oxygen  Beam 

The  most  common  way  of  generating  atomic  oxygen  pre¬ 
viously  was  by  means  of  an  electrodeless  discharge.  This 

includes  both  radio  frequency  ('v-lO^  Hz)  and  microwave 

g 

(^10  Hz)  discharges.  Apart  from  the  usually  expensive 
and  elaborate  equipment  necessary  for  such  a  system,  there 
are  several  drawbacks.  One  highly  undesirable  feature  is 
that  ionization  and  electronic  excitation  of  the  gas  in¬ 
variably  take  place  along  with  dissociation.  This  may  have 
a  significant  effect  on  the  reaction  process.  Due  to  the 
nature  of  the  discharge,  it  is  impossible  to  control  the 
temperature  of  the  gas  to  cover  a  broad  range  (e.g.,  from 

room  temperature  to  target  surface  temperature) .  It  has 

75 

been  shown  by  Rose  and  Brown  that  for  a  microwave  dis¬ 
charge  the  optimal  breakdown  pressure  for  oxygen  is  about 
4  torr,  in  which  case  recombination  of  atoms  downstream 
of  the  discharge  cavity  is  significant  and  essentially 
reduces  the  atom-to-molecule  ratio  at  the  reaction  surface. 
To  increase  dissociation,  various  carrier  gases  such  as 
water  vapor,  nitrogen,  helium  and  argon  nave  been  added  to 
discharge  cavities.  Such  a  trial-and-error  technique  for  a 
particular  system  depends  mainly  on  the  ingenuity  of  the 
experimentalist.  Though  varying  degree  of  success  has  been 
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achieved  (e.g.,  17%  dissociation  achieved  by  Lam  using 
double  discharge  cavities  with  the  addition  of  water  vapor) , 
the  introduction  of  foreign  gases  in  the  reaction  contributes 
another  factor  of  doubt.  The  same  problems  plague  other  0 
generation  techniques  such  as  high  voltage  DC  discharge  and 
chemical  reaction  (e.g.,  titration  of  N  with  NO). 

Direct  equilibrium  thermal  dissociation  of  oxygen  appears 
to  be  an  attractive  procedure,  since  the  present  design  need 
incorporates  both  the  heating  of  02  as  well  as  the  generation 
of  0.  The  composition  and  state  of  the  gas  is  then  known 
quite  accurately  and  the  risk  of  contaminants  is  minimal. 

An  investigation  into  such  a  scheme  revealed,  however,  some 
associated  problems. 

For  a  symmetrical  diatomic  gas  like  oxygen,  the  dis¬ 
sociation-recombination  reaction  can  be  represented  by 

o2  ^  0  +  0 

The  degree  of  dissociation ,  e  v  can  be  defined  by 

No  _  moNo  _  Mass  of  dissociated  0- atoms 

„o  ..o'  Total  mass  of:  gas 

n  m  n  3 

o 

where  N  =  number  of  dissociated  0-atoms 
o 

N°  *  total  number  of  0-atoms  present  in  the  mixture 

and  m  ~  mass  of  an  0-atom. 
o 
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From  the  law  of  mass  action,  a  at  equilibrium  is  given  by 


i  uraraiA*. 
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2N 


exp(-0d/T) 


(2.1) 


Eq.  (2.1)  can  be  rewritten  in  terms  of  the  equilibrium 
constant,  K(T),  as 


a2  _  K  (T) 
1-g  4p 


or 


a  = 


I 


K  (T) 

K  (T)  +4p 


(2.2) 


where  p  =  total  pressure  of  the  gas  mixture. 

7  8 

From  a  tabulation  of  K(T)  ,  we  can  obtain  the  degree 
of  dissociation  as  functions  of  temperatures  and  pressures. 
Figure  1  shows  such  a  set  of  curves.  For  our  purposes  the 
0-to-C>2  pressure  ratio,  6,  defined  by 


8 


(2.3) 


provides  a  similar  but  more  useful  set  of  curves  shown  in 
Fig.  2.  Clearly,  there  in  practicably  no  dissociation  of 
oxygen  at  temperatures  below  2000°K  and  the  degree  of 
dissociation  is  a  strong  function  of  pressure  at  a  given 
temperature.  For  example,  at  a  source  temperature  of 
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2200CK  -  presumably  a  reasonable  operating  temperature  for 

a  source  fabricated  from  a  refractory .material  in  an 

oxygen  atmosphere  -  6  is  equal  to  11%  at  1  torr  whereas 

when  p  is  decreased  to  0.1  torr,  3  is  about  30%.  Though 

the  ah "olute  flux  of  0-atoms  may  be  higher  for  1  torr 

than  a  corresponding  lower  total  pressure,  the  increase  in 

0-to-02  flux  ratio  is  more  attractive  despite  the  lower 

pressure  level,  since  a  primary  concern  is  a  separating 

out  of  the  effects  of  0  and  0 2  for  the  oxidation  reaction. 

In  principle,  the  problem  of  a  mixture  of  0  and  02  can  be 

eliminated  by  operating  a  source  at  a  very  low  pressure, 

e.g.,  pure  atomic  oxygen  can  be  obtained  at  a  source 
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pressure  of  about  10  torr  and  a  temperature  of  2300°K. 

In  practice,  however,  this  would  mean  a  beam  of  extremely 
low  intensity  and  provide  an  insurmountable  detection 
problem,  at  least,  by  present  day  techniques.  The  related 
problem  of  minimum  detectable  signals  will  be  considered 
in  later  sections.  At  present,  suffice  to  say,  for  a 
maximum  source  operating  temperature  of  2300°K  the  source 
pressure  should  be  less  than  1  torr  for  appreciable 
dissociation  of  oxygen. 

For  such  low  source  pressure  levels,  the  use  of  a 
Knudsen  or  effusive  source  immediately  suggests  itself. 
Though  generally  of  lower  intensity  and  characteristically 
having  a  wider  distribution  of  velocities  than  a  nozzle  or 
Kan trowitz -Grey  source,  the  effusive  source  does  have  some 


advantages  for  present  purposes.  Recombination  of  atoms 
during  the  free  jet  expansion  process  associated  with  a 
nozzle  source  is  thereby  eliminated,  since  there  is  a 
virtual  absence  o^  gas-gas  collisions  for  an  effusive 
source  expansion.  The  properties  of  an  effusive  beam, 
such  as  velocity • distribution,  intensity  and  composition 
can  be  predicted  easily  by  applying  elementary  kinetic 
theory  once  the  source  pressure  and  temperature  are  known. 

In  contrast,  for  a  nozzle  beam,  the  precise  "freezing" 
behavior  is  not  known,  especially  for  the  ill-defined 
nozzles  customarily  employed  in  beam  systems.  Additionally, 
a  skimmer  or  collimator  of  the  nozzle  beam  generally  has  a 
disturbing  effect  on  the  flow  field.  Thus,  in  order  to 
characterize  the  state  of  such  a  beam,  an  experimental 
measurement  such  as  the  time-of-flight  technique  for 
velocity  distribution  is  almost  unavoidable .  Although 
the  total  gas  flow  in  a  nozzle  beam  system  is  orders  of 
magnitude  greater  than  that  for  typical  effusive  beam 
apparatus,  most  of  the  aas  has  to  be  removed  upstream  of 
the  skimmer.  As  it  turns  out,  such  dumping  of  oxygen  into 
the  source  chamber  further  results  in  an  intolerable  background 
oxygen  pressure  that  le~ '  to  rapid  oxidation  of  any  heating 
element.  This  could  be  averted  only  after  major  modifica¬ 
tions  of  the  experimental  setup.  The  reader  is  referred  to 

79 

an  article  by  Anderson  et  al  for  detailed  discussions  of 
the  difference  in  properties  and  operation  of  effusive  and 


nozzle  beams. 
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One  of  the  criteria  of  a  Knudsen  source  is  that  the 
Knudsen  number,  A/D,  where  A  is  the  mean-free-path  of  the 
gas  and  D  is  a  characteristic  dimension  (e.g. ,  the  diameter) 
of  the  orifice,  is  of  an  order  greater  than  unity. 
Necessarily,  the  orifice  dimension  also  should  be  small 
compared  to  the  size  of  the  plenum  chamber.  Physically, 
the  particles  then  will  wander  through  the  orifice  without 
any  collision  during  the  transit  and  without  disturbing  the 
state  of  the  gas  in  the  source.  Under  these  circumstances, 
for  a  circular  orifice  of  negligible  wall  thickness ,  the 
flux  density  or  intensity  of  particles  at  a  distance,  £, 
from  the  orifice  along  a  ray  forming  an  angle,  6,  with 
the  axis  normal  to  the  orifice  area  is 


I  = 


,ns  CS\  ,itD2x  cos  9 
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(2.4) 


As  an  estimate  for  the  appropriate  orifice  size,  consider 
the  case  for  maximum  intensity,  i.e.,  Kn  =  1.  This 
implies 


D  =  A  = 


kT 


/2nm5  /2p  tt6  “ 

»~8 


(2.5) 


For  0^,  6  =  3.0  x  10  cm  and  at  a  sample  source  operating 
condition  of  2200°K  and  1  torr,  D  =  0.058  cm.  This  results 
in  an  0-to-02  intensity  ratio. 
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^  =  0.11, 

°2 

or  about  one  O-atom  for  every  nine  0^ -molecules  in  the 
beam,  a  ratio  which  would  be  too  small  for  separating 
the  individual  effects  of  0  and  C^.  A  one-to-three  ratio 
of  0-to-02  flux  at  a  source  temperature  of  2200°K  requires 
a  source  pressure  of  0.2  torr.  This  means,  however,  a 
five-fold  reduction  in  the  total  beam  intensity.  An 
enlargement  of  the  source  diameter  may  appear  to  solve 
this  particular  problem,  but  it  was  found  during  the  course 
of  the  experiment  that  the  most  troublesome  feature  of  a 
thermal  dissociation  0-at.om  source  proved  to  be  the  rapid 
oxidation  rate  of  the  heating  element  (and  not  the  signal 
level) .  Since  the  total  flux  from  the  source  is  pro¬ 
portional  to  the  square  of  the  source  diameter,  the  bulk 
of  which  must  be  evacuated  by  the  source  chamber  pump,  such 
an  approach  to  enhancement  of  signal  intensity  results  in  a 
drastic  shortening  of  filament  l_re.  Of  course,  it  is 
impossible  to  predict  an  optimal  operating  condition  a  priori 
since  the  optimization  process  involves  many  unknown 
quantities,  only  two  of  which  are  the  oxidation  rate  of 
the  heating  element  and  the  minimum  detectable  signal.  It 
was  only  during  the  course  of  the  experiment  that  an 
ultimate  operating  condition  was  arrived  at  that  satisfied 
all  of  the  practical  requirements. 
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In  practice,  the  orifice  wall  is  of  finite  thickness. 

Equation  (2.4)  for  the  beam  intensity  still  holds  for  the 

beam  axis  (i.e.,  the  6=0  direction),  provided  that  either 

the  wall  thickness,  L,  is  less  than  the  orifice  diameter, 

or  the  Knudsen  number,  A/l,  is  greater  than  unity.  However, 

the  total  flux  out  of  the  orifice  is  reduced;  the  correction 
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factor  was  worked  out  by  Clausing  in  terms  of  the 
dimensionless  variable  L/D. 

The  presence  of  oxygen  at  high  temperature  severely 
limits  the  choice  of  source  and  heater  materials.  Ceramics 
appear  to  be  the  only  choice  for  a  source  chamber  material 
and  thorium  oxide  (thoria)  was  selected  mainly  for  its 
inertness  and  high  melting  temperature  (3300°K) .  Neverthe¬ 
less,  it  does  have  drawbacks  in  connection  with  poor  thermal 
conductivity  and  poor  thermal  shock  resistance.  Extreme 
care  must  be  exercised  in  heating  and  cooling  the  thoria 
to  prevent  cracking,  and  no  significant  load  can  be  applied 
on  the  material  as  it  tends  to  soften  at  temperatures  above 
2300°K.  Unfortunately,  no  extensive  quantitative  data  is 
available  for  the  material,  even  from  the  manufacturer,  and 
knowledge  about  the  physical  properties  must  be  acquired 
from  experience  for  a  particular  application.  Tungsten  was 
chosen  as  the  heater  material,  again  for  its  high  melting 
temperature  (3655°K)  and  its  relative  inertness  in  an  oxygen 
atmosphere.  However,  iti  non-negligible  oxidation  rate  at 
high  temperatures  and  its  reductive  reaction  when  in  contact 
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with  thoria  do  pose  the  problem  of  a  finite  heater  life. 

The  final  design  of  the  source  and  its  performance  will 
be  described  in  Chapter  III. 

It  should  be  noted  that  various  methods  to  increase 
beam  intensity,  such  as  multiparallel  channel  or  slit 
sources,  were  considered  and  rejected  because  of  difficulty 
in  fabrication  and/or  an  implied  intolerable  rise  in 
source  chamber  pressure  level. 

B.  Detection  Scheme 

The  principal  disadvantage  of  molecular  bea;r  ’  lies 
of  heterogeneous  chemical  reactions  is  that  related  to 
working  with  a  small  signal  embedded  in  a  relatively  high 
noise  background.  Basically,  this  stems  from  the  fact  that 
usually  only  a  small  portion  of  the  signal  flux  is  sampled 
directly  by  the  detector.  The  signal-to-noise  ratio 
generally  depends  on  the  mass  of  the  product  measured,  the 
quality  of  the  background  pressure  attainable,  the  incident 
beam  intensity  and,  in  this  instance,  the  reaction  prob¬ 
ability.  By  modulating  the  signal  flux  and  applying  some 

form  of  phase-lock  detection,  the  weak  signal  may  be 

69 

extracted  from  the  noise.  In  recent  years,  sophisticated 
commercial  analog  lock-in  amplifiers  have  become  available 
and  are  almost  indispensable  in  modulated  molecular  beam 
techniques.  The  analog  lock-in  amplifier,  however,  is  not 
without  limitations  and  is  not  necessarily  adequate  for 
this  (Knudsen)  experiment  in  view  of  the  huge  decrease  in 
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incident  beam  intensity  (an  approximate  two  orders  of 
magnitude  drop  from  that  of  a  typical  nozzle  beam)  . 

Theoretically,  any  desired  S/N  ratio  could  be 
obtained  from  a  lock-in  amplifier  output  provided  that  the 
integration  time  is  long  enough;  in  practice,  only  a  finite 
integration  time  is  possible  even  for  the  most  sophisticated 
analog  lock-in  amplifiers ,  mainly  for  reasons  of  capacitor 
leakage  in  the  RC  type  integrator  circuits  used  as  well 
as  circuit  drift  and  instability.  An  integration  time 
constant  of  100  sec.  for  the  Princeton  Applied  Research 
Mcdel  HR-8  lock-in  amplifier  is  typical.  In  addition,  the 
absolute  results  obtained  in  an  analog  detection  scheme  may 
not  be  too  reliable,  since  the  current  or  voltage  type 
output  is  often  not  repeatable  due  to  ^luc  *n..  is  in  both 
electron  multiplier  and  amplifier  gain.'  ic 

calibrations  are  necessary  but  may  prove  too 

cumbersome  in  certain  experimental  setups.  In  view  of 
these  shortcomings ,  a  pulse  counting  type  detection  scheme 
was  investigated. 

For  extremely  weak  signals  (low  in  count  rates) , 
digital  or  pulse  counting  techniques  are  ideal  in  principle. 
By  counting  each  individual  ion  from  a  mass  spectrometer  as 
a  discrete  pulse,  one  essentially  eliminates  the  uncertainty 
factor  in  the  various  gains  of  the  amplifiers.  In  other 
words,  each  ion  is  given  equal  weight  in  computing  the  signal 
level  independent  of  such  electron  multiplier  and  amplifier 
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gain  variations.  By  setting  a  suitable  threshold  level 
in  a  discriminator  (i.e.,  using  a  pulse  height  discrimina¬ 
tion  method),  certain  noise  levels  (e.g.,  white  noises  in 
electron  multiplier  and  amplifiers)  can  be  filtered  out, 
which  is  not  possible  with  an  analog  scheme.  White  noise 
in  an  analog  detection  system  can  be  reduced  only  by 
integration  The  use  of  memory  cells  for  signal  storage 
in  a  pulse  counting  scheme  also  means  indefinite  integration 
time,  limited  only  by  the  capacity  of  the  memory  units. 
However,  the  pulse  counting  technique  has  not  been  fully 
exploited  in  molecular  beam  studies  mainly  because  the 
existing  signal  processing  rate  has  been  inadequate.  The 
highest  count  rates  to  be  handled  with  negligible  counter 
deadtime  are  about  105  counts/sec,  an  ’  '  -»nal  loss  due  to 

pulse-pair  resolution  under  such  con  . -  s  is  significant 
for  most  scattering  experiments.  It  is  due  to  the  recent 
availability  of  a  high  speed,  digital  lock-in  amplifier 
(SSR  Instrument  Co.  Model  1110  Digital  Synchronous  Computer, 
DSC)  together  with  a  compatible  amplifier/discriminator, 

A/D  (SSR  Instrument  Co.  Model  1120)  that  pulse  counting 
becomes  extremely  attractive  for  the  present  application. 
Important  features  of  the  DSC  are:  maximum  count  rate  of 
85  MHz;  pulse-pair  resolution  of  12  nsec;  dual  channels 
for  integration;  and  channel  capac  7  of  10  counts.  The 
A/D,  in  addition  to  its  role  as  a  noise  filter,  also 
generates  standard  size  pulses,  as  required  by  the  DSC, 
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according  to  the  input  pulse  rate.  Due  to  the  extremely 
wide  band  and  high  sensitivity  of  the  A/D,  transients 
from  switches,  running  motors,  signals  from  TV  and  radio 
broadcasts,  and  so  forth  can  all  be  anticipated  to  be 
unavoidable  pickups  by  the  amplifier  and  thus  introduce 
spurious  counts.  Meticulous  care  is,  therefore,  necessary 
to  provide  good  radio  frequency  shielding  and  grounding 
practice . 

Whenever  a  small  signal  in  the  presence  of  a  large 
background  (as  in  the  present  case  to  be  shown  in 
Section  II-C)  has  to  be  retrieved,  signal  modulation 
technique  can  be  applied.  In  this  case,  the  DSC  is 
operated  in  the  so-called  "chopper"  mode.  Trigger  signals 
from  the  modulator  or  chopper  serve  to  initiate  two  counting 
periods  per  chopper  cycle:  one  when  the  chopper  is  open, 
the  other  when  the  chopper  is  closed,  as  shown  in  the  timing 
diagram,  Fig.  3.  The  corresponding  signal  and  noise  pulses 
are  directed  into  two  high  speed  counters,  "A"  and  "B", 
by  means  of  digital  gates.  Pulses  gated  to  counter  "A" 
consist  of  both  the  modulated  signal  and  DC  background  noise 
pulses  whereas  pulses  gated  to  counter  "B"  consist  only  of 
background  noise  pulses.  This  process  continues  until  a 
preset  number  of  chopper  cycles  has  been  achieved.  When 
the  accumulated  pulses  in  "A"  and  "B"  are  subtracted,  the 
resulting  difference  is  that  portion  of  the  signal 
originally  modulated  by  the  chopper.  This  number,  when 
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divided  by  the  number  of  chopper  cycles  and  the  counter 
sampling  time,  is  the  signal  count  rate.  The  counter 
sampling  time  during  each  of  the  two  phases  of  the  operation 
is  generally  set  to  give  a  maximum  duty  cycle,  whereas  the 
selection  of  the  total  number  of  counting  cycles  depends 
mainly  on  the  natural  signal-to-noise  ratio  and  the  accuracy 
of  the  measurement  required.  At  the  completion  of  a 
number  of  preset  chopper  cycles,  the  DSC  can  be  set  to 
either  repeat  the  "same"  measurement  or  hold  the  data. 

For  the  repetitive  operation,  the  set  of  data  is  available 
until  the  end  of  the  next  counting  period. 

To  estimate  the  total  counting  time  required  to 
achieve  a  given  accuracy  of  measurement,  an  analysis  will 
be  given  assuming  the  errors  introduced  by  the  processes 
of  totaling  the  counts  and  measuring  the  elapsed  time  are 
negligible  compared  with  the  statistical  error.  This 
assumption  is  true  for  most  counting  equipment.  Further¬ 
more,  all  fluctuations  of  both  the  signal  and  noise  level 
are  assumed  to  follow  the  Poisson  distribution  (i.e.,  all 
events  are  random  and  uncorrelated) .  The  standard 
deviation,  o,  will  be  used  for  the  precision  index  in  **he 
following  discussion. 

If  N  is  the  number  of  counts  summed  over  a  ti. 
interval,  t,  the  counting  rate  R  is 

R  =  |  (2.6) 
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This  value  with  its  standard  deviation  may  be  stated  as 
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R  *  aR 


R  ±  (|)1/2 


(2.7) 


The  error  is  then 

Op  -l 

4  =  - i-T72  =  -T7J  (2.8) 

R  (R  t)WZ  Ni/2 

Eq.  (2.8)  is  a  formal  statement  of  the  very  important 
conclusion  that  the  percentage  error  of  a  counting  measure¬ 
ment  is  determined  entirely  by  the  total  number  of  counts 
accumulated. 

When  several  quantities  having  standard  deviations 
o^r  o0 i • • • »  on  are  combined  by  either  addition  or  subtrac¬ 
tion,  the  standard  deviation  of  the  result  is 

cj  —  ( t  ...  t  0^  )  (2.9) 

For  the  present  case  of  a  background  or  noise  counting 
rate  of  R^j  ±  oN  and  a  total  counting  rate,  due  to  both 
signal  and  the  background,  of  R^  ±  oT,  the  counting  rate, 
Rs,  due  to  the  signal  only  is 

Rs  ±  cs  =  (RT  -  Rjj)  ±  (oT2  +  cjn2)1/2  (2.10) 

The  error  for  the  signal  counting  rate  is  then 
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Since  the  sampling  times  in  channels  "A"  and  "B"  are  the 
same,  Eq,  (2.11)  becomes 
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(2.12) 


Eq.  (2.12)  is  essentially  the  .eciprocal  of  the 
s.i gnal-to-noise  ratio  after  integration,  namely. 


int 


R  t1'2 
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<2En  *  v1' 


(2.13) 


Rq/^M  signal— to-noise  ratio  of  the  system 

before  the  integration  scheme  and  its  estimation  is  given 
in  Section  II-C.  With  typical  values  of 


Rg  =  10  cps  and  Rg/^  =  1/5000 

an  integration  time  of  10^  sec.  (approximately  3  hours)  is 
required  for  a  maximum  error  of  10%  in  the  measured  signal. 
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The  total  number  of  counts  accumulated  in  this  case  is 
9 

5  x  10  ,  vh.v'.h  will  have  exceeded  the  existing  counter 
capacity.  To  alleviate  this  problem,  the  integration 
scheme  may  simply  be  carried  out  in  a  series  of  shorter 
integration  periods  within  which  the  total  pulse  count 
is  less  than  the  saturation  limit  of  the  counters.  This 
burden  of  recording  many  sets  of  data  has  the  advantage, 
however,  of  providing  a  means  to  eliminate  some  of  the 
spurious  counts  without  jeopardizing  the  whole  measure¬ 
ment  as  will  be  discussed  in  Section  II-E. 

The  pulse-pair  resolution  limit  of  the  counting 

equipment  is  a  measurement  of  the  counter  dead  time.  A 
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first  order  correction  for  the  dead  time  is  given  by 


R  =  Ro/(l  -  ;io  t) 


(2.14) 


where  R  is  the  true  count  rate 

R  is  the  observed  count  rate 

O 

and  t  is  the  pulse-pair  resolution  time. 

_  8 

For  a  value  of  t  of  10  sec,  the  correction  is  approxi- 

5 

mately  0.1%  at  a  count  rate  of  10  pulses  per  second. 


which  can  be  safely  neglected. 

C.  Signal  Output  and  Signal-to-Noise  Considerations 


A  flow  chart  for  the  signal  processing  scheme  is  shown 


in  Fig.  4.  The  signal  flux  from  the  target  or  the  incident 
beam  was  first  modulated  by  the  tuning  fork  chopper  before 
its  mass  spectrum  was  analyzed  by  the  mass  spectrometer. 
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The  signal,  in  the  form  of  molecules  or  atoms,  is  con¬ 
verted  into  ions  in  the  mass  spectrometer  ionizer.  The 
resulting  ion  number  flux  is  given  by 


F. 

ion 


n.  a. 
ion  ion 


a  i 
e  e 


(2.15) 


since  the  mass  spectrometer  ionizer  is  a  density  sensitive 
detector.  In  Eq.  (2.15),  it  is  assumed  that  the  particles 
are  singly  ionized  and  that  the  signal  flux  is  constant 
across  the  length  of  the  ionizer  so  that  the  local, 
instantaneous  number  density  in  the  ionizer  is  uniform. 

Only  a  portion  of  the  ions  produced  in  the  ionizer 
will  pass  through  the  spectrometer  and  arrive  at  the  first 
dynode  surface  of  the  electron  multiplier,  their  fraction 
depending  upon  the  ion  drawout  efficiency  of  the  mass 
spectrometer  and  the  mass  number  to  which  the  mass  spectro¬ 
meter  is  tuned.  This  fraction,  essentially  the  mass 
spectrometer  transmission  factor,  t  ,  depends  on  the 
design  of  the  mass  spectrometer  and  is  generally  quite 

small.  For  a  double-slit,  magnetic  deflection  type  mass 
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spectrometer,  a  value  of  10  is  typical. 

The  total  gain,  Gem,  of  an  electron  multiplier  used 
as  an  ion  detector  is  given  by 
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Gem  =  gl  g2 


(2.16) 


where  g^  is  the  first  dynode  conversion  efficiency  in 
electrons/ion 

g ^  is  the  secondary  emission  coefficient  of  the 
remaining  dynodes 

and  a  is  the  number  of  remaining  dynode  stages, 
g^  and  g ^  will/  of  course,  vary  with  incident  ion  energy, 
dynode  surface  material  and  electron  multiplier  voltage. 
The  electron  current  output  of  the  mass  spectrometer, 
imp ,  is  then  given  by 


q  F.  t  G 
M  ion  ms  em 


n.  a.  Hit  g,  q  * 
ion  ion  e  e  ms  ^1  y2 


(2.17) 


The  performance  of  the  mass  spectrometer,  usually  expressed 

in  terms  of  a  sensitivity,  ,  can  now  be  evaluated.  The 

ms 

mass  spectrometer  sensitivity,  in  torr  ^ ,  is  given  by 


i  0 .  H  t  q. 

ms  _  ion  e  ms  yl  *2 

ie  p  zr 


(2.18) 


where  p  =  n.  kT  and  T  =  300°K 
r  ion  a  a 

In  a  pulse  counting  scheme,  the  absolute  value  of 
the  electron  multiplier  gain  is  unimportant  as  compared  to 
the  overall  ion  detection  efficiency,  Eem,  in  detected  ions/ 
incident  ion.  Thus,  the  pulse  frequency  output  of  the  mass 


#*v-J  ’-,  ^  •  fv  v  - 


spectrometer,  F  ,  can  be  written  as 


n .  c .  £  i  t  E 

F  =  F.  t  E  =  -jSSL-iffl—g.  e.  JBL- *§  (2.19) 

ms  xon  ms  em  q 


For  incident  ion  energies  between  2.5  and  5  KeV,  the 
efficiency  remains  substantially  constant  at  a  value  of 
about  unity.  As  before,  we  can  define  another  mass 
spectrometer  sensitivity,  S^g,  in  terms  of  the  pulse 
frequency  output: 


.1  _  ms 

’ms  ~  ~ I  p 
e  * 


a  ■  £  t  E 

xon  e  ms  em 


(2.20) 


The  units  for  are  pulses/seo-anp-torr. 


From  Eqs.  (2.18)  and  (2.20), 


c  /  em  ) 
ms  q  G* 
u  em 


(2.21) 


This  result  can  be  used  to  estimate  the  sensitivity  of 

the  present  pulse  counting  scheme  from  results  obtained  in 
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the  analog  system  used  by  Shih.  Choosing  typical  values 
as  follows. 


S  =  300  (torr) 
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results  in. 


1.87  x  1015 


pulses 

sec-amp-torr 


(2.22) 


A  convenient  expression  for  estimating  the  signal  output 
pu. se  frequency  will  be  one  that  relates  the  number  of 
ion  output  pulses  to  the  number  density  at  the  ionizer 
as  given  by  Eq.  (2.19).  In  terms  of  ,  this  can  be 
written  as 


F 


ms 


n. 

ion 


n.  S* 
ion 


(2.23) 


S*  is  simply  a  modified  mass  spectrometer  sensitivity  that 

converts  the  number  density  at  the  ionizer  directly  to 

output  pulse  frequency.  Its  magnitude,  like  that  of  S*s, 

depends  on  the  detected  gas  species.  However,  for  the 

gases  of  interest  in  this  experiment,  the  variation  in 
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magnitude  of  S*  will  be  less  than  25  per  cent.  At  an 
electron  emission  current  of  1  ma  and  a  mass  spectrometer 

sensitivity  as  given  in  Eq.  (2.22),  a  number  density  of 

3  -2 

1  part/cm  at  the  ionizer  will  give  rise  to  5.8  x  10 

pulses/sec  at  the  output  of  the  electron  multiplier. 

An  estimation  of  the  signal-to-noise  ratio  will  now 

proceed. 

The  incident  beam  intensity  at  the  target  surface, 
assuming  that  the  target  is  in  the  beam  axis  and  that  it 
subtends  a  small  solid  angle  as  viewed  from  the  source 
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orifice,  is  given  by  Eq.  (2.4): 


IT 


(~r~)  (-V) 


(2.24) 


total  particle  flux  on  the  target,  Ftf,  is 


ft  =  IT  A^,  cos  ei 


(2.25) 


where  6  ^  is  the  beam  incident  angle  on  the  target  relative 
to  the  surface  normal.  Assuming  all  products  of  reactions 
are  diffusively  emitted  from  the  target,  and  that  the 
target  can  be  treated  as  a  point  source  as  viewed  by  the 
detector,  then  the  total  reaction  product  flux  of  the  jth 
species  entering  the  detector,  i.e.,  the  mass  spectrometer 
ionizer,  is 


(FjJo  =  Ft 


p,  cos  eD  AD 


(2.26) 


where  is  the  probability  of  the  formation  and  desorption 
of  a  jth  reaction  product  species  for  every  incident 
particle.  (A  consideration  of  the  correction  of  the  particle 
flux  due  to  the  actual  target-detector  geometry  is  given  in 
Appendix  A.)  P^  is  related  to  the  reaction  probability,  e  , 
defined  as  the  ratio  of  the  carbon  atom  removal  rate  to  the 


incident  particle  flux  by, 


«* xt*-**  Yi: 


-  33P..' 


Carbon  atom  removal  rate 
Incident  particle  flux 


ZCj  pj  ft 


J  Ci  Pj 


(2.27) 


where  Cj  is  the  number  of  carbon  atoms  in  the  jth  reaction 
product. 

Now,  the  total  jth  reaction  product  flux  entering 
into  the  ionizer  can  be  written  as 


<Fj ) D  (nj*D  Vj  AD 


(2.28) 


where  Vj  is  the  average  velocity  of  the  jth  reaction 
product.  If  the  products  of  reaction  are  in  thermal 
equilibrium  with  the  target  surface  and  they  desorb  from 
the  surface  with  a  Maxwellian  velocity  distribution,  then 


l^kTr, 


V.  = 

3 


(2.29) 


and  therefore. 


(n.)_  =  (F . )  •  _ 1 

3  D  3  D  J  9TrfcT 


(2.30) 


Eq.  (2.30)  represents  the  signal  level  of  the  jth  reaction 
product  that  would  be  detected  by  the  mass  spectrometer. 

An  actual  pulse  frequency  output  from  the  electron 
multiplier  can  be  obtained  from  Eq.  (2.23) . 


il.  3 
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In  a  molecular  beam  experiment,  the  signal  is 
usually  buried  in  a  high  background  noise  level.  For 
a  pulse  counting  scheme,  the  noise  that  needs  to  be 
considered  is  that  corresponding  to  ambient  gases  in  the 
test  chamber  which  have  the  same  mass  number  as  the 
signal  being  detected.  Electronic  noises,  such  as 
electron  multiplier  dark  current  and  input  resistor  shot 
noise,  can  all  be  eliminated  by  using  the  previously 
mentioned  pulse  height  discrimination  technique.  The 
ambient  gases  are  due  to  two  sources  primarily.  One 
is  the  residual  gas  mixture  in  the  vacuum  system,  the 
amount  of  which  depends  on  the  system  design  and  the 
pumping  speed.  For  an  oil  diffusion  pumped  and  liquid 
nitrogen  cooled  vacuum  system,  the  most  common  background 
gases  are  hydrogen,  helium,  water  vapor,  carbon  monoxide 
and  nitrogen.  The  other  source  for  ambient  gases  is  that 
of  products  of  reaction  from  the  target  surface. 

The  total  reaction  product  flux  of  the  jth  species 
leaving  the  target  surface  is 

Fj  =  Ft  P.  (2.31) 

This  will  result  in  an  ambient  partial  pressure  of  the 
jth  gas  component  in  the  test  chamber, 


F  P . 


(2.32) 
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where  is  the  test  chamber  pimping  speed.  If  we  denote 
the  partial  pressure  of  the  jth  gas  component  due  to  the 
residual  gases  as  (p.)_,,  then  the  total  ambient  partial 
pressure  of  the  jth  gas  component  is 

<Pj>a-  (Pj>T+  (Pj>R=^i+  <P:»R  (2-33» 

In  terms  of  number  density ,  equivalently 


a 


(Pj)R] 


(2.34) 


where  T  is  usually  equal  to  the  room  temperature.  The 

cl 

signal-to-noise,  S/N,  ratio  for  the  jth  reaction  product 
is  then. 


S  ln-i )  n 

<l>j  -  (2-35) 

J  a 

where  (n.)D  and  (n.)  are  given  by  Eqs.  (2.30)  and  (2.34) 
3  u  3  a 

respectively. 

For  an  estimation  of  the  magnitude  of  the  S/N  ratio, 
the  following  realistic  parameters  may  be  used: 


Jr 


F'.vw  >- ~y+>  -  %*r ' 
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1  torr. 

T  =  2200  °K 
s 

0.025  cm, 

iST  =  9  cm 

.  2 

1  cm  , 

A^  =  1  cm2 

—  4  S  0 

D  4-  ' 

HTD  =  18  cm 

1300°K , 

T  =  300°K 
a 

-2 

10 

(Pj>R  ■  10" 

T 

V  9D 

CT  =  1 

?.  =  1 
3 


V^  =  2000  litres/sec 


Assuming  CO  is  the  reaction  product,  the  resulting  S/N 
ratio  is 


(§.)  -  x.1.0  « 

n)c0_  3,22  x  107  "  2700 


(2.36) 


This  is  a  very  small  S/N  ratio  indeed,  and  without  any 
S/N  enhancement  scheme  the  signal  would  be  impossible  to 
detect.  For  C02  detection,  the  S/N  ratio  is  much  improved 
since  there  is  usually  little  or  no  residual  CC>2  gas 
present,  and  thus  the  noise  is  contributed  entirely  by 
the  reaction  product  flux  from  the  target  surface.  In 
this  case,  assuming  (Pc0  )R  to  be  zero  and  using  similar 
parameters,  the  resulting  S/N  ratio  is 


{-) 

Vco. 


=  1-2]  x  10  a  1_ 
9.50  x  105  80 


(2.37) 
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It  is  worthwhile  to  point  out  the  advantage  of  having 
a  huge  pumping  speed,  as  illustrated  by  Eq.  (2.34),  using 
the  above  parameters  for  CO  detection: 


(n.) 


j  a  kT 


F  P  . 

~  [~f— 1  +  (Pj)^ 


a  V, 


or 


-9 

(nCO)  a  =  k?“  {2*95  x  1'3~11+10"9)==  ^—(0.0295+1)  (2.38) 

cl  SL 


We  can  see  that  the  ambient  noise  is  contributed  almost 
entirely  by  the  residual  gases  with  only  about  3%  being 
due  to  the  reaction  product  flux  in  this  case.  The 
advantage  of  a  high  pumping  speed  is  even  more  dramatically 
demonstrated  for  the  case  of  C02  detection,  as  shown  by 
the  relatively  high  S/N  ratio. 

In  estimating  the  signal  output  level,  it  was 
assumed  that  the  particle  beam  is  well  collimated  so  that 
all  particles  effusing  from  the  source  orifice,  except  for 
the  portion  directed  at  the  target  surface,  are  pumped 
away  in  the  source  chamber.  This,  in  fact,  was  the  case 
for  most  of  the  present  experimental  runs.  It  was  found, 
however,  that  for  source  temperatures  above  2200 °K,  an 
efficient  way  of  prolonging  the  heater  life  was  to  move 
the  source  orifice  close  to  the  source  collimator  so  that 
ssentially  all  particles  actually  are  effused  into  the 
.>■  ;  chamber.  The  drawback  of  this  solution  to  the  heater 


t 
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life  problem  is  that  the  resulting  increase  in  oxygen 
pressure  in  the  test  chamber  may  contribute  a  significant 
portion  to  the  detected  signal. 

The  total  flux  out  of  an  effusive  channel  source  is 


(2.39) 


where  K  is  the  Clausing  factor.80  The  resulting  increase 
in  ambient  pressure  in  the  test  chamber  is 

<F  >.  =  ^  (2.40) 

s  a 

t 

and  the  incident  particle  flux  on  the  target  due  to  this 
pressure  is  given  by 


(ft) 


a 


(ps>a 

727rmkTa 


*  *T 


(2.41) 


A  ratio  of  this  incident  flux  and  the  flux  directly  from 
the  source  orifice  onto  the  target  (Eq.  (2.25))  furnishes 
an  estimate  of  the  correction  that  is  necessary  to  arrive 
at  the  reaction  probability.  This  ratio,  after  simplifica¬ 
tion,  is 

^FT^  a  _  K7tX,ST 

FT  V.  /2umkT  cos  0 . 


(2.42) 
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With  the  same  parameters  previously  used  for  estimating 
the  S/N  ratio  and  the  additional  parameter  of  L  =  0.080  cm, 
Eq.  (2.42)  results  in 

<Va 

— — —  =  0.15  (2.43) 

FT 

a  correction  which  is  not  negligible. 

D.  Data  Reduction  i Scheme 

The  reaction  probability,  as  defined  by  Eq.  (2.27) , 
is  obtained  experimentally  by  measuring  the  incident 
beam  flux  and  the  reaction  product  angular  flux  distribu¬ 
tion  from  the  target  surface  in  the  principal  plane, 
i.e. ,  the  plane  defined  by  the  beam  axis  and  the  normal 
to  the  target.  The  reaction  product  angular  flux 
uistribution  is  then  assumed  to  be  axisymmetric  with 
respect  to  the  surface  normal.  This  assumption  appears 
to  be  reasonable  in  view  of  the  actual  measured  angular 
flux  distribution,  which  shows  that  the  reaction  product 
desorbs  in  a  lobular  fashion  towards  the  surface  normal. 

The  incident  beam  flux  is  measured  by  positioning  the 
mass  spectrometer  ionizer  directly  on  the  beam  axis.  In 
this  case,  the  beam  intensity  at  the  detector  is 


(2.44) 


i/V'fo  yA  !•*■ 


**?*»v*t 


Thus,  the  beam  intensity  at  the  target  is  obtained  by 


5.  2 

*T  '  ‘Vi  <I^> 


(2.45) 


For  a  Maxwellian  beam,  the  average  beam  velocity  is 


9  tpcT, 


VB  = 


(2.46) 


Therefore,  the  number  density  at  tlie  detector  for  the 
incident  beam  is 


ni  ^D^i  j9  7rkT. 


(2.47) 


The  pulse  frequency  output  from  the  mass  spectrometer 
for  the  incident  beam  is  then 


(F  ) •  =  n.  S. 
ms  i  li 


(2.48) 


where  denotes  the  mass  spectrometer  conversion  factor 
for  the  incident  gas.  From  Eq.  (2.25),  the  total  particle 
flux  on  the  target,  FT,  in  terms  of  the  detected  pulse 
frequency  of  the  incident  beam  is 


A  cos  0  /9TikT  8,  2 

F  =  (F  ).  — - i.  / - - 

T  ms  i  s*  J  &ns  St 


(2.49) 


Consider  the  case  when  the  products  of  reaction  are 
diffusively  desorbed  from  the  target  surface  and  let  the 


intensity  of  the  jth  species  measured  at  the  normal 
direction  (i.e.,  @D  =  0)  be  then  the  total  flux 

of  the  jth  reaction  product  desirbed  from  the  surface, 
(F^,  is 


2tt  tt/2 

/  (I.n)  .  cos  6  dm  =  /  /  (I1n)  .  cos  0  sin  0 

fl  1U  3  «f>=0  0=0  3 

Ti(I10)j  (2.50) 


If  (FjHg)  j  denotes  the  measured  pulse  frequency  output  from 
the  mass  spectrometer  for  the  jth  reaction  product  at 
0Q  =  O,  then  using  Eqs.  (2.30)  and  (2.50)/  (F^) .  is  given 
by 


(F  ), 
ms  j 


(2.51) 


Therefore,  the  reaction  probability  is 


reaction  products  desorb  in  a  lobular  manner  in  the  normal 

direction  and  the  flux  distribution  can  be  described 

empirically  by  cosn0,  then  the  total  flux  of  the  jth 

reaction  product  desorbed  from  the  surface,  (F  ).,  is 

n  3 


<Vj  ’  /  (InO)j  cosne  a“ 
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(2.53) 


where  (1^)^  is  the  measured  intercity  of  jth  species  at 
the  normal  direction. 

The  reaction  probability,  in  this  case,  is  simply 
obtained  by  multiplying  Eq.  (2.52)  by  the  factor  2/(n+l). 

For  a  dissociated  incident  beam,  the  reaction 
probabilities  of  0  and  0 2  have  been  determined  by  applying 
the  method  of  superposition.  In  principle,  this  is  only 
applicable  for  linear  surface  processes.  Although  non¬ 
linear  surface  processes  are  quite  common  (e.g.,  bimolecular 
recombination) ,  the  method  is  justified  in  this  experiment 
in  view  of  the  extremely  low  incident  flux  on  the  surface 
and  the  observed  first  order  pressure  dependence  of  the 
oxidation  rate.  The  reaction  probability  of  0^  at  the  dis¬ 
sociated  beam  temperature  was  obtained  by  extrapolation  from 
results  of  lower  beam  temperatures,  at  which  the  beam  con¬ 


sists,  essentially,  of  pure  C^. 
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Up  to  this  point#  the  analysis  has  not  considered 
the  molecular  transit  effects.  Due  to  the  finite 
molecular  velocities,  there  will  be  a  time  lag  between 
the  time  when  the  molecules  first  pass  the  chopper  and 
when  they  are  detected  at  the  mass  spectrometer  ionizer. 

This  time  lag  is  usually  expressed  as  a  phase  lag  in  terms 
of  the  chopper  frequency.  Previously,  it  was  assumed 
that  the  molecules  travel  with  an  average  velocity  of  the 
molecular  beam,  making  it  quite  straight  forward  in  de¬ 
termining  the  phase  lag  for  detection  process.  Moreover, 
the  wave  shape  of  the  chopped  beam  will  be  exactly  similar 
to  the  chopper  opening  shape.  However,  since  the  beam 
contains  the  whole  Maxwellian  velocity  spectrum,  the 
chopped  wave  shape  will  be  distorted  at  the  detector, 
depending  on  the  beam  temperature  and  flight  path.  In 
other  words,  the  chopped  wave  shape  will  vary  as  the 
faster  molecules  move  ahead  and  the  slower  molecules  fall 
behind  from  the  average  speed  molecules.  In  general,  this 
results  in  a  decrease  in  the  wave  amplitude  and  the  effect 
is  known  as  dispersion  attenuation. 

The  dispersion  attentuation  and  p  se  shift  for 
Maxwellian  molecular  beams  have  been  computed  and  tabulated 
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for  different  beam  properties  by  Harrison,  Hummer  and  Fite.  ' 
The  magnitude  of  these  corrections  depends  on  the  chopper 
frequency,  fc,  the  chopper-to-detector  distance,  S-CD,  and 
the  gas  temperature,  T,  as  reflected  in  the  most  probable 
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1/2 

speed,  V  =  (2kT/m)  .  The  amplitude  dispersion 

attenuation  and  phase  shift  are  tabulated  as  functions 
of  the  dimensionless  parameter,  X,  defined  by 


X  = 


(2.54) 


For  f  =  377  Hz,  &CD  =  6  cm,  and  a  molecular  oxygen  beam 
of  300°K,  X  =  0.363  giving  an  amplitude  dispersion 
attenuation  of  4%  and  a  phase  lag  of  22°,  as  compared  to  a 
16°  lag  from  the  simple-minded  calculation  of  using  the 
average  beam  speed,  V,  in 


* 


(2.55) 


4>  being  the  phase  lag  in  degrees. 

For  gas  at  higher  temperatures ,  the  amplitude 
attenuation  is  even  less  and  the  difference  between  the 
L base  lag  from  the  exact  analysis  and  Eq.  (2.55)  will 
become  almost  negligible.  Curves  of  the  amplitude  attenua¬ 
tion  and  phase  lag  for  range  of  X  of  interest  in  this 
experiment  are  available  in  Ref.  29  and  they  are  reproduced 
in  Fig.  5.  It  can  be  seen  that  for  minimal  corrections, 
short  flight  paths  are  applicable  for  a  given  chopper 
frequency.  However,  long  flight  paths  are  desirable  when 

the  transit  dispersion  effects  are  measured  for  determining 
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the  speed  spectrum  of  a  molecular  beam. 
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E.  Statistical  Treatment  of  Experimental  Data 

In  this  section  the  rather  controversial  question  of 
data  rejection  will  be  considered.  It  concerns  the 
problem  of  what  to  do  if,  among  a  set  of  observations, 
one  or  more  have  deviations  so  large  as  to  seem  unreasonable. 
Such  an  observ  on  creates  an  awkward  situation  for  the 
experimenter.  If  one  retains  the  questionable  observation, 
it  can  have  quite  a  large  effect  on  the  mean.  If,  on  the 
other  hand,  it  is  discarded,  one  runs  the  risk  of  throwing 
away  information  which  might  lead  to  discovery  of  some 
unexpected  phenomenon  in  the  experiment.  Important  dis¬ 
coveries  have  resulted  from  apparently  anomalous  data. 

Or.e  could  take  the  point  of  view  that  there  is  never  any 
justification  for  discarding  data  on  purely  statistical 
grounds  unless  there  is  a  definite  reason  for  suspecting 
that  a  particular  observation  is  invalid.  In  that  case, 
there  is  nothing  more  to  say  except  to  advocate  taking 
sufficient  data  such  that  the  results  are  "unaffected"  by 
the  questionable  observations. 

Due  to  the  extremely  wide  band  and  sensitive  amplifier 
used  in  the  detection  scheme,  strong  transient  electronic 
noises  from  the  surroundings  are  liable  to  be  picked  up 
as  signals  by  the  counter,  even  though  the  input  end  of 
the  amplifier  has  been  subjected  to  meticulous  shielding. 
These  transient  noises  (e.g.,  arcing  caused  by  bad  contacts 
in  electrical  switches)  are  usually  high  frequency  bursts 
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with  very  short  durations  that  cannot  be  averaged  out 
by  the  integration  scheme.  Due  to  the  low  signal 
frequency  output/  taking  additional  data  to  minimize  the 
effects  of  these  noises  will  require  an  unreasonable 
observation  time.  It  is  precisely  with  this  problem  in 
mind  that  a  systematic  treatment  of  the  experimental  data 
was  considered. 
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Chauvenet's  criterion  is  frequently  used  for  the 
purpose  of  rejecting  outlying  data.  The  criterion  states 
that  any  one  in  a  series  of  N  readings  will  be  rejected 
when  its  deviation  from  the  mean  of  the  series  is  such  that 
the  probability  of  the  occurrence  of  such  a  deviation  being 
as  large  or  larger  is  less  than  1/2N. 

In  general,  the  probability  for  a  measurement  to  occur 
in  an  interval  within  ko  of  the  mean  is 


P  (k)  =  / 

V2tt  -k 


k  -xV2 


(2.56) 


R 


1 

N 


(2.58) 


£q.  (2.56)  is  the  integral  of  the  Gauss  or  normal  error 

function  and  numerical  values  of  this  integral  are 
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extensively  tabulated. 

The  probability  for  a  measurement  to  occur  in  an 
interval  larger  than  ko  of  the  mean  is  [1-P(k)3  and 
according  to  Chauvenet's  criterion,  a  reading  vrll  be 
rejected  if 

l-P(k)  <  ijj-  (2.59) 

For  example,  if  there  are  ten  readings  for  one  experi¬ 
mental  measurement,  i.e.,  N=10,  then 

P (k)  ^  1-0.05  =  0.95 
and  from  Ref.  37,  k  1.96 

This  means  that  a  reading  should  be  rejected  if  its 
deviation  from  the  mean,  |R^-Rj,  is  greater  than  1.96 
times  the  standard  deviation,  a.  Table  II-l  gives  the 
magnitude  of  this  .'aviation  in  terms  of  multiples  k  of  a 
for  several  values  of  N. 


Table  II-l 

Maximum  Acceptable  Deviations  in  Accordance 
with  Chaurvnet’s  Criterion 


Number  of  readings,  N 

2 

5 

7  10  15  20  25 

Ratio  of  deviation  to 
standard  deviation,  k 

1.15 

1.65 

1.80  1.96  2.13  2.24  2.33 
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This  method  of  "filtering"  the  experimental  data  has 
been  adopted  in  the  present  work.  After  the  "anomalous" 
observations  are  eliminated  in  accordance  with  Chauvenet’s 
criterion,  the  mean  and  the  standard  deviation  are  re¬ 
computed  using  the  remaining  observations.  The  criterion 
for  discarding  experimental  data  is  carried  out  once  more 
and,  if  necessary,  the  whole  process  is  repeated  until  all 
observations  are  within  the  maximum  acceptable  deviation. 
For  widely  scattered  data,  the  repeated  applications  of 
Chauvenet's  criterion  may  lead  to  the  elimination  of  most 
or  all  of  the  data.  Fortunately,  in  the  present  work,  it 
has  not  been  necessary  to  carry  out  the  procedure  for  more 
than  once  and  in  no  case  the  elimination  of  more  than  two 
observations  is  required. 


Preceding  page  blank 
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CHAPTER  III 

EXPERIMENTAL  DETAILS 

A.  Description  of  Apparatus 
1.  General  Features 

29 

The  present  experiment  and  the  work  of  Shih 

represent  the  first  generation  of  molecular  beam  studies 

of  heterogeneous  chemical  reactions  in  the  Molecular  Beam 

Laboratory  of  the  M.I.T.  Department  of  Aeronautics  and 

Astronautics.  Although  the  ideal  design  of  an  experimental 

setup  involves  the  capability  of  controlling  a  maximum 

number  of  experimental  parameters  for  studying  heterogeneous 

reactions,  the  complexities  introduced  with  each  variable 

generally  make  it  impractical  to  do  so.  An  early  decision 

was  made  to  emphasize  a  system  with  a  high  detection 

sensitivity  (i.e. ,  capable  of  studying  heterogeneous 

-3  -4 

otems  with  reaction  probabilities  as  low  as  10  to  10  ) 

together  with  reasonable  flexibility  (i.e.  ,  parameters  such 
as  incident  beam  angle,  reaction  product  detector  angle, 
detector  angular  resolution,  to  be  readily  varied) .  No 
velocity  analyzing  scheme,  such  as  mechanical  velocity 
selection  or  time-of-f light  technique,  was  incorporated  into 
the  design  because  such  a  scheme  usually  implies  a  severe 
reduction  in  signal  levels.  In  view  of  the  high  operating 
temperature  range  of  the  source  and  target,  key  components 
of  the  system  were  all  designed  to  be  water-cooled  and 
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materials  in  the  vacuum  chambers  were  chosen  to  have 
minimum  outgassing  rates  wherever  possible.  New  parts, 
such  as  braided  fibar  glass  insulators,  were  always 
degassed  in  a  separate  oven  at  temperatures  of  about 
400 °C  before  they  were  placed  in  the  vacuum  system, 
whereas  metal  components  were  subjected  to  repeated 
cleansing  in  solutions  of  trichloroethylene  and  acetone. 
An  ultrasonic  cleaner  was  also  used  to  clean  parts  which 
could  not  be  done  by  hand  easily.  Fabrication  of  the 
vacuum  chambers  made  use  of  heliarc  welding  to  exclude 
atmospheric  gases  and  prevent  the  formation  of  nitrides 
and  oxides  on  or  in  welds.  Interior  surfaces  as  well  as 
all  internal  components  of  the  vacuum  chambers  were 
mechanically  polished  whereas  the  exterior  surfaces  had 
sand-blast  finishes.  An  overall  view  of  the  experimental 
setup  is  shown  in  Fig.  6. 

2.  Vacuum  System 

The  vacuum  system  consisted  of  two  separately 
pumped  chambers  —  an  18"  diameter  by  25"  test  chamber 
housing  the  target-detector  assembly  and  a  12 ‘  diameter 
by  20"  source  chamber  which  housed  the  beam  generating 
assembly.  The  two  chambers  were  connected  by  a  12" 
diameter  flanged  opening,  as  shown  in  Fig.  7.  Both  were 
fabricated  from  stainless  steel  SS  304. 

The  test  chamber,  in  the  shape  of  a  bell  jar,  had  a 
12"  diameter  flanged  opening  in  its  side,  opposite  to  the 
flanged  connection  to  the  source  chamber,  for  easy  access 
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to  the  interior  of  the  chamber  without  necessitating 
lifting  up  of  the  whole  chamber.  It  was  evacuated  by  a 
National  Research  Corporation  10,000  liters/sec,  16"  oil 
diffusion  pump  (Model  HS16)  using  Silicone  D.C.  705  oil. 

Due  to  its  relatively  high  conductance,  a  copper,  chevron- 
type  baffle,  cooled  by  either  liquid  nitrogen  or  Freon  12 
refrigerant,  was  used  to  prevent  backstreaming  of  the 
pump  oil  vapor  and  to  trap  condensible  gases.  On  top  of 
the  baffle  was  a  4"  high,  aluminum  feedthrough  ring  with 
eight  equally  spaced  3"  diameter  ports.  One  of  the  ports 
was  connected  to  a  mechanical  pump  for  vacuum  rough-down 
purpose  while  the  others  were  an  assortment  of  electrical, 
mechanical  and  water  feedthroughs.  The  water  feedthroughs 
were  mostly  made  in  the  laboratory,  whereas  the  electrical 
feedthroughs  were  generally  high  voltage,  hermetically 
sealed  and  of  a  commercial  kind.  The  test  chamber  was 
placed  directly  above  the  feedthrough  ring,  whose  inside 
diameter  was  equal  to  the  chamber  diameter,  so  as  to 
achieve  a  minimum  loss  in  pumping  conductance.  A  stainless 
steel  wire  mesh  was  placed  in  between  the  baffle  and  the 
feedthrough  ring  —  the  purpose  of  whi ch  was  to  trap  any 
dropped  objects.  Though  the  pumping  speed  then  was  reduced 
somewhat,  the  saving  in  system  downtime  well  compensated 
the  minor  drawback.  All  vacuum  surfaces  were  sealed  by 
single  Viton-A  "0"  rings  with  Apiezon-L  vacuum  grease 
except  the  12"  diameter,  aluminum  cover  plate  for  the 
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access  opening,  which  had  a  double  "0"  ring  seal,  since  it 
was  subjected  to  repeated  cycling.  Two  glass  view  ports 
were  provided  for  monitoring  the  test  chamber  conditions: 
one  (3”  diameter)  in  the  side  of  the  chamber  facing 
directly  at  the  target  surface  was  used  to  measure  the 
latter's  temperature  with  an  optical  pyrometer,  and  the 
other  one  (5"  diameter)  in  the  center  of  the  cover  plate 
of  the  access  opening  served  to  view  mechanical  alignment 
and  allow  source  tip  temperature  measurement. 

An  interior  view  of  the  test  chamber  with  the  bell 
jar  removed  is  shown  in  Fig.  8.  The  photograph  was  taken 
at  an  angle  from  the  top  to  show  all  components  and 
support  structures.  The  main  frame  consisted  of  four  1” 
diameter  stainless  steel  rods  that  were  mounted  vertically 
on  a  ring  fixed  to  the  inside  of  the  feedthrough  ring. 

Two  platforms  were  supported  by  main  frame  cross  bars 
that  also  provided  rigidity  to  the  structure.  The  upp_r 
platform  supported  the  mass  spectrometer  analyzer  tube 
and  its  magnet  as  well  as  the  drive  train  for  their  rota¬ 
tion  and  translation  motion.  The  drive  train  consisted  of 
gears  and  linkages  driven  by  Globe  AC  motors,  whose 
planetary  gear  boxes  were  completely  degreased  and  then 
lubricated  with  Silicone  D.C.  705  pump  oil.  The  rotation 
motion  of  the  mass  spectrometer  was  pivoted  about  the  center- 
line  axis  of  the  target,  thereby  allowing  measurements  of 
the  re-emitted  particle  flux  angular  distribution  from  the 
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targ ..  v.  .  ehe  principal  plane  defined  by  the  beam  axis 

and  caxget  normal.  The  position  of  the  source  housing, 

as  shown  in  Fig.  9  at  its  location  relative  to  the 

target  cetector  assembly,  nevertheless  limited  the  mass 

spectrometer  rotation  range.  For  example,  at  0^  =  45°, 

0°  ^  0D  ^  90°,  whereas  at  0^  =  35°,  10°  $  @D  £  90°. 

The  translation  notion  of  the  mass  spectrometer  allowed 

for  adjustment  of  the  target-to-detector  distance.  As 

shown  in  Section  II-C,  the  signal  level  is  inversely 

2 

jroportional  to  ®-TD-  A  shortaning  of  fcTD  leans  to  an 
increased  signal  level  and  thus  also  signal-to-noise  ratio. 
However,  the  price  paid  for  this  is  reduced  angular 
resolution. 

The  position  of  the  mass  spectrometer  was  monitored 
by  two  potentiometers,  turning  simultaneously  with  the 
rotation  and  translation  motors  and  forming  one  branch 
of  a  Wheatstone  bridge  completed  by  two  identical 
potentiometers  Iocs ted  external  to  the  vacuum  chamber. 

The  imbalance  of  the  Wheatstone  bridge  induced  by  the  pot 
rotation  when  the  mass  spectrometer  was  moved  was  manually 
nulled  by  the  external  pots,  thereby  establishing  the 
position  of  the  mass  spectrometer  in  accord  with  previous 
calibration. 

The  lower  platform  (Fig.  8)  supported  the  far  it 
assembly,  which  could  be  rotated  as  well  as  moved  vertically. 
The  up-and-down  motion  was  necessary  to  provide  the  target 
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with  one  (down)  position  which  did  not  interfere  with 
incident  beam  measurements.  The  beam  incident  angle 
was  varied  by  rotating  the  target  while  in  the  up  position, 
in  which  case  a  geared  potentiometer  was  engaged  by  a  pin 
in  the  target  support  shaft.  The  rc nation  angle  was  then 
measured  in  a  scheme  similar  to  that  for  the  mass  spectro¬ 
meter.  ei  could  be  varied  from  0°  to  90°.  The  motion  of 
the  tarc-et  was  actuated  manually  from  outside  of  the  vacuum 
system  by  means  of  a  Consolidated  Vacuum  Corporation  linear¬ 
rotary  feedthrough  (Model  SP-50)  located  at  one  of  the 
ports  of  the  feedthrough  ring.  A  6"  long  flexible  metal 
shaft  (Stow  Manufacturing  Co.  Model  13FC) ,  encased  in  a 
3/4"  diameter  tubular  housing,  served  as  the  coupling 
between  the  linear-rotary  feedthrough  and  the  target 
support  shaft. 

The  chopper  and  its  collimators  were  also  mounted  on 
the  lower  platform  between  the  target  and  the  mass 
spectrometer  detector.  Rotation  of  the  assembly  about 
the  target  centerline  was  provided  by  a  Globe  AC  motor. 
Alignment  of  tl.o  chopper  with  the  mass  spectrometer 
detector  was  accomp  ished  by  visually  matching  two  wire 
pointers  —  one  on  each  component  —  within  an  accuracy  of 
+  0.5°. 

A  schematic  diagram  of  the  test  chamber  interior 
arrangement  is  shown  in  Fig.  10. 

Fig.  11  shows  a  view  of  the  vacuum  system  from  the 
source  chamber  side.  The  source  chamber,  in  the  shape  of 
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a  tee,  was  evacuated  by  an  Edwards  High  Vacuum  Inc. 

300  liters/sec,  6"  vapor  booster  pump  (Model  9B3) 
using  Silicone  D.C.  704  oil  together  with  a  Freon  12 
cooled  copper  chevron  baffle.  One  side  of  the  tee 
served  as  an  access  port  whereas  the  other  side  was 
connected  to  an  aluminum  feedthrough  ring  similar  in 
design  and  function  to  that  of  the  test  chamber. 

Bolted  onto  the  feedthrough  **ing  was  a  rectangular 
stainless  teel  source  housing  that  protruded  into  the 
test  chamber.  The  source  support  assembly  could  be 
moved  vertically  as  well  as  in-and-cut  to  provide  both 
alignment  and  variation  of  source-to-collimator  distance. 

It  was  driven  manually  external  to  the  source  chamber  through 

rotary  feedchrougns  with  flexible  couplings  and  its 
position  was  calibrated  against  the  number  of  turns  of 
the  crank  handle.  The  source  collimator,  which  consisted 
of  a  0.003"  thick  tungsten  foil  with  a  3/32"  diameter  hole 
mounted  on  a  beryllium-copper  holder,  was  bolted  to  the 
center  of  the  end  plate  of  the  rectangular  housing,  which 
was  water-cooled  due  to  the  promixity  of  the  hot  source 
and  target.  The  source  chamber  was  completely  sealed  by 
Viton~A"0"  rings  from  the  test  chamber,  excepting  for  the 
collimator  hole.  A  3"  diameter  glass  view  port  was  pro¬ 
vided  in  the  center  of  the  cover  plate  for  the  access 
opening.  An  interior  view  of  the  sou:  :e  chamber  is  shown 
in  Fig,  12. 
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The  initial  pump-down  of  both  chambers  from  atmosphere 

was  done  by  a  NRC  100  CFM,  2-stage  rotary  gas  ballast  pump 

(Model  100-C) ,  which  also  served  as  the  backing  pump  for 

the  two  oil  diffusion  pumps.  The  exact  pumping  speeds  o_ 

the  diffusion  pumps  pertaining  to  the  present  experimental 

29 

setup  were  determined  by  Shih.  A  summary  of  the  results 
is  presented  in  Appendix  B.  Pressures  in  the  vacuum  system 
were  measured  by  two  General  Electric  hot  filamen  *  miniature 
ionization  gauges  (Model  22GT115)  with  Granville-Phillips 
gauge  controllers  (Model  236).  Both  gauges  were  mounted 
in  the  respective  feedthrough  rings  of  the  test  and  source 
chambers.  The  ion  gauge  calibration  for  different  gases 
is  given  in  Appendix  C. 

3.  Target  and  Its  Support 

The  target  used  in  this  experiment  was  fabricated 
from  ultra  pure  isotropic  graphite  sheets  (Grade  DFP-2) 
obtained  from  POCO  Graphite,  Inc.  The  graphite  sheets 
(0.0P""  thick,  4"  x  6")  were  shipped  individually  sealed 
in  air-tight  plastic  envelopes  to  avoid  cor  'lamination  and 
were  not  opened  until  ready  for  use.  Though  the  process 
of  purification  was  not  known,  the  total  impurity  content 
of  the  graphite  was  stated  to  be  less  than  5  ppm  with  only 
traces  of  Fe,  Si,  A1  and  Mg.  Except  for  degassing  under 
vacuum,  no  special  treatment  of  the  graphite  was  performed. 
Some  of  its  physical  properties  are  shown  in  Table  III-l. 
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The  target,  1  cm  x  2.6  cm,  was  clamped  in  the  holder 
as  shown  in  Fig.  13.  The  holder  consisted  of  two  water- 
cooled  copper  blocks,  each  with  a  screwed-on  flat  clamp, 
connected  by  an  aluminum  oxide  plate  which  also  served  to 
is*.  _te  the  two  blocks  electrically.  The  unit  was  mounted 
on  the  movable  shaft  with  another  Al203  plate  in  between  for 
electrical  insulation.  The  centerline  of  the  target 
coincided  with  the  axis  of  rotation  for  the  mass  spectro¬ 
meter  and  chopper.  A  spring  tension  was  applied  to  the 
target  to  prevent  it  from  warping  due  to  thermal  expansion; 
this  was  done  by  1'astening  the  upper  copper  block  to  the 
Al203  plate  through  a  spring  plate.  A  silver-plated 
beryllium  copper  sh  *■  was  connected  to  each  of  the  copper 
blocks  at  the  back  side  of  the  target.  The  sheets  served 
as  both  radiation  heat  shields  for  the  A12C>2  plate  and 
thermal  energy  reflectors. 
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The  target  was  heated  resistively  with  a  Hewlett- 

Packard  DC  power  supply  (Model  6269A) .  The  water  lines 

for  cooling  the  copper  blocks  (1/16"  diameter  flexible 

copper  tubing)  served  as  the  electrical  conductors  and 

they  were  covered  by  braided  fiber  glass  insulators.  The 

temperature  of  the  target  was  monitored  by  an  optical 

pyrometer  through  the  3"  diameter  view  port.  A  spectral 

emis-ivity  of  0.90  at  a  wavelength  of  0.65y  was  used  for 

the  graphite  temperature  calibration  though  the  actual 

emissivity  of  a  finely  machined  graphite  surface  varies 
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from  0.90  to  0.95  with  temperature  and  surface  finish. 

The  resulting  error  at  an  observed  temperature  of  1400°K 
was  about  10 °K  and  in  no  case  was  the  error  greater  than 
1%.  A  correction  for  sighting  through  a  clear  glass  window 
was  also  made  according  to  the  pyrometer  manufacturer's 
recommendation.  At  a  graphite  surface  temperature  of 
about  1600°K,  a  total  power  input  of  180  watts  from  the 
power  supply  was  required. 


4.  Chopper  and  Its  Collimator 


Instead  of  a  conventional  slotted  disc  mounted  on 


a  motor  shaft,  a  tuning  fork  chopper  was  used  for  modulating 


the  molecular  beam  and  re-emitted  particle  flux  from  the 


target.  There  are  several  advantages.  The  compact  size 
makes  it  relatively  easy  for  mounting  and  positioning. 


especially  in  the  present  case  of  rather  limited  space 
between  the  target  and  detector.  No  special  sealing  is 
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necessary  as  in  the  case  for  a  high  speed  motor,  since 
there  are  no  parts  that  need  lubrication.  The  driving 
voltage  of  the  tuning  fork  can  also  serve  as  a  reference 
point  for  the  signal  processing  scheme;  for  a  motor 
drive  chopper,  a  separate  setup  such  as  a  light-photocell 
combination  would  have  been  necessary.  The  tuning  fork 
chopper  has,  however,  the  drawback  of  a  fixed  frequency. 

It  cannot  be  readily  used  in  experiments  that  require  a 

2 

variable  chopper  frequency  as  in  the  case  of  Olander  et  al. 

The  tuning  fork  chopper  was  ordered  from  American  Time 
Products  Dept,  of  Bulova  Watch  Co.  The  unit  (Model  L40  HHB) 
with  frequency  of  400  Hz  was  supposed  to  operate  in  vacuum 
and  up  to  a  temperature  of  250 °C.  As  it  turned  out,  the 
unit  was  designed  only  for  a  bakeout,  not  operating, 
temperature  of  250°C.  Eventually,  an  almost  completely 
custom-made  unit  with  a  frequency  of  377  Hz  was  used.  A 
picture  of  the  chopper  is  shown  in  Fig.  14.  This  consisted 
of  a  fork  driven  by  two  electromagnets  located  at  about 
the  middle  of  each  prong.  Two  stainless  steel  blades, 
soldered  on  top  of  the  prongs  by  a  special  low  temperature- 
high  strength  silver-based  alloy  (Gemini  SA  alloy  from 
Nassau  Alloys  and  Research  Corp.),  served  as  the  chopper. 
They  were  placed  such  that  the  blades  just  overlapped  in 
the  stationary  mode,  giving  equal  open  and  shut  times 
(i.e.  ,  a  duty  cycle  of  0.5)  when  in  operation.  The  blades 
were  polished  on  the  side  facing  the  hot  target  to  minimize 
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heat  absorption  and  had  a  sandpapered  finish  on  the  back 
side  so  that  any  absorbed  heat  could  be  more  readily 
radiated  away  due  to  the  relatively  higher  emissivity. 

The  tuning  fork  was  bolted  to  an  L-shape  copper  bracket 
with  water-cooled  coils  silver-soldered  on  the  surface 
facing  the  target.  The  bracket#  excepting  for  a  small 
opening  for  the  passage  of  gas  molecules,  shielded  the 
entire  chopper  assembly  as  well  as  a  portion  of  the  mass 
spectrometer  ionizer.  A  tungsten  sh<“»t  of  0.012  cm 
thickness  with  a  0.508  cm  square  opening  was  mounted 
between  the  chopper  blades  and  copper  bracket  and  served 
as  the  chopper  collimator.  The  whole  assembly  sat  on  an 
aluminum  plate  that  rotated  about  the  target  centerline 
as  described  in  Section  III-A-2. 

The  electromagnets  were  driven  by  an  oscillator 
circuit  (ATP  Model  5A  driver)  and  the  exact  frequency  was 
determined  by  a  General  Radio  Company  Model  1191  digital 
counter.  The  chopper  amplitude,  viewed  with  a  GRC  Strobotac 
shining  through  the  6"  diameter  view  port,  was  adjusted  to 
cover  the  whole  collimator  opening.  The  shutter  function, 
expected  to  be  a  half  sine  wave,  was  measured  by  a  light- 
photocell  scheme.  Fig.  15  shows  the  measured  shutter 
function  together  with  the  driver  voltage  which  was  used 
as  a  reference  point  for  signal  processing.  It  can  be  seen 
that  the  chopper  lags  the  driver  voltage  by  approximately 
90°  (an  exact  phase  lag  of  90°  was  measured  by  a  lock-in 
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amplifier) .  This  phase  lag  was  not  unpredicted  duo  to 
the  inductive  nature  of  the  driving  force.  The  shutter 
function  was  confirmed  to  be  a  half  sine  wave  by  com¬ 
paring  the  magnitude  of  the  photocell  signal  as  displayed 
on  the  oscilloscope  with  that  of  the  output  of  the  same 
signal  as  processed  by  a  lock-in  amplifier  (which  measures 
the  root-mean-square  value  of  the  fundamental  component 
of  the  Fourier  expansion  of  an  input  signal) .  The  calcu¬ 
lated  value  based  on  the  maximum  amplitude  of  the  scope 
trace  agreed  with  the  lock-in  amplifier  output  to  within  2%. 

Tnough  the  chopper  assembly  was  completely  shielded 
from  the  hot  target  except  for  the  chopper  blades,  the  heat 
flux  still  had  an  effect  on  the  chopper  operation  neverthe¬ 
less.  The  chopper  frequency  and  amplitude  were  very  stable 
up  to  a  target  temperature  of  1700°K.  A  further  increase 
in  target  temperature  invariably  resulted  in  an  erratic 
behavior  of  the  chopper,  such  as  fluctuation  of  frequency 
and  drastic  reduction  in  amplitude.  Subsequent  reduction 
in  target  temperature  alleviated  this  problem.  No  definite 
explanation  could  '  e  found;  possibly  softening  of  the 
solder  led  to  this  effect  since  the  chopper  performance, 
especiaJly  the  amplitude,  depended  heavily  on  the  delicate 
operation  of  soldering  the  blades  onto  the  prongs.  Accord¬ 
ingly,  the  maximum  target  temperature  was  limited  to  1700 °K 
for  all  measurements  except  when  degassing  the  target. 


5.  Mass  Spectrometer 


The  mass  spectrometer  was  a  General  Electric  unit 
that  consisted  of  a  stainless  steel  analyzer  tube  (Model 
22PT120)  and  a  3  kilogauss  permanent  magnet  (Model  22PM103) 
for  ion  focussing  and  deflection.  The  analyzer  tube  con¬ 
sisted  of  an  exposed  Nier-type,  electron -bombardment  ion 
source;  a  90-degree  sector,  double-slit  magnetic  analyzer 
with  a  5  cm  radius  of  curvature;  and  a  10-stage  electron- 
multiplier  ion  detector.  Originally,  the  analyzer  tube 
was  designed  to  mount  directly  on  a  vacuum  system  through 
the  use  of  a  4  3/4"  O.D.  flange  at  the  ion  source.  In  the 
present  application  the  analyzer  was  modified  to  be  used 
inside  the  vacuum  system  and  supported  through  the  magnet. 

The  flange  was  removed  along  with  rerouting  of  some 
electrical  connections  to  permit  operation  as  a  "flow 
through"  ionizer.  The  original  electron  multiplier  was 
replaced  by  an  EMI  Gencom  Division  18-stage,  venetian- 
blind  type  multiplier  with  BeCu  dynode  surfaces  (Model  D1AT) 
for  higher  gain  and  stability.  Resistors  (glass-encased 
metal  film  type)  for  the  dynode  chain  were  all  spot-welded 
on  the  dynode  stages  and  decoupling  capacitors  were  installed 
between  the  last  three  stages  and  ground  to  avoid  the  need 
for  large  chain  currents  when  measuring  large  amplitude 
pulses.  The  current  output  end  was  redesigned  to  provide 
complete  shielding  of  the  multiplier  and  electrical  leads. 

The  stainless  steel  housing  of  the  multiplier  was  first 


-  -*■  -  .  -- - •  aa> i. 


ru. 


*i«.i  aa 


7frZ,~ZfF*mv** KTWSTr^-rxK  V*<^ 


95 

lined  by  a  layer  of  0.025  cm  thick  y-metal  to  minimise 
magnetic  pickup  and  then  followed  by  a  0.14  cm  thick 
copper  enclosure  with  two  BNC  connectors  at  the  end  for 
signet-  and  high  voltage  leads.  The  signal  lead  was  a 
special  coaxial  cable  with  two  braids  of  flat  strips  of 
silver-covered  copper  (Stripflex  cable  from  Times  Wire  and 
Cable  Company)  for  greater  shielding  efficiency,  and  it  was 
connected  to  the  input  of  the  amplifier/discriminator 
through  a  one-pin  BNC  feedthrough  to  provide  as  much 
isolation  as  possible  from  other  wirings.  An  overall  viev; 
of  the  mass  spectrometer  is  shown  in  Fig.  16. 

The  ion  source  consisted  of  an  ion  cage,  two  thoria 
coated  iridium  filaments  (one  being  a  spare) ,  draw-out  and 
deflection  electrodes  and  the  defining  slits  (entrance  and 
collimating) .  An  electron  current  was  drawn  from  the 
filament  to  the  ion  cage.  Inside  the  ion  cage  the  neutral 
atoms  and  molecules  diffusing  into  the  cage  were  ionized 
by  electron  bombardment.  The  ions  were  drawn  from  the 
ion  cage  by  a  constant  drawout  potential  and  accelerated 
through  ?  particular  ion  focussing  potential  (accelerating 
potential) .  The  beam  of  ions  thus  formed  was  focussed  by 
the  deflection  electrodes  and  passed  into  the  magnetic 
analyzer  through  the  defining  slits. 

The  ion  beam  leaving  the  ion  source  consisted  of 
particles  having  approximately  identical  energies  and  direc 
tion.  The  kinetic  energy  (K.E.)  of  the  ions  in  the  beam  is 
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defined  by: 


K.E. 


qv 


acc 


(3.1) 


When  a  charged  particle  moves  through  a  magnetic  field, 
it  is  subjected  to  a  force  acting  perpendicular  to  both 
the  field  and  the  direction  of  the  motion.  This  deflecting 
force  has  a  magnitude  given  by: 


F 


mag 


qvB 

c 


(3.2) 


and  equal  to  the  centrifugal  force  on  the  particle;  thus 


F 

mag 


(3.3; 


Combining  Eqs.  (3.1)  and  (3.3)  to  eliminate  v  yields: 


m  =  _1_  r2B2 
q  2c2  Vacc 


(3.4) 


Thus,  for  a  given  magnetic  field  strength,  the  mass  of 
those  particles  which  are  deflected  by  the  correct  amount 
so  as  to  strike  the  ion  collector  is  inversely  proportional 
to  the  ’on  accelerating  voltage.  The  mass  range  is  scanned, 
therefore,  by  varying  the  accelerating  potential. 

With  the  3  kilogauss  nagnet,  the  mass  spectrometer  was 
able  to  cover  a  mass  range  of  2-200  AMU  with  unit  resolution 
of  m/q  up  to  50. 
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The  electronic  controller  of  the  mass  spectrometer 
(GE  Model  22PC110)  consisted  of  a  chart  recorder  and 
three  separate  modules:  electron  emission  control, 
signal  readout  and  ion  sweep  control.  The  experiment 
demanded  the  viewing  of  a  particular  mass  peak  for  long 
periods  of  time,  and  the  circuit  drift  in  the  ion  sweep 
control  module  was  found  to  be  unacceptable  over  the 
necessary  time  intervals.  Subsequently,  only  the  electron 
emission  control  module  was  used  in  the  experiment.  This 
contained  the  emission  regulated  filament  power  supply  and 
the  controls  to  adjust  the  draw-out  and  two  deflectirm 
voltages,  which  formed  and  collimated  the  ion  beam.  The 
ion  sweep  control  module,  which  provided  the  variable  ion 
accelerating  voltage  for  electrostatically  scanning  the 
mass  range,  was  replaced  by  a  John  Fluke  Manufacturing  Co. 

DC  power  supply  (Model  407DR)  with  the  voltage  measured  by 
a  digital  voltm-  ter.  The  electron  multiplier  voltage 
generator  was  also  replaced  by  a  Hewlett-Packard  high 
voltage  DC  power  supply  (Model  6516A)  for  stability.  Al_ 
leads  were  connected  through  a  GE  8-wire,  high  voltage 
ceramic  feedthrough  (Model  22HN020).  A  picture  of  the 
instrument  panels  is  shown  in  Fig.  17. 

Under  electron  bombardment  in  the  ion  source,  some 
of  the  gas  molecules  are  dissociated  into  their  constituent 
atoms.  This  gives  rise  in  the  mass  spectrum  to  a  mass  peak 
for  each  fragment  of  the  parent  molecule.  These  peaks  along 


with  the  peak  resulting  from  the  parent  molecule  are  known  l 

as  the  "cracking  pattern"  for  that  species.  This  cracking  5 

j 

pattern  is  unique  for  each  gas  under  a  set  of  mass  spectro-  \ 

\ 

meter  operating  conditions  and  is,  in  effect,  used  to  j 

determine  the  relative  amounts  of  two  species ,  both  of 

which  make  a  major  contribution  to  the  spectrum  at  identical  ' 

masses.  An  overall  calibration  of  the  mass  spectrometer  for 
different  gas  species  is  given  in  Appendix  C. 

6 .  Amplifier/Discriminator 

The  output  current  from  the  mass  spectrometer  was 
fed  into  a  SSR  Instrument  Co.  amplifiar/discriminator 
(Model  1120) .  A  table  of  its  specifications  is  given  in 
Table  III-2. 


Table  III-2 


SSR  Model  1120  Amplif ier/Discrin. 

Input  Impedance 
Overall  Gain 

Discriminator  Sensitivity 
Rise  Time 

Pulse-Pair  Resolution 

Maximum  Input  Pulse  Amplitude 
for  Linear  Operation 

Minimum  Detectable  Pulse  Amplitude 

Input  RMS  Noise  Equiv.  (50  ohm) 

Pulse  Output,  Current  Source  Output 

Output  Pulse  Width 


nator  Specifications 

50  ohms 
2300  min. 

25  mV  to  500  mV,  internal 
screwdriver  adjustment 

6  asec 

10  nsec 

500  yV 
50  yv 

10  yv  (max) 

20  mA,  MECL  compatible 
10  nsec  between  10%  points 


Jr,  UtoriftA©  i 


The  amplifier  consisted  of  a  low  noise,  current-mode 
stage  having  a  50  ohm  input  impedance  followed  by  two 
current-mcde  stages  and  a  base-grounded  amplifier  having 
a  combined  amplification  of  350.  The  discriminator, 
preceded  by  a  discriminator  driver  having  an  amplification 
of  7,  was  an  internally  adjustable  tunnel  diode  discriminator 
having  an  adjustable  threshold  from  25  to  500  mV.  It  was 
followed  by  a  line  driver  to  permit  connection  over  long 
lengths  of  coaxial  cable.  The  A/D  output  could  be  used 
directly  to  drive  high  speed  counting  equipment  requ.ring 
1  volt  or  less  signal  amplitude. 

The  A/D  was  contained  within  a  heavy  walled  aluminum 
housing  with  multiple  internal  shields.  The  input  connector 
was  a  panel  BNC  jack,  and  all  output  and  power  connections 
were  made  through  an  Amphenol  connector  incorporating  both 
coaxial  connecting  pins  for  signa.  =s  and  normal  pins  for 
power.  The  unit  should  be  mounted  in  as  close  proximity  to 
the  electron  multiplier  anode  connection  as  is  feasible. 

This  was  done  by  strapping  the  A/D  immediately  adjacent  to 
the  signal  output  feedthrough  and  using  a  6”  double-shielded 
interconnecting  cable.  Careful  filtering  and  decoupling  of 
the  electron  multiplier  high  voltage  power  supply  is  mandatory 
to  avoid  coupling  spurious  pulses,  ripple  and  noise  into  riie 
electron  multiplier  anode  circuit  from  the  supply.  This 
was  accomplished  by  by-passing  the  last  three  dy nodes  to 
ground  by  high  frequency,  metal-encased  capacitors  of  0.01  uF. 
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The  power  requir  1  for  the  A/D  was  supplied  by  the  * 

3 

Digital  Synchronous  Computer.  1 

7 •  Digital  Synchronous  Computer  i 

The  sign?.!  pulses  were  processed  by  a  SSR  Digital  \ 

Synchronous  Ccv.puter  (Model  1110) .  It  is  a  highly  versatile 
digital  courting  and  computing  instrument  that  ccr.sists  of: 

A.  Two  High  Speed  8-digit  scalers. 

B.  A  crystal  controlled  Master  Clock  for  accurate 

timing.  ■ 

C.  A..  Arithmetic  Processor  which  coi  verts  the 
counter  outputs,  their  sum  and  their  difference 
into  "floating  point'1  form  (3  significant  digits 
and  exponent) . 

D.  A  Preset  Counter  which  preselects  chopper  cycles, 
time,  or  input  events,  depending  on  the  mode  of 
operation. 

E.  A  Preset  Counter  for  chopper  sampling  time. 

F.  Front  panel  controls  and  display. 

Its  specifications  are  given  in  Table  III-3. 


Table  III-3 

SSR  Model  1110  Digital  Synchronous  Computer  Specifications 


Pulse  Pair  Resolution  12  nsec 

Maximum  Count  Rate  85  MHz 

Crystal  Clock  Accuracy  .01%  (long  term) 

1  part  in  10 8  (short  term) 

O 

Gating  Symmetry  (Chopper  Mode)  30  nsec  or  1  part  in  10  , 

whichever  is  greater 

Dead time  between  Successive 

Main  Cycles  67-85  usee 

Capacity  of  "A”  and  "B"  Counters  10®  counts  BCD 


67-85  usee 


Digital  Output 


Preset  Scale  Factor 


Input  Signal  Pulses  Negative  1  V  min  (5  V  max) 

Duration  (FWHM)  10  nsec  (min) 

Chopper  Reference  Signal  ±  50  mV  min  (5  V  max) 

(or  Sync  Trigger)  Adjustable.  Inputs  may  be 

operated  in  parallel  or 
individually. 

Digital  Output  Four  16-bit  full  parallel 

words . 

Preset  Scale  Factor  Three  significant  figures 

plus  exponent  for  time  (usee) , 
counts,  or  chopper  cycles. 
Range:  001  to  999  x  10^. 

The  DSC  was  operated  in  the  "chopper"  mode  as  described 
in  Section  II-B.  An  external  "chopper  trigger"  was  used  to 
create  two  counting  periods  per  chopper  cycle.  During  each 
counting  period,  counts  were  stored  either  in  counter  "A"  or 
counter  "B".  Usually  when  the  chopper  was  open,  counts 
were  stored  in  counter  "A"  (Signal  +  Background)  and  when 
it  was  closed,  counts  were  stored  in  counter  "B"  (Background) . 

The  total  counting  time  was  determined  by  selecting  the 

9 

number  of  chopper  cycles  (up  to  10  )  on  the  front  panel  preset 
switch  whereas  the  counting  duration  for  each  chopper  cycle 
was  set  on  the  rear  panel  preset  switch  labeled  "Chopper 
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Sampling  Time"  (1  Msec  to  1  sec).  Upon  finishing  a  counting 
cycle,  the  results  (A,  B,  A+B,  A-B)  were  selectively  dis¬ 
played  in  three  significant  figures,  plus  exponent,  plus  sign. 
An  automatic  mode  was  available,  in  which  case  each  time  a 
counting  cycle  ended  and  the  display  was  updated,  a  new 
counting  cycle  would  begin. 

Since  the  driving  voltage  of  the  tuning  fork  chopper 
was  used  to  trigger  the  DSC  sampling  gates,  a  phase  adjusting 
scheme  was  necessary  to  allow  for  the  difference  in  phase  of 
the  chopper  opening  and  the  driver  voltage,  and  particle 
transit  time.  Originally,  a  delay  circuit  was  supposed  to 
be  incorporated  in  the  DSC  and  could  be  activated  by  the 
addition  of  a  potentiometer  at  appropriate  points  in  the 
circuit.  It  was  found,  however,  that  such  an  approach 
brought  about  an  instability  in  the  circuit  and  caused  an 
erratic  behavior  of  the  unit.  This  phenomenon  was  confirmed 
by  SSR  designers ,  who  conceded  that  the  phase  delay  mode  of 
the  unit  was  not  yet  fully  debugged.  Subsequently,  a 
Rutherford  Electronics  Co.  pulse  generator  (Model  B7B)  was 
used  to  adjust  the  opening  times  of  the  DSC  counter  gates. 

A  timing  diagram  for  the  phase  adjusting  scheme  is  shown 
in  Fig.  18.  The  driving  voltage  of  the  tuning  fork  was  used 
to  trigger  the  pulse  generator,  which  in  turn  produced  a 
rectangular  pulse  of  adjustible  amplitude,  duration  and  delay 
The  pulse  was  adjusted  to  equal  on-off  time  and  fed  into  the 
DSC  as  trigger  with  the  delay  set  to  maximize  signal  output. 
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The  exact  delay  time  was  monitored  continuously  by  a 
General  Radio  interval  counter.  A  schematic  diagram  of 
the  signal  processing  procedure  is  shown  previously  in 
Fig.  4. 

8.  Beam  Source 

A  cross-section  diagram  of  the  beam  source  is 

shown  in  Fig.  19.  The  source  plenum  chamber  consisted 

of  a  closed  end  1/4"  O.D.  thorium  oxide  (thoria)  tube  of 

1/32"  wall  thickness  and  2  1/2"  in  length.  A  0.010" 

diameter  hole  in  the  closed  end,  drilled  by  means  of  an 

ultrasonic  wire  drilling  technique  in  t;  e  M.I.T.  Ceramics 

Laboratory,  served  as  the  source  orifice.  The  front  1  1/4" 

portion  of  the  tube  was  heated  by  a  helically  wound  tungsten 

ribbon  filament  (0.010"  x  0.025"}  of  approximately  12  coils. 

The  whole  heating  section  was  protected  by  a  spiral  radiation 

shield  for  its  compactness,  effectiveness  and  convenience  of 
89 

construction.  The  shield  was  made  by  rolling  a  sheet  of 
1  1/2"  x  15"  x  0.005"  thick  tungsten  into  a  spiral  of 
approximately  six  layers  each  separated  by  two  crinkled 
tungsten  strips.  The  shield,  about  3/4"  in  diameter,  was 
supported  by  two  alumina  washers  fitted  over  the  thoria  tube. 
The  thoria  tube  was  adapted  to  a  water-cooled  BeCu  block 
with  two  Swagelok  fittings  and  a  pressure  tap.  The  metal- 
to-thoria  seal  in  the  source  fitting  was  achieved  by  two 
gold  "0"  rings  with  an  indium  "0"  ring  sandwiched  in  between. 
This  solution  was  arrived  at  after  repeated  experimentation. 
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At  high  source  operating  temperatures,  the  indium  "0"  ring 
had  actually  melted  but  was  retained  in  the  two  gold  "O" 
rings  gap  by  surface  tension,  thereby  maintaining  an  effec¬ 
tive  seal.  The  gas  inlet  line  was  connected  to  the  other 
Swagelok  fitting  and  the  source  holder  was  bolted  onto  an 
all  stainless  steel  source  support  carriage  assembly  that 
provided  fore-aft  and  up-down  motion  with  respect  to  the 
source  collimator  as  described  in  Section  III-A-2.  Two 
views  of  the  source  as  mounted  on  the  carriage  assembly  are 
shown  in  Fig.  20. 

The  tungsten  ribbon  filament  was  heated  by  a  homemade 
AC  power  supply  capable  of  putting  out  5000  watts.  All 
electric  power  connections  in  the  vacuum  chamber  were  1/4" 
diameter  solid  or  stranded  copper  leads  for  good  heat 
conduction,  since  there  were  no  provisions  for  water  cooling. 


9.  Gas  Inlet  System 


A  schematic  diagram  of  the  source  gas  inlet  system 


is  shown  in  Fig.  21.  Ultra  pure  grade  molecular  oxygen  was 


supplied  by  AIRCO  that  needed  no  further  purification  for 
use  in  this  experiment.  A  lot  analysis  of  the  gas  impurities 
is  shown  in  Table  III-4. 
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Table  III-4 

Impurity  Contents  of  0^, 

Impurity  PPM 

Kr  10.5 

N2  3.8 

Xe  0.276 

Hydrocarbons  <  0.5 

H20  0.8 

Gas  from  the  cylinder  first  passed  through  a  cold 
trap  that  was  employed  to  condense  water  vapor  and  other 
condensibles  from  regular  grade  oxygen  that  was  used 
earlier.  The  gas  flow  rate  was  then  controlled  by  a 
Granviile-Phillips  "variable  leak"  diaphragm  valve 
(series  203)  with  conductance  continuously  variable  from 
100  cc/sec  to  less  than  10  cc/sec.  By  setting  the  up¬ 
stream  supply  pressure  at  10  psig,  the  valve  was  found  to 

be  quite  satisfactory  in  maintaining  a  typical  flow  rate  of 
-4 

2  x  10  cc/sec.  Any  variation  of  the  flow  rate,  however, 
required  an  adjustment  time  of  up  to  half  an  hour.  The 
source  tube  pressure  was  measured,  at  the  point  of  the 
pressure  tap  in  the  source  holder  block,  by  a  Wallace- 
Tiernan  6"  dial  absolute  pressure  gauge  (0.1  -  20  torr) , 
which  was  calibrated  against  a  trapped  McLeod  gauge.  A 
calculation  showed  that  less  than  a  1%  correction  was 
necessary  to  allow  for  the  pressure  drop  from  the  pressure 
tap  point  to  the  source  tip.  The  gas  supply  line  was 
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evacuated  by  a  Welch  mechanical  pump  with  a  cold  trap  to 
prevent  backs treaming  of  pump  oil.  All  gas  lines  were 
soft  copper  tubings,  which  had  been  thoroughly  cleaned  by 
trichloroethylene  and  acetone.  The  valves  were  Nupro 
vacuum  bellows  valves. 

10.  Optical  Pyrometer 

The  temperatures  of  the  source  and  target  were 
monitored  by  an  optical  pyrometer  made  by  the  Pyrometer 
Instrument  Co. ,  Inc.  The  "Pyro  Micro-Optical  Pyrometer" 
was  a  telescopic,  disappearing  filament  type  with  vernier 
horizontal  and  vertical  adjusting  gears.  It  was  capable 
of  measuring  temperatures  of  targets  smaller  than  0.001" 
in  diameter  and  had  a  direct  reading  galvanometer  with 
overlapping  black -body  scales  from  700  -  3200 °C.  The  gal¬ 
vanometer  was  compensated  against  variations  in  ambient 
temperatures  and  had  a  suppressed  zero  in  order  to  permit 
maximum  use  of  the  full  length  of  the  scale.  Temperature 
measurements  were  reproducible  within  two  or  three  degrees 
Centigrade. 

An  optical  pyrometer,  however,  determines  only  the 

brightness  temperature  of  an  incandescent  body  at  the  mean 

effective  wavelength  between  the  brightness  temperature 

and  the  color  temperature  of  the  body.  If  the  spectral 

emissivity  of  the  body  is  known,  the  temperature  of  the 

body  may  be  determined  from  the  following  equation  using 

90 

Wien’r  radiation  approximation  : 
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T 


An  ex  (T) 
e 


(3.5) 


where  X is  the  mean  effective  wavelength 

(T)  is  the  spectral  emissivity  at  Xe 

C0  is  the  second  radiation  constant,  1.432  cm-°K 
an^  Tg  is  the  brightness  temperature  of  the  body  at  A^ 


A  tabxe  listing  the  true  temperatures  corresponding  to 
various  emissivities  and  brightness  terperatures  for 
Ae  =•  0.65m  was  provided  by  the  pyrometer  manufacturer. 

Beam  Properties  and  Source  Performance 

The  oxygen  beam  was  formed  by  collimating  the  effusive 
flux  from  the  circular  source  orifice.  To  ensure  that 
Knudsen  flow  existed,  both  the  source  pressure  and  tempera 
ture  were  maintained  at  levels  such  that  the  source  gas 
mean-free-path  was  greater  than  or  equal  to  the  critical 
orifice  dimension  —  in  this  case  the  wall  thickness,  L, 
since  it  was  greater  than  the  orifice  diameter,  D.  Namely 


For  0,  the  criterion  for  Knudsen  flow  (Eq.  (3.6))  is  less 

stringent/  since  6  <  6  .  Accordingly,  the  following 

o  o2 

discussion  is  restricted  to  molecular  oxygen  only.  Through¬ 
out  the  experiment,-  Eq.  (3.7)  was  nearly  always  adhered  to 
except  in  some  cases  when  the  source  was  operating  at  room 
temperature.  In  order  to  achieve  X  /L  ^  1  at  room  tempera- 
ture,  ps  ^  0.10  torr.  The  drawback  for  this  low  source 
pressure,  of  course,  is  a  low  centerline  intensity,  which 
means  a  low  signal  level.  Since  the  calculation  of  the 
absolute  reaction  probability  (as  given  in  Eq.  (2.52)) 
depends  partly  on  the  assumption  of  a  Knudsen  source  to 
interpret  the  measured  incident  beam  signal,  the  only  data 
that  were  used  to  calculate  such  quantities  were  those  that 
satisfied  Eq.  (3.7).  However,  for  measurements  that  did 
not  require  a  precise  knowledge  of  the  incident  beam  (such 
as  the  relative  angular  flux  distribution  from  the  target) , 
the  source  was  operated  at  higher  pressures  to  save, 
primarily,  signal  integration  time.  Such  an  approach  was 
justified  by  comparing  the  angular  distributions,  measured 
with  the  source  operated  at  Kn  =  1  and  at  Kn  as  low  as  1/4; 
no  difference  of  the  angular  distributions  was  observed. 

The  lower  limit  of  Kn  =  1/1  was  chosen  to  avoid  both  an 
excessive  increase  in  the  background  0 2  pressure  and  a 
radica^  departure  from  the  ideal  channel  source.  The 


departure  from  an  ideal  source  at  low  Kn  is  illustrated  in 
Fig.  22  that  shows  the  detected  incident  beam  centerline 
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signal  versus  source  pressure  for  a  room  temperature  source. 
The  straight  line  is  the  theoretical  centerline  intensity 
for  Knudsen  flow,  which  is  directly  proportional  to  the 
source  pressure  as  can  be  seen  in  Eq.  (2.4).  As  expected, 
the  measured  signal  at  higher  source  pressure  fell  short  of 
the  theoretical  value  due  to  the  departure  from  Knudsen 
flow.  This  method  of  increasing  centerline  intensity, 
though  resulting  in  a  small  loss  in  beam  directivity  or 
"peaking",  was  frequently  used  in  practical  application  of 
a  channel  source  to  fully  utilize  the  source  chamber  pumping 
capacity.  Theory  has  been  developed  to  predict  the 
variation  of  beam  intensity  with  source  pressure. 

Interested  readers  are  referred  to,  for  examples.  Refs.  91 
and  92. 

The  direct  experimental  verification  of  the  centerline 
intensity  of  a  channel  source  requires  a  difficult  absolute 
calibration.  Instead,  an  indirect  method  of  checking  the 
total  leak  rate  and  determining  the  angular  flux  pattern 
was  used  to  evaluate  the  source  performance. 

The  total  flux  out  of  an  effusive  source  is 
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if  the  Clausing  factor,  K,  remains  unchanged  as  will  be 
the  case  when  the  Knudsen  flow  criterion  (i.e.,  Kn  >  1) 
is  satisfied  at  all  times. 

The  relationship  in  (3.9)  was  tested  by  providing 

a  fixed  leak  through  the  Granville-Phillips  “variable 

leak"  valve  l i.e.  ,  ensuring  a  constant  Fg)  with  the  source 

at  initial  conditions  of  p  =0.1  torr  and  T  =  300°K 

s  s 

(which  satisfied  the  Knudsen  flow  criterion) .  The  source 
temperature  was  then  raised  and  the  corresponding  source 
pressure  was  recorded,  with  ample  waiting  time  in  between 
change  of  source  temperatures  to  ensure  equilibrium  condi¬ 
tions.  The  results  are  shown  ir.  Fig.  23  and  they  clearly 
verified  the  relationship  in  (3.9).  The  implication  of 
such  a  result  was  that  the  source  functioned  as  an  effusive 
source  within  the  pressure  and  temperature  range  tested 
(which,  incidently,  satisfied  Kn  >  X  at  all  times).  It 
should  be  pointed  out  that  there  was  practically  no 
dissociation  of  the  oxygen  under  these  test  conditions. 
Otherwise,  the  result  could  be  interpreted  differently, 
as  the  presence  of  0  could  lead  to  a  different  leak  rate. 

One  could  further  check  the  Clausing  factor  by 
measuring  the  total  leak  rate  out  of  the  source.  This 
could  be  done  quite  easily  by  feeding  the  source  from  a 
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constant  pressure  reservoir  and  measuring  the  volumetric 
change  of  the  reservoir.  However,  due  to  the  extremely 
low  flow  rate  of  the  source  and  the  requirement  of  an 
upstream  pressure  of  more  than  one  atmosphere  for  satis¬ 
factory  operation  of  the  GP  "variable  leak"  valve,  a 

reasonable  volumetric  change  of  the  reservoir  (e.g. ,  a 
3 

change  of  10  cm  )  would  require  a  flow  time  of  more  than 
20  hours.  Since  an  accurate  knowledge  of  the  total  flux 
rate  from  the  source  was  not  important  as  compared  to  the 
assurance  of  an  effusive  source  in  this  experiment,  this 
measurement  was  by-passed. 

The  source  flux  angular  distribution  at  different 
source  temperature  is  shown  in  Fig.  24.  No  absolute 
measurement  of  the  intensity  is  necessary  in  this  case. 

The  angular  distribution  was  measured  by  rotating  the 
mass  spectrometer  with  respect  to  the  target  axis.  The 
geometry  i?  also  shown  in  Fig.  24.  The  angle  of  rotation 
of  the  detector,  4>D ,  and  the  angle  measured  from  the  source 
tip,  4>s,  are  related  by 


sin  4>d 

tan  =  “os  <}>D  +  *ST/*TD 


(3.10) 


The  detector  width  was  set  by  the  amplitude  of  the  chopper. 

At  <p  -  0,  the  detector  width  corresponded  to  a  A<t>  of  0.6° 
s  s 

Due  to  the  geometry  of  the  detector  with  respect  to  the 
source,  <f>g  was  limited  to  cover  the  range  from  -3  1/2°  to  7 
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Though  a  far  from  complete  pattern  of  the  source  emission 
v»as  obtained,  the  4>  range  measured  was  sufficient  to 
evaluate  the  uniformity  of  the  incident  flux  on  the  target, 
since  the  target  subtended  only  a  source  angle  of  less 
than  6°.  The  angular  flux  distribution  for  an  ideal  cosine 
source  is  plotted  in  Fig.  24  for  comparison.  As  can  be 
seen,  the  angular  flux  distribution  of  the  source  at  the 
temperature  range  tested  (300°  to  2200°K)  deviated  little 
from  the  ideal  Knudsen  source;  it  was  more  lobular,  if  any, 
as  would  be  expected  for  a  channel  source.  In  taking  the 
data,  the  O-  mass  peak  (^  =  32)  was  monitored  for  T  ^  1800 °K 
whereas  at  Tg  =  22G0°K  the  0  mass  peak  (~  =  16)  was 
monitored.  Though  it  cannot  be  said  that  the  signal  level 
at  j!p  -  16  represents  solely  the  atomic  oxygen  flux  in  the 
beam  (as  will  be  discussed  later) ,  the  measurement  nonethe¬ 
less  is  an  indication  of  the  0  level.  In  any  case,  for 
the  small  source  angle  subtended  by  the  target,  the  non¬ 
uniformity  of  the  incident  flux  on  the  target  can  be  quite 
safely  neglected. 

It  is  quite  difficult  to  separate  a  recorded  mass  peak 

signal  level  in  a  mass  spectrograph  into  fractions  that 

are  contributed  by  different  mass  species.  For  example,  at 

m/q  =  16,  even  with  the  presence  of  only  0  and  Oj,  the 

signal  level  represents  o  nstituents  of  0  and  0  ,  the 

2 

latter  of  which  arises  from  both  O  and  the  dissociative 
ionization  of  Oj,  namely 
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0  +  e  •*■  0  +  2e 


02  +  e  -»■  0  +  0  +  2e 


(3.11) 


Additional  contribution  to  the  signal  level  is  possible 
with  the  presence  of  other  species  such  as  water  vapor. 

To  avoid  such  complication  in  interpreting  the  mass 
spectrograph,  the  ratio  of  0-to-C>2  signal  level  (more 
precisely  the  ratio  of  m/q  =  16  and  m/q  =  32  mass  peak 
signal  levels)  was  used  to  determine  the  degree  of  oxygen 
dissociation  in  the  beam.  In  this  scheme,  vhe  0/0 2  ratio 
was  first  obtained  for  a  pure  02  beam,  which  established 
the  "cracking  pattern"  of  02  under  the  operating  conditions 
of  the  present  mass  spectrometer.  The  0/02  ratio  was  then 
measured  for  a  dissociated  oxygen  beam.  Comparisons  were 
made  between  the  two  0/0 2  ratios  to  infer  the  degree  of 
dissociation  of  the  dissociated  oxygen  beam. 

The  0/02  ratio  for  pure  02  was  found  to  be  equal  to 

0.087.  This  value  was  constant  (within  3%)  for  both 

"stagnant"  02  as  occurred  in  02  background  gas  and  "directed" 

02  as  in  the  02  beam  flux  for  source  temperatures  from  300° 

to  1800°K.  This  result  was  contrary  to  that  obtained  by 
76 

Lam  ,  who  found  that  monoenergetic  molecular  oxygen  beams 
generally  gave  a  smaller  0/C>2  ratio  than  that  of  a  stagnant 
system  of  the  same  composition.  That  discrepancy  was 
attributed  to  the  present  case  of  a  relatively  large  mass 
spectrometer  ion  draw-out  energy  (approximately  200  eV) 
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compared  to  the  beam  energy  (maximum  0.4  eV)  and  the 
Maxwellian  velocity  distribution  in  the  beam.  After 
ionization  and  being  subjected  to  the  draw-out  potential, 
a  particle  has  practically  no  memory  of  its  previous  energy 
and  direction,  making  them  irrelevant  to  the  detected  0/C>2 
ratio. 

In  principle,  one  could  calculate  the  atomic  oxygen 
level  in  the  beam  once  the  cracking  pattern  of  02  was 
established  and  provided  that  such  parameters  as  ionization 
cross-section,  mass  spectrometer  ion  draw-out  efficiency 
and  electron  multiplier  dynode  ion-to-electron  conversion 
efficiency  were  accurately  known.  However,  such  parameters 
could  only  be  obtained  through  calibration,  which  in  the 
present  case  was  not  possible.  Instead,  a  calculation  of 
the  0/0 2  ratio  was  made  for  a  dissociated  beam  at  a 
particular  source  pressure  and  temperature  assuming  a 
mass  spectrometer  sensitivity  for  the  atomic  oxygen,  as 
defined  in  Eq.  (2.23).  This  assumed  mass  spectrometer 
sensitivity  for  0,  a r  educated  guess  to  begin  with,  was 
used  to  calculate  subsequent  0/02  ratios  for  different 
beam  compositions  as  generated  by  different  source  condi¬ 
tions.  Comparisons  were  made  between  the  "calculated"  and 
measured  0/02  ratios  while  the  mass  spectrometer  0 
sensitivity  was  continuously  adjusted  to  provide  a  most 
consistent  fit.  Such  a  table  showing  the  calculated  and 
measured  0/02  ratios  that  covers  the  entire  range  of  source 
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conditions  that  was  used  to  generate  atomic  oxygen  is 
given  in  Table  III-5. 


Table  III-5 

Dissociation  Level  as  Shown  by  0/0  ^  Ratio 


(torr) 

TS(°K) 

F/F 
°  °2 

°I°2- 

Calculated 

Ratio 

Measured 

0.10 

2220 

0.49 

0.1311 

0.128 

0.20 

2205 

0.31 

0.1139 

0.112 

0.24 

2200 

0.27 

0.1091 

0.111 

0.26 

2110 

0.14 

C . 1031 

0.100 

*  * 

For  the  calculated  0/0  _  ratio  in  the  table,  an  S  /S  ratio 
of  1/8  was  used.  This  value  appears  quite  reasonable  in 
view  of  the  reduced  ionization  cross-section  of  0  (approxi¬ 
mately  equal  to  1/4  of  the  ionization  cross-section  of  02;  , 
and  the  fact  that  the  mass  spectrometer  output  was  maximized 
on  the  02  rather  than  the  0  peak.  In  measuring  the  0/02 
ratio,  the  mass  spectrometer  was  always  "tuned"  to  the  02 
peak.  This  provided  a  consistent  basis  for  the  measurement, 
since  it  was  possible  to  arrive  at  a  different  0/02  ratio 
for  the  same  beam  composition  if  the  mass  spectrometer  was 
"tuned"  to  the  0  peak  instead.  The  consistency  of  the 
measured  and  calculated  0/02  ratios  over  the  entire  source 
operating  range,  as  shown  in  Table  III-5,  clearly  indicates 
the  satisfactory  operation  of  the  source  in  regard  to  the 
equilibrium  dissociation  level  of  oxygen. 


The  phase  lag  of  the  incident  beam  due  to  the  particle 

transit  time  was  measured  for  different  source  t.  _  jratures 

and  is  shown  in  Fig.  25.  The  results  agreed  quite  well 

83 

with  the  theoretical  predictions  of  Harrison  et  al. 

However,  due  to  the  circuit  drift  of  the  pulse  generator 
and  the  rather  insensitive  nature  oh  the  phase  measurement, 
the  measured  phase  shift  was  only  accurate  to  within  +  3". 
Accordingly,  it  was  not  possible  to  separate  out  the  0  and 
C>2  cracking  signals  at  m/q  =  16  for  a  dissociated  beam  by 
a  phase  lag  discrimination  technique,  since  at  a  source 
temperature  of  2200 °K  a  phase  lag  difference  of  only  4° 
existed  for  monoenercetic  O  and  02  beams.  With  the 
amplitude  dispersion  effect  for  a  Maxwellian  beam,  it  was 
impossible  to  distinguish  0  and  C>2  by  the  phase  lag  method. 

Up  to  a  source  temperature  of  1800°K,  the  source 
functioned  as  expected  and  the  only  necessary  precaution 
was  an  allowance  for  ample  time  for  the  source  to  reach 
thermal  uniformity  during  the  heat-up  and  cool-down 
processes.  A  four-hour  heat-up  time  was  required  to  avoid 
cracking  of  the  thoria  tube  due  to  thermal  shock,  while 
a  two-hour  period  was  sufficient  for  the  cool-down  process. 
When  the  source  was  operated  in  the  oxygen  dissociative 
regime  (i.e.,  source  temperature  greater  than  2100°K) , 
rapid  oxidation  of  the  tungsten  filament  by  the  background 
oxygen  in  the  source  chamber  invariably  took  place.  At 
a  source  temperature  of  2200°K  and  source  pressure  of  0.20 
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torr,  the  filament  burnt  out  after  approximately  10  hours 
of  continuous  operation.  Efforts  spent  to  prolong  the 
filament  life  by  using  heftier  tungsten  filament  and 
shielding  it  from  direct  impingement  of  the  reflected 
oxygen  from  the  source  collimator  were  of  little  avail; 
a  lengthening  of  only  a  few  hours  of  the  filament  life 
was  achieved.  This  problem  of  limited  filament  life  at 
£  2100°K  was  finally  remedied  by  positioning  the  source 
close  to  the  source  collimator  so  that  virtually  all  oxygen 
particles  effused  into  the  test  chamber.  This  remedy 
would  ordinarily  be  impossible,  since  the  increase  in 
noise  level  in  the  test  chamber  due  to  the  background 
oxygen  would  most  likely  be  intolerable.  Fortunately, 
due  to  a  combination  of  extremely  low  pressure  of  the 
source  (0.1  ^  0.25  torr)  and  a  huge  pumping  speed  in  the 
test  chamber  (10,000  liters/sec),  only  a  longer  integration 
time  and  a  slight  correction  had  to  be  made,  as  discussed 
previously  in  Section  II-C.  The  source  has  been  operated 
under  this  condition  for  more  them  100  hours  and  exhibited 
neither  a  visible  change  nor  change  of  electrical  resistance 
of  the  tungsten  filament.  The  minor  drawbacks  of  longer 
integration  time  and  a  correction  in  data  reduction  well 
compensated  the  formidable  task  of  replacing  the  filament, 
which  was  both  uneconomical  and  time  consuming. 
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C.  Experimental  Program  and  Procedure 
1.  Experimental  Program 

From  the  initiation  of  the  present  study,  it  had 
been  realized  that  a  major  effort  was  to  be  spent  on  the 
experimental  hardware  —  in  particular,  the  high  temperature 
source  and  the  pulse  counting  equipment.  Before  the 
present  vacuum  system  was  available,  source  tests  had  been 
conducted  in  a  small  vacuum  system  in  regard  to  the  maximum 
operating  temperature,  temperature  gradient  between  the 
heating  element  and  thoria,  filament  life,  metal-to-thoria 
seal,  and  so  on.  The  small  vacuum  system  was  quite  useful 
in  this  respect  —  in  particular,  due  to  its  short  vacuum 
cycling  time,  many  ideas  could  be  tried  out  rather  readily. 
Though  an  understanding  of  the*  source  performance  was  far 
from  complete,  the  results  obtained  in  that  period  proved 
to  be  quite  valuable  in  arriving  at  the  final  design  that 
was  used  in  the  experiment. 

At  the  time  when  the  vacuum  system  was  ready  for  use, 
the  most  promising  source  configuration  consisted  of  two 
concentric  thoria  tubings  forming  the  annular  Knudsen 
cavity  with  a  tungsten  rod  inside  the  smaller  tubing 
serving  as  the  heating  element.  This  source  was  installed 
in  the  source  chamber  and  much  of  the  room  temperature  beam 
data  was  obtained  with  that  arrangement.  At  high  source 
temperatures,  cracking  of  the  tubing  occurred  as  a  .asult 
of  differential  thermal  expansions.  Eventually,  the  source 
described  in  Section  III-A  was  developed. 
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In  view  of  the  low  beam  intensity,  transient  studies, 

29 

as  carried  out  by  Shih,  were  not  attempted.  Instead, 
the  emphasis  was  placed  on  experimental  accuracy  and  the 
dependence  of  the  reaction  on  relevant  parameters.  In 
order  to  separate  out  the  individual  effects  of  0  and  O 2 
on  the  oxidation  rate  when  using  a  mixed  beam,  the  oxidation 
rate  of  02  had  to  be  accurately  determined.  In  fact,  this 
was  the  first  objective  when  the  pulse  counting  system  was 
installed.  Comparison  of  room  temperature  C>2  results  with 
those  of  Shih  was  used  to  evaluate  the  detection  system. 

Hot  02  beams  were  used  after  the  detection  system  proved 
to  be  satisfactory.  The  dissociated  beam  measurements  were 
carried  out  during  the  latter  part  of  the  experiment. 

Although  "fast"  time-response  surface  phenomena  could  noc 
be  detected,  a  time-dependent  nature  of  the  oxidation  rate 
for  both  0  and  C>2  was  observed. 

While  it  was  planned  earlier  to  also  study  the 
oxidation  rate  of  pyrolytic  graphite  (made  by  depositing 
benzene  vapor  on  a  tungsten  or  isotropic  graphite  substrate) , 
an  expected  drop  in  signal  level  (which  meant  a  higher 
sensitivity  requirement  than  the  present  detection  system) 
precluded  this  study. 

2.  Mechanical  Alignment 

The  mechanical  alignment  of  the  experiment  consisted 
of  two  requirements.  The  first  requirement  was  to  position 
the  beam  axis,  the  target  normal  and  the  center  points  of 
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the  chopper  and  mass  spectrometer  ionizer  all  in  the  same 
horizontal  plane.  This  was  partially  accomplished  by 
adjusting  the  platforms  that  supported  the  mass  spectro¬ 
meter,  target  and  chopper,  as  seen  in  Pig.  8,  by  means  of 
level  and  height  gauges  with  the  test  chamber  removed. 

The  target  support  platform  served  as  a  datum.  The  second 
alignment  criterion  was  that  the  beam  axis,  target  surface, 
center  points  of  the  chopper  and  ionizer  should  also  lie 
in  a  common  vertical  plane.  Again,  the  target,  chopper 
and  ionizer  were  aligned  with  the  test  chamber  removed 
and  with  the  aid  of  two  plumb  lines  suspended  from  the 
mass  spectrometer  platform.  The  axes  of  rotation  of  the 
mass  spectrometer  and  target  should  coincide;  this  was 
accomplished  by  conti  mously  adjusting  the  platforms' 
support  so  that  the  edge  view  of  the  target  was  always 
viewed  as  a  line  through  the  centerline  of  the  ionizer. 
While  the  procedure  may  sound  complicated,  many  of  the 
alignments  were  incorporated  into  the  design.  With  the 
source  and  test  chambers  assembled,  the  alignment  of  the 
source  with  the  rest  of  the  components  completed  the 
procedure.  With  the  aid  of  a  telescope,  the  source  tip 
orifice  was  fi*st  positioned  at  the  center  of  the  then- 
aligned  ionizer;  the  source  collimator  was  then  centered 
around  the  source  orifice.  A  picture  of  the  aligned 
assembly  as  viewed  through  the  test  chamber  access  opening 
was  shown  previously  in  Fig.  9.  The  mass  spectrometer 
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translation  and  rotation  potentiometers  were  calibrated 
against  the  known  gear  ratios  once  the  alignment  was 
completed  whereas  the  target  rotation  potentiometer  was 
calibrated  against  the  mass  spectrometer. 

3.  Typical  Operating  Procedure 

A  set  operating  procedure  was  always  followed 
to  ensure  both  data  consistency  and  the  attainment  of 
steady  state  conditions  of  the  target  in  view  of  the 
foreseen  dependence  of  the  reaction  r  fce  on  previous 
target  hist  cy. 

The  system  was  maintained  in  vacuum  whenever  possible 

to  minimize  contamination  of  the  target  and  electron 

multiplier  dynode  surfaces.  When  the  target  was  first 

—7 

installed,  it  was  degassed  under  vacuum  (<  10  torr)  at 
1800°K  for  over  24  hours.  The  oxygen  beam  was  then  turned 
on  with  the  target  at  about  1300°K  (the  temperature  for 
maximum  reaction  probability)  for  at  least  48  hours  to 
activate  the  surface  after  annealing  and  before  any 
measurements  were  carried  cut.  This  procedure  of  degassing 
and  activation  was  repeated  whenever  vacuum  was  broken  and 
after  the  subsequent  pump-down. 

A  typical  day  of  experimenting  began  by  switching  the 
test  chamber  baffle  from  Freon  12  to  .liquid  nitrogen  cooled 
and  turning  on  all  electronic  equipment,  in  particular  the 
mass  spectrometer  controller  and  pulse  generator,  for 
warm-up.  This  was  followed  by  the  slow  process  of 


simultaneously  heating  up  the  source  and  adjusting  the 
source  leak  for  a  hot  beam  experiment.  The  source  inner 
wall  temperature  was  continuously  monitored  through  the 
orifice  by  the  optical  pyrometer  viewing  through  the  6!' 
diameter  view  port  on  the  test  chamber  cover  plate.  A 
spectral  emissivity  of  0.20  at  =  0.65y  for  the  thoria, 
as  furnished  by  the  manufacturer,  was  used  for  its 
temperature  calibration.  A  source  heat-up  time  of  approxi¬ 
mately  4  hours  was  required,  after  which  all  other  equip¬ 
ment  had  well  reached  their  steady  state  <  P'  rating 
conditions. 

The  incident  beam  was  first  monito-  td  to  check  all 
equipment  as  well  as  to  record  its  intensity  and  0/C>2 
ratio.  The  planned  experiment  was  then  carried  out  with 
the  Bayard-Aipert  imitation  gauges  turned  off  about 
15  minutes  prior  to  the  data  taking;  this  step  was  taken 
since  it  was  found  that  the  interaction  of  background 
oxygen  with  the  hot  filament  gave  rise  to  an  increased 
m/q  =  28  background  signal,  apparently  due  to  the  CO 
produced  by  the  reaction  of  the  impurities  in  the  tungsten 
filament  with  the  oxygen.  During  the  experiment,  the 
functioning  of  all  equipment  was  continuously  checked  and 
the  quality  of  vacuum  was  monitored  by  the  background 
count  rates  from  the  DSC.  At  the  end  of  an  experiment 
(which  at  times  ran  for  more  than  24  hours  in  a  stretch) , 
the  source  was  slowly  cooled  down  (over  a  time  span  of 


2  hours)  and  the  target  was  kept  at  1300 °K  with  the  beam 
on.  The  baffles  were  switched  back  t'  Freon  12  cooling 
and  the  whole  system  remained  in  vacuum. 

4.  Data  Processing 

The  DSC  has  the  canability  (with  the  addition  of 
a  plug-in  Serial  Converter  circuit  board)  to  transmit 
data  to  any  serial  digital  device  such  as  a  Teletype 
printer,  a  tape  punch,  or  a  time-shared  computer  for  fast 
or  real-time  pre-programmed  data  processing.  Due  to 
economics  and  the  relatively  small  amount  of  data,  all 
data  were  recorded  by  hand.  Much  of  the  data  reduction 
was  done  with  the  aid  of  a  Wang  Laboratories,  Inc.  Model 
380  programmable  calculator  and  a  Panasonic  Model  850 


electronic  calculator. 
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CHAPTER  IV 

RESULTS  AND  DISCUSSIONS 

A.  Introduction 

The  experiment  was  carried  out  chronologically  :.n  two 

phases:  the  first  phase  concerned  wi/th  the  C-C>2  reaction 

whereas  the  latter  was  aimed  at  the  C-0  reaction.  Though 

the  C-C2  reaction  had  been  studied  previously  using  a 

molecular  beam  technique/  '  the  present  experimental 

set-up  with  its  unique  capability  of  varying  several 

important  experimental  parameters  made  it  worthwhile  to 

further  investigate  the  reaction.  In  particular,  only 

29 

one  previous  reference  considered  the  reaction  product 
angular  distribution  and  therefore  a  check  on  that  work 
is  worthwhile.  The  effects  of  beam  energies  and  incident 
angles  on  the  oxidation  rates  may  be  helpfuj  in  understanding 
the  reaction  mechanisms.  In  addition,  precise  and  consistent 
C>2  data  are  required  for  the  specific  subtraction  process 
intended  to  deduce  the  0  results,  as  discussed  in  Section  II-D. 

B.  C>2  Scattering  from  Graphite 

The  angular  distribution  of  scattered  02  from  the 
graphite  target  was  the  first  measurement  made  for  this 
gas-surface  interaction.  One  of  the  purposes  of  this 
measurement  was  to  evaluate  the  detection  scheme  and,  in 
particular,  to  test  such  theoretical  prediction  as  the 
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signal-to-noise  ratio  after  integration  with  respect  to 
integration  time,  as  shown  in  Eq.  (2.13).  The  system 
behaved  quite  as  expected.  For  example,  a  signal 
fluctuation  level  of  approximately  25  percent  was  ob¬ 
tained  after  each  integration  period  of  5  minutes  for  a 
given  operating  condition  (i.e.,  at  given  Rg  and  Rj^  count 
rates)  whereas  prolonging  the  integration  period  to  30 
minutes  resulted  in  a  signal  fluctuation  level  of  approxi¬ 
mately  10  percent,  thereby  proving  the  square  root 
dependence  of  the  (S/N)^nt  ratio  on  the  integration  time. 

Results  of  the  angular  distribution  of  scattered  c>2 
are  shown  in  Figs.  26-28.  Theoretical  distributions  based 
on  a  cosine  flux  distribution  are  plotted  for  comparison. 

It  should  be  pointed  out  that  the  theoretical  cosine  dis¬ 
tribution  deviates  from  that  of  both  a  point  and  infinite 
target;  this  is  due  to  the  target  edge  corrections  for 
the  present  target-to-detector  geometry.  An  exact  analysis 
of  the  target  edge  corrections  is  presented  in  Appendix  A. 
In  any  case,  all  theoretical  angular  distributions  are 
plotted  out  with  the  present  target-to-detector  geometry 
in  mind.  Figs.  26  and  27  show  the  angular  distributions  of 
scattered  room  temperature  0n  at  incident  angles  of  45°  and 
60°  respectively.  The  target  temperatures  were  essentially 
all  about  1400 °K,  As  can  be  seen,  the  scattering  patterns 
resembled  those  of  diffuse  scattering  independent  of  the 
incident  angles  tested.  Similar  results  were  obtained  by 
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Shih,  who  used  nearly  monoenergetic  beams  at  an  incident 
angle  of  45°  with  target  temperatures  varying  from  300°  to 
1760°K.  In  particular/  Shih  also  found  no  apparent  change 
of  the  angular  distribution  due  to  either  target  tempera¬ 
ture  or  target  history.  This  behavior  can  be  attributed 
to  the  roughness  of  the  target  surface  used  in  this 
experiment  (incidentally/  exactly  the  same  type  of  isotropic 
graphite  was  used  by  Shih) .  If  a  molecule  makes  several 
collisions  with  a  roughened  surface  before  being  re-emitted, 
the  angular  distribution  of  the  re-emitted  flux  should-*  be 
both  insensitive  to  the  incident  angle  and  correspond  to  a 
random  pattern  (i.e.,  a  cosine  distribution).  Similar 
arguments  hold  for  the  insensitivity  of  the  scattered 
distribution  in  regard  to  the  incident  beam  temperature. 

This  is  shown  in  Fig.  28  where  a  beam  temperature  of 
1370°K  was  used;  again,  no  apparent  dependence  on  the 
incident  angle  was  observed. 

The  "Maxwellian"  beam  temperature  of  the  re-emitted 
C>2  at  Tt  =  1370°K  was  measured  to  be  about  1000°K  (i.e., 

<J>  =  12°)  regardless  of  detector  angle.  This  implied  an 
incomplete  energy  accommodation  of  0^  on  hot  graphite. 
However,  such  an  interpretation  should  be  taken  with  some 
reservation  due  to  the  inherent  insensitive  nature  of  the 
phase  shift  measurement  and  the  experimental  accuracy,  as 
discussed  previously  in  Sections  II-D  and  III-B.  At  the 
operating  target  temperature  of  1400°K,  an  inaccuracy  in 
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the  phase  shift  measurement  of  only  2°  could  raise  the 
possibility  of  a  fully  accommodated  oxygen  molecule. 
Nonetheless,  results  of  Shih  indicate  that  the  energy 
accommodation  for  C>2  adsorption  on  graphite  is  less 
than  unity  for  a  hot  target.  For  a  surface  temperature 
range  of  950°  to  1750°K,  Shih  found  that  the  "Maxwellian" 
beam  temper ature  of  the  re -emitted  0 2  was  insensitive  to 
the  target  temperature  and  remained  at  a  constant  value 
of  about  900°K  with  a  fluctuation  of  +  150°K.  The 
interpretation  of  the  measured  phase  lag  was  based, 
however,  on  the  assumption  of  a  Maxwellian  velocity 
distribution  of  the  re-emitted  flux  and  that  by  no  means 
was  confirmed. 

Previous  investigations  in  the  thermal  accommodation 

of  different  gases  on  graphite  revealed  similar  trends. 

93 

Meyer,  in  his  study  of  the  energy  exchange  between 

argon  and  a  hot  graphite  surface,  noted  the  existence  of 

a  critical  surface  temperature  above  which  the  temperature 

of  the  argon  leaving  the  surface  remained  constant. 

94 

Recently,  Siekhaus  et  al  studied  the  energy  of  02 
scattered  from  pyrolytic  graphite  using  molecu  .ar  beams. 

In  the  target  temperature  range  of  1000°  to  2200 °K,  they 
found  a  corresponding  re-emitted  02  Maxwellian  temperature 
range  of  700°  to  900°K,  again  indicating  an  incomplete 
energy  accommodation  and  its  weak  target,  temperature 
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dependence.  Furthermore,  the  beam  temperature  (from 
300°  to  1300 °K)  was  found  to  have  no  effect  on  the 
observed  phase  shifts.  This  was  attributed  to  trapping 
of  the  gas  on  the  surface  for  a  sufficiently  long  period 
such  that  the  initial  state  is  forgotten;  the  result  was 
apparently  in  agreement  with  the  present  measured  02 
scattering  patterns  in  regard  to  both  incident  beam 
angles  and  temperatures. 


Reaction  Product  and  Its  Angular  Distribution 


For  both  O  and  02  reactions,  CO  was  found  to  be  the 
only  detectable  reaction  product.  The  C02  signal, 
barely  observable  in  the  background  count  rate,  was 


estimated  to  be  at  least  two  orders  of  magnitude  lower 


than  that  for  CO  and  a  quantitative  measurement  of  the 


signal  level  was  simply  beyond  the  sensitivity  of  the 


present  detection  scheme.  This  finding  was  in  agreement 

with  that  of  Shih,  who  used  only  molecular  oxygen  with  a 

beam  intensity  almost  two  orders  of  magnitude  higher  than 

the  present  case.  Both  CO  and  C02  have  been  observed  by 

previous  investigators  as  products  of  graphite  oxidation. 

However,  as  the  target  temperature  is  increased  to  above 

1000 °K  for  the  C-02  reaction  (the  temperature  range  of 

the  present  experiment  covers  from  1000°  to  17CC°K) ,  the 

ratio  of  CO  to  C02  is  drastically  increased  with  CO  as  the 

dominant  product  (at  least  one  order  of  magnitude  higher 
26  2  8 

than  C02 ) .  '  C02  detection  from  previous  experiments. 
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a  number  of  which  were  reviewed  in  Section  I-B,  can  be 
partially  explained  by  the  presence  of  gas-gas  collisions 
where  CO  in  the  gas  phase  is  oxidized  to  CO2  by  the  back¬ 
ground  O2.  Such  a  mechanism  is  clearly  ruled  out  in  the 
present  case  as  well  as  similar  experiments  '  in  which 

the  background  pressure  is  of  the  order  10  torr. 

95 

Batty  and  Stickney  have  analyzed  the  chemical 
reactions  of  gaseous  02  with  solid  surfaces  under  conditions 
of  high  temperature  and  low  pressure.  Using  equilibrium 
considerations  and  without  assuming  a  detailed  kinetic 
model,  they  formed  a  thermodynamic  model  of  reactive  gas- 
surface  scattering  that  correlates  many  features  of  existing 
experimental  data.  In  particular,  for  the  C-02  reaction, 

CO  was  predicted  to  be  the  dominant  product  for  target 
temperatures  above  800 °K  with  a  O^-to-CO  ratio  of  less 
than  10 

There  is  scarcely  any  data,  in  particular  the  reaction 
product,  on  the  C-0  reaction  at  the  temperature  range 
covered  in  this  experiment;  however,  those  at  lower  target 
temperatures^  indicated  that  CO  was  also  the  principal 
product.  Subsequently,  all  future  referrals  to  products 
of  reaction  are  concerned  solely  with  CO. 

The  angular  distributions  for  CO  desorption  off  the 
target  are  shown  in  Figs.  29  and  30  for  different  source 
temperatures  and  beam  incident  angles.  Again,  the  target 
temperature  was  maintained  at  about  1400 °K  for  all  runs. 
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This  temperature  was  chosen  principally  because  the 

reaction  probability  was  found  to  exhibit  a  maximum 

(thereby  contributing  to  maximum  CO  signals  and  shortening 

29 

of  integration  time)  at  TT  -  1400°K.  In  addition,  Shih 
found  that  the  angular  distributions  were  independent  of 
target  temperatures.  Since  a  typical  set  of  angular 
measurements  required  an  operating  time  of  up  to  10  hours, 
variation  of  the  target  temperature  was  not  considered. 

Fig.  29  shows  the  angular  distributions  of  the  CO 
flux  from  the  target  for  a  room  temperature  02  incident 
beam  at  incident  angles  of  45°  and  60°.  The  results  also 
covered  a  four-fold  increase  in  incident  intensity  (as 
discussed  in  Section  III-B) .  No  apparent  dependence  of 
the  angular  distribution  on  the  incident  angles  can  be 
seen.  It  must  be  conceded,  however,  that  the  incident 
angles  used  did  not  encompass  a  wide  range;  this  was  due  to 
both  geometric  constraints  and  a  severe  drop  in  signal 
levels  at  large  incident  angles  (since  the  incident  flux 
on  the  target  varies  as  cos  0^).  Nevertheless,  the 
incident  angle  range  did  represent  a  41%  change  of  incident 
beam  normal  momentum  and  a  22%  change  of  tangential  momentum. 
The  insensitive  nature  of  the  CO  angular  distributions  with 
regard  to  incident  angles  —  equivalently,  normal  and 
tangential  momenta  —  may  be  related  to  the  fact  that  the 
first  step  necessary  for  a  reaction  to  occur  at  the  surface 
is  adsorption  of  the  gas  molecules.  Once  adsorbed,  a 
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molecule  tends  to  have  a  lesser  "memory"  of  its  past 
history.  While  it  may  havr  an  effect  on  the  reaction 
probability  (in  the  sense  of  affecting  the  sticking 
probability) ,  the  incident  angle  should  have  a  minor  or 
no  influence  on  the  desorption  patterns  of  the  reaction 
products,  as  indeed  was  found  to  be  the  case. 

Similar  results  with  re  spec,  t  to  the  CO  desorption 
patterns  corresponding  to  different  beam  temperatures 
were  obtained,  as  shown  in  Fig.  30.  It  should  be  pointed 
out  that  two  of  the  source  temperatures  —  2200°  and 
2220°K  —  represented  dissociated  beams  with  0-to-02 
flux  ratios  of  approximately  one-to-two,  but  apparently 
no  effects  on  the  desorption  patterns  could  be  observed. 

The  CO  signals  in  these  two  cases  were  contributed  almost 
entirely  by  the  C-0  reaction,  since  O  was  found  to  be 
more  than  ten  times  as  reactive  as  02,  as  will  be  discussed 
later.  The  arguments  for  the  unchanged  CO  desorption 
patterns  suggested  earlier  for  the  case  of  room  temperature 
02  beams  may  also  be  applied  here.  However,  to  draw  the 
conclusion  that  the  reaction  sites  for  both  0  and  02  are 
similar  based  on  the  similarity  of  the  CO  desorption 
patterns  is  certainly  premature;  the  influence  of  the 
rough  surface  structure  of  the  present  target  could  very 
well  dominate  over  the  difference  in  desorotion  patterns, 
if  any,  arising  from  CO  formed  at  different  sites.  The 
larger  fluctuation  of  data  observed  at  0D  =  60°  (Fig.  30) 
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choice  of  isotropic  graphite  as  the  target  material 

(primarily  for  its  higher  reaction  pi  obability  than,  say, 

pyrolytic  graphite)  has  the  drawback  of  a  rather  rough 

and  uncharacteristic  surface,  resulting  in  both  a  possible 

broadening  of  desorption  peaks  and  an  insensitivity  of  the 

desorption  patterns  with  regard  to  incident  beam  parameters 

(such  as  beam  energy  and  incident  angle).  Clearly, 

additional  studies  must  be  made  before  an  understanding  of 

such  a  desorption  phenomenon  can  be  expected.  The  result 

does  have  an  immediate  significance  in  the  sense  that 

experimentally  measured  reaction  probabilities  (such  as 

2  8 

those  measured  by  Olander  et  al  ) ,  which  assumed  a  diffuse 
desorption  pattern  may  br  questionable;  the  corrections  of 
such  results,  if  necessary,  will  be  quite  straight-forward 
once  the  angular  desorption  distributions  for  those  carbon 
surfaces  are  known. 

It  should  be  pointed  out  that  the  measured  "Maxwellian" 
beam  temperature  of  the  desorbed  CO  corresponded  quite  closely 
to  the  target  temperature  within  the  accuracy  of  the  phase 
shift  measurement.  This  beh .-vior  of  a  fully  accommodated 
CO  particle  was  in  contrast  to  the  partially  accommodaced 
scattered  02*  Such  a  difference  is  quite  reasonable,  since 
the  formation  of  a  reaction  product  involves  a  "completely 
trapped"  oxygen  particle  and  most  likely  the  reaction  product 
desorbs  after  equilibration  with  the  target  surface. 
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D .  C-0^  Reaction 

Before  the  results  are  presented,  some  pertinent 
points  about  obtaining  such  measurements  will  be  mentioned. 
The  preliminary  treatment  of  the  target  before  any  measure¬ 
ments  was  discussed  previously  in  Section  III-C-3.  In 
addition,  a  waiting  period  of  up  to  one  hour  was  needed 
before  data  taking  for  each  change  of  target  temperatures. 
This  step  was  taken  to  ensure:  (1)  the  attainment  of 
thermal  equilibrium  of  the  mass  spectrometer  ionizer  com¬ 
ponents  due  to  the  change  of  heat  flux  as  the  target 
temperature  was  varied;  (2)  a  steady  state  condition  of 
the  reaction  surface. 

It  was  found  that  the  mass  spectrometer  ionizer 
performance  (such  as  the  controls  for  tuning  to  a  peak) 
was  affected  by  the  heat  load  on  it.  This  was  apparently 
due  to  the  differential  thermal  expansions  of  the  different 
materials  making  up  the  ionizer  components.  The  drift  in 
tuning  was  particularly  noticeable  when  the  target  was 
first  heated  up.  However,  once  the  ionizer  attained  thermal 
equilibrium  (approximately  one  hour  after  target  heat-up) , 
any  subsequent  variation  of  the  target  temperatures  (in 
the  range  of  1000°  to  1700 °K)  had  a  minimal  effect  on  its 
performance  —  an  adjustment  time  of  only  about  5  minutes 
was  required. 
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Shih  found  that  whenever  the  target  temperature 

was  varied,  a  finite  adjustment  time  was  necessary  in 

order  for  the  oxidation  rate  to  reach  a  steady  ouate  value. 

The  adjustment  time  depended  on  both  the  target  temperature 

range  and  the  direction  of  temperature  changes.  Such 

27  28 

behavior  has  been  reported  previously  '  as  a  "hysteresis" 
effect  when  a  cyclic  measurement  was  carried  out  over  times 
that  were  short  compared  to  the  adjustment  times.  Basing 
the  procedure  on  the  transient  measurement  results  of  Shih, 
the  experiment  for  steady  state  oxidation  rates  was  carried 
out  in  such  a  way  so  as  to  avoid  the  hysteresis  effect. 
Typically,  the  target  temperature  was  first  set  a  1000 °K 
and  then  increased  at  100°K  intervals  (up  to  TT  =  1700 °K) 
after  each  measurement  was  made.  For  1000°K  <  T^,  <  1300°K, 
a  waiting  period  of  30  minutes  was  observed  while  for 
Tt  >  1300 °K,  the  waiting  period  lasted  for  1  hour.  Such 
time  spans  between  measurements  seem  adequate  as  there 
was  little  change  in  oxidation  rates  with  longer  waiting 
periods,  -.astly,  once  a  set  of  results  was  obtained,  the 
target  was  brought  back  to  about  1400 °K  with  the  beam  on 
for  at  least  36  hours  to  "reactivate"  the  surface  before 
another  set  of  measurements  was  carried  out. 

The  reaction  probabilities  of  C>2  with  graphite  at 
different  beam  temperatures  and  incident  angles  are  shown 
in  Table  IV-1.  The  reaction  probability  is  previously 
defined  by  Eq.  (2.27). 


Table  IV-1 


Reaction  Probabilities  of  0_  at  Different 
Beam  Temperatures  and  Incident  Angles 


0. 

1 

tb(°k) 

V (cm/sec) 

V  COS  0^ 

(cm/sec) 

S 

60° 

300 

5.23  x  104 

2.62  x  104 

1.56 

x  10"2 

60° 

1430 

1.14  x  105 

5.70  x  104 

1.40 

x  10“2 

60° 

1875 

1.31  x  195 

6.55  x  104 

1.62 

x  10**2 

45° 

300 

5.23  x  104 

3.70  x  104 

1.66 

x  10-2 

45° 

1370 

1.12  x  105 

7.92  x  104 

1.50 

x  10-2 

Throughout  the  above  measurements#  the  target  temperature 

was  maintained  at  1385°K  for  comparison  purposes  (the 

variations  of  ^  with  target  temperatures  wili  be  discussed 

later) .  Within  a  5-fold  increase  in  beam  temperatures 

(i.e.#  beam  energies)  and  a  2  1/2-fold  increase  in  normal 

be  aim  velocity  components  (shown  in  column  4  of  the  table)  , 

the  reaction  probabilities  remain  virtually  constant  within 

experimental  accuracy  (approximately  15%  fluctuations)  and 

show  no  definite  dependence  on  either  parameters.  Olander 
28 

et  al,  using  pyrolytic  graphite#  also  found  the  reaction 
rate  to  be  insensitive  to  the  temperature  of  the  incident 
beam  (300°K  <  TD  <  1200°K)  for  the  basal  plane  reaction. 

D 

The  constant  reaction  probability  with  respect  to  beam 
temperature  implies  one  or  both  of  two  possibilities. 
Assuming  that  a  beam  molecule  makes  only  one  collision  with 
the  surface  before  reacting  or  scattering#  then  the  results 
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suggest  that  the  surface  reaction  kinetics  are  truly 
independent  of  the  thermodynamic  state  of  the  incident 
gas  molecules  such  as  translational  energy,  rotational 
and  vibrational  excitations.  On  the  other  hand,  if  a  beam 
molecule  made  several  collisions  with  the  surface  before 
reacting  or  being  re-emitted  (most  likely  in  the  case  of 
a  rough  and  uncharacteristic  surface) ,  then  the  molecule 
would  be  preheated  to  (or  nearly  equal  to,  depending  on 
the  thermal  accommodation  coefficient)  the  target  tempera¬ 
ture  regardless  of  the  incident  beam  temperature  and  a 
beam  heating  experiment  would  be  inconclusive.  Similar 
possibilities  hold  for  the  dependence  of  the  reaction 
probability  as  regard  to  angle  of  incidence.  In  view  of 
the  rough  surface  nature  of  the  present  target,  the  latter 
case  probaLly  prevails.  The  roughness  of  the  target 
surface,  apart  from  examination  under  a  microscope,  is 
further  confirmed  by  the  measured  re-emitted  oxygen  flux 
distribution,  which  is  a  diffuse  pattern  regardless  of 
beam  and  target  tempera-cures  as  well  as  angle  of  incidence, 
as  discussed  previously  in  Section  IV-B. 

Though  Olander  et  al  started  with  a  polished  pyrolytic 
graphite  sample,  they  noted  definite  texturing  of  the 
surface  due  to  oxidation  after  removal  from  the  vacuum 
system  following  several  months  of  experimentation.  Scanning 
electron  micrographs  of  the  target  before  and  after  the 
experimental  program  showed  that  the  relatively  high 
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oxidation  rates  represent  conditions  far  too  chaotic  to 
permit  formation  of  the  nearly  hexagonal  etch  pits  ob¬ 
served  in  the  low  temperature  oxidation  of  single  crystal 

45 

flakes  cf  natural  graphite.  The  actual  area  of  the 
surface  is  undoubtedly  considerably  greater  than  the 
superficial  area.  Accordingly,  it  is  difficult  to  contend 
that  the  reaction  occurs  upon  perfect  basal  planes  of  the 
graphite  lattice,  despite  the  fact  that  a  specimen  of  this 
orientation  was  utilized  —  a  result  that  could  lead  to 
similar  findings  with  the  present  work  in  regard  to  the 
insensitiveness  of  the  reaction  probability  with  respect 
to  beam  temperatures. 

The  reaction  probabilities  of  C>2  at  different  target 
and  beam  temperatures  are  shown  in  Fig.  31  in  the  Arrhenius 
plot.  It  can  be  seen  readily  that  the  independent  nature 
of  the  reaction  probabilities  with  respect  to  beam  tempera¬ 
tures  extends  over  the  entire  target  temperature  range 
covered  (1000°  to  1700°K) .  This  result  should  be  expected 
in  view  of  the  arguments  given  previously.  A  typical  set 
of  reaction  probabilities  measured  by  Shih,  using  the  same 
target  material  but  with  a  nozzle  beam  source  and  an  analog 
lock-in  detection  ocheme,  is  included  for  comparison.  It 
is  satisfying  to  see  that  the  results  agree  with  those  of 
the  present  case  even  though  the  incident  beam  flux  of  the 
two  independent  experiments  differs  by  almost  two  orders  of 
magnitude.  Such  a  first  order  dependence  of  the  reaction 
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rate  on  the  incident  pressure  (implied  by  th-u  constant 
reaction  probability  with  varying  incident  flux)  within 
the  measured  pressure  range  is  quite  reasonable  in  view 
of  the  lew  incident  flux  and  the  low  reaction  probabilities, 
which  suggest  that  the  low  reaction  probability  is  due  to 
a  surface  rsciction  limited  step  rather  than  some  other 
mechanisms  such  as  gas-to-surf ace  diffusion  limited  reaction. 
Olander  et  al  found  the  same  first  order  pressure  dependence 
for  the  reaction  probability  in  a  comparable  range  of  beam 
intensities.  It  should  be  pointed  out  that  the  results  of 
Shih  vary  about  30-40%  in  absolute  magnitude  in  comparison 
to  a  maximum  of  15%  fluctuations  in  the  present  case.  The 
present,  relatively  consistent  02  oxidation  data  is  required 
to  extrapolate,  with  confidence,  the  C>2  oxidation  rate  to 
higher  beam  temperatures  for  the  interpretation  of  dis¬ 
sociated  beam  results  (in  spite  of  the  feict  that  e  turns 

c2 

out  to  be  independent  of  TR  in  this  case) . 

Fig.  31  also  shows  that  the  reaction  probabilities  do 
not  display  simple  Arrhenius  behavior  (i.e.  ,  a  constant 
activation  energy  represents  by  a  constant  slope  in  the 
Arrhenius  plot)  and  exhibit  a  maximum  at  T^,  =  1350°K  regard¬ 
less  of  beam  temperatures.  The  existence  of  a  maximum  in 
reaction  probabilities  with  varying  target  temperatures  is 
perhaps  the  most  consistent  feature  of  various  previous 
data  for  graphite  oxidation  involving  both  different  forms 
of  graphite  and  incident  pressure  levels.  The  present  work, 
however,  is  probably  the  first  to  show  such  a  behavior  for 
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high  temperature  0 2,  since  the  previous  experiments  deal 
almost  exclusively  with  room  temperature  oxygen.  In 
general/  it  is  the  competition  between  generation  of  active 
sites  by  0 2  impingement  and  their  deactivation  by  thermal 
annealing  that  results  in  the  observed  maximum  in  the 
reaction  probability  with  varying  target  temperatures. 
Though  previous  results  show  a  qualitative  agreement  for 
the  shape  of  the  reaction  probability  curve  in  an  Arrhenius 
plot,  the  magnitude  and  location  of  the  maximum  reaction 
rate  is  frequently  inconsistent  even  among  graphites  of 
the  same  type  —  some  examples  were  given  in  Chapter  I 
in  the  review  of  past  graphite  oxidation  investigations. 

Perhaps  the  most  widely  accepted  model  for  correlating 

the  C-C>2  reaction  rates  at  various  02  pressures  and  carbon 

temperatures  is  the  semi  empirical ,  two-site  model  of  Nagle 

37 

and  Strick land-Cons table.  Without  considering  detailed 

surface  mechanisms  and  nature  of  reaction  sites ,  they 
arbitrar  .  Lv  classified  the  surface  reaction  sites  into  two 
groups:  a  more  reactive  A-site  and  a  less  reactive  B-site. 
The  surface  reaction  .aeohanisms  are  postulated  as  follows: 

k.  * 

1.  (a)  A  +  02  — A 

*  k7 

(b)  A  — a  +  2co(g)  (4-D 

ko 

2.  B  +  c2  — **  A  +  2C0  (g)  (4.2) 

kT 

3.  A  — B  (4.3) 


where  A  represents  a  surface  complex/ 

A/B  represents  a  more/less  reactive  site, 
and  k's  are  the  respective  rate  constants. 


(4.1)  provides  the  formation  of  a  surface  oxide 
complex  and  its  subsequent  desorption  as  CO,  without  a  net 
generation  or  depletion  of  active  sites.  The  apparent 
increase  in  reactivity  with  continuous  impingement  of  O 2 
on  the  surface  is  represented  by  the  generation  of  A-sit.es 
from  B-sites  in  (4.2).  The  deactivation  of  a  surface  due 
to  thermal  annealing  is  provided  by  (4.3).  It  is  the 
competition  between  the  last  two  processes  that  results  in 
the  observed  maximum  in  the  reaction  probability  with  vary¬ 
ing  target  temperatures.  Though  somewhat  crude,  the 
classification  of  the  surface  sites  into  two  groups  is  not 
unreasonable  physically.  It  is  will  documented  that  grain 
boundaries  and  lattice  defect  points  such  as  vacancies, 

edge  and  screw  dislocations  are  the  more  reactive  sites 
45 

for  oxidation. 

At  steady  state  surface  conditions,  the  Nagle  and 
Strickland-Cons table  model  results  in  a  specific  surface 
oxidation  rate  given  by 


-  kA  Po2  -| 

=  >-*z  Poj  X  '  ^  P°2  a'XJ 


2_ 

cm  sec 


(4.4) 


where  x  = 


kT  i 

a  +  r1- — > 


(4.5) 


=  the  fraction  of  A-sites  on  the  surface. 


144 


The  first  and  second  terms  in  Eq.  (4.4)  represent  the 
oxidation  rates  at  A  and  B  sites  respectively.  Nagle  and 
Strickland-Cons table  then  tried  to  correlate,  and  with  some 
success,  the  available  data  of  oxidation  rates  of  pyrolytic 
graphite  with  the  following  set  of  rate  constants: 

kA  *  20  exp (-15 ,100/T)  gm/cm^-sec-atm 

-3  2 

kfi  =  4.46  x  10  exp (-7,640/T)  gm/cm  -sec-atm 

kT  =  1.51  x  105  exp(-48,800/T)  gm/cm2-sec 

and  kz  =  21.3  exp(2,060/T)  atm-1  (4.6) 

58 

Radcliffe  and  Appleton  are  able  to  correlate  more 
extensive  soot  and  pyrolytic  graphite  oxidation  data  using 
the  two-site  model  of  Nagle  and  Strickland-Constable  with 
the  same  parameters.  In  particular,  agreement  is  good  for 
an  oxygen  pressure  range  of  2.5  x  10  ^  to  1.0  atm.  Recently, 
Park  and  Appleton'*'0^  have  measured  the  surface  oxidation 
rates  of  soot  in  a  shock-tube  over  the  range  of  temperature 
1700°  -  4000°K  and  of  pressure  0.05  -  13  atm  of  oxygen. 

Again,  .he  results  can  be  approximately  correlated  (within 
a  factor  of  two)  by  the  semi-empirical  formula  (Eq.  (4.4)) 
proposed  by  Nagle  and  Stiickland-Constable.  Furthermore, 
surface  temperatures  above  approximately  2800°K,  the 
reaction  rate  tends  to  again  increase,  after  reaching  a 
local  maximum  at  a  lower  temperature.  This  behavior  is 
predicted  in  the  two-site  model,  apparently  due  to  an 
increase  in  reactivity  of  B-sites  at  high  temperatures.  In 


Fig.  31,  there  is  an  apparent  flattening  of  the  reaction 

probability  curve  at  temperatures  above  15C0°K  that  may 

indicate  the  possibility  of  an  increase  in  reactivity 

with  further  rise  in  temperatures.  Unfortunately ,  an 

increase  in  surface  temperatures  to  above  1700 °K  is 

beyond  the  capabilities  of  the  present  experimental  set-up. 

Clearly,  it  would  be  interesting  to  see  whether  such 

behavior  actually  takes  place.  Strickland-Constable *s 
26 

earlier  work  in  measuring  the  oxidation  rates  of 
graphite  at  C>2  pressure  level  of  approximately  2.5  x  10  atm, 
however,  does  not  show  an  increase  in  reactivity  at  high 
temperatures  (>  2000°K)  as  predi  ced  by  the  two-site  model 
even  though  it  fits  quite  well  with  the  rest  of  his  data, 
including  the  location  of  the  maximum  reactivity  and  its 
absolute  magnitude.  One  other  possibility  that  may  lead 
to  the  apparent  increase  in  reactivity  observed  by  Park  and 
Appleton101  is  the  direct  sublimation  of  graphite, 

4 

significant  at  temperatures  above  3000°K  . 

While  the  Nagle  and  Strickland-Cc ns table  model  with 

their  chosen  rate  constants  seems  to  correlate  quite  well 

(within  a  factor  of  two)  for  high  pressure  measurements, 

substantial  disagreements  are  present  for  the  low  pressure 

data.  The  semi-empirical  formula  tends  to  underestimate 

-7  -9 

the  reaction  rates  for  02  pressure  range  of  10  to  10  atm 
by  approximately  two  orders  of  magnitude.  Furthermore,  it 
no  longer  predicts  the  characteristic  local  maximum  for 
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the  reaction  rate  at  p  <  10  atm,  as  shown  by  Shih. 

°2 

This  is  in  direct  contradiction  with  the  present  results 
-9 

(pQ  =  10  atm)  as  well  as  with  other  low  pressure 
2 

28  —7 

measurements,  such  as  those  of  Olander  et  al  (p  =  10 

°2 

atm)  ,  Shih  (p  »  io  ®  atm)  and  Duval^  (p  =  10“®  atm) , 

°2  2 
which  all  show  a  local  maximum  in  the  reaction  rate  even 

though  different  forms  of  carbon  were  used. 

One  may  take  the  view  that  the  apparent  differences 
between  the  Nagle  and  Strickland-Constable  predictions  and 
the  present  available  experimental  results  of  low  oxygen 
pressure  may  be  attributed  to  the  large  number  of  parameters 
(eight  in  all  —  one  pre-exponential  factor  and  one 
activation  energy  for  each  of  the  four  rate  constants) 
that  were  originally  chosen  to  fit  principally  the  Strickland 
Constable  data.  Obviously,  such  a  set  of  parameters  may 
not  be  the  optimum  for  correlating  results  that  vary  over 
ten  orders  of  magnitude  in  oxygen  pressure  levels.  Further¬ 
more,  recent  results  of  Park  and  Appleton at  relatively 
high  pressures  (0.05  atm  <  PQ  <13  atm)  indicate  that  the 
Nagle  and  Strickland-Constable  formula  also  tends  to  under¬ 
estimate  the  oxidation  rates  at  high  temperatures  (i.e. , 

Tt  >,  2200°K).  However,  in  view  of  the  current  lack  in 
fundamental  understanding  of  the  surface  oxidation  mechanism, 
an  attempt  to  adjust  the  parameters  in  order  to  obtain  a 
better  fit  for  all  the  experimental  results  will  not  be  made 


here. 


On  the  other  hand,  one  may  agree  that  though  the 
Nagle  and  Strickland-Cons table  formulation  of  the  oxidation 
reaction  is  basically  sound,  as  proved  by  the  experiments 
(such  as  those  of  Shin  and  Olander  et  al)  of  thermal 
annealing  and  oxidation-induced  generation  of  active  sites, 
it  nevertheless  neglects  the  detailed  dynamics  of  the 
surface  processes.  In  spite  of  the  fact  that  the  formula¬ 
tion  does  provide  a  fractional  order  dependence  of  ft  on 

p  (as  shown  by  the  first  term  in  Eq.  (4.4),  which  arises 
°2 

because  the  degree  of  coverage  of  the  A-sites  by  surface 
oxides  increases  with  increasing  oxygen  pressure  in 
accordance  with  the  Langmuir  isotherm) ,  which  was  experi¬ 
mentally  confirmed10^"  by  the  changing  of  the  reaction 
order  from  first  to  zeroth  as  the  oxygen  pressure  was 
increased,  at  low  pressures  a  somewhat  more  detailed 
account  of  the  surface  processes  (such  as  oxygen 
dissociation-recombination,  migration  and  diffusion)  may 
be  necessary.  The  success  of  the  Nagle  and  Strickland- 
Constable  model  in  fact  mey  be  due  to  the  relatively  fast 
reaction  rates,  especially  in  cases  of  high  pressure  levels, 
which  tend  to  mask  whatever  surface  phenomena  that  may  take 
place  in  a  time  scale  -t  compared  to  the  experimental 
observation  time.  However,  at  low  pressure  levels,  surfaces 
processes  are  all,  if  not  equally,  important  and  a  correct 
prediction  of  the  oxidation  rate  should  require  that  account 
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be  taken  of  all  such  surface  phenomena.  Since  the  Nagle  and 
Strickland-Constable  model  does  not  consider  any  detailed 
reaction  mechanisms/  any  prediction  of  the  oxidation  rate 
in  such  a  low  pressure  regime  may  well  be  expected  to  be 
in  error. 

It  is  impossible  to  formulate  a  detailed  reaction 
model  from  the  present  measurements  of  steady  state  oxida¬ 
tion  rates.  For  such  a  purpose,  it  may  be  helpful  if  the 
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incident  beam  chopping  scheme  of  Olander  et  al  is 
incorporated  into  the  design,  in  which  case  information 
such  as  residence  time  may  be  obtained.  However,  the 
resulting  experimental  complexities  and  corresponding 
reduction  in  signal  levels  may  be  prohibitive  in  practice. 

In  the  present  case,  the  two-site  model  of  Nagle  and 
Strickland-Constable  seems  to  be  adequate  to  explaining 
most  of  the  observed  phenomena  as  well  as  illustrating  some 
of  the  major  differences  between  0  and  0 ^  oxidation  as  will 
be  discussed  in  the  next  section. 

E.  C-0  Reaction 

1.  Reaction  Probability 

The  reaction  probability  of  0  at  different  target 
temperature  is  shown  in  Fig.  32  in  an  Arrhenius  plot.  The 
corresponding  0 2  reaction  probability  is  included  for 
comparison.  The  0  oxidation  rate  was  inferred  by  subtracting 
the  CO  signal  due  to  0 2  from  the  CO  flux  measured  for  a 
dissociated  beam  with  an  0-to-02  flux  ratio  of  approximately 


one-to-three.  An  average  maximum  02  reaction  probability 
_2 

of  1.6  x  10  was  used  for  this  purpose  and  the  waiting 
period  after  each  target  temperature  change  was  identical 
to  that  described  in  the  previous  section  for  the  measure¬ 
ment  of  the  C-O2  reaction  probability. 

Some  distinct  features  of  the  0  reaction  probability 

are  apparent  in  Fig.  32.  First  and  upmost  is  the  fact 

that  e  is  more  than  an  order  of  magnitude  higher  than  e 

o  o2 

at  corresponding  target  temperatures.  The  ratio  of  eQ  to 
cQ  varies  from  a  high  of  approximately  37  at  TT  =  1000 CK 
to  a  low  of  approximately  16  at  TT  =  1350 °K  within  the  target 
temperature  range  of  1000°  to  1700°K.  The  O  reaction 
probability  also  exhibits  a  non-Arrhenius  behavior  (i.e., 
a  variable  slope  in  the  Arrhenius  plot,  implying  a  variable 
activation  energy)  with  the  maximum  reactivity  at  T^,  =  1400°K, 
v:hich  is  within  50 °K  of  the  target  temperature  for  maximum 
C>2  reactivity.  However,  it  is  reasonable  to  state  that 
the  O  and  02  reactivities  peak  at  a  c r  mon  target  tempera¬ 
ture  within  the  experimental  accuracy.  The  slope  of  the  eQ 
curve  is  also  less  steep  than  that  of  eQ  ,  implying  that 
the  0  reaction  probability  is  only  slightly  target  tempera¬ 
ture  dependent.  In  fact,  wituin  the  target  temperature 
range,  eQ  varies  by  only  45%  as  compared  to  75%  for  . 

In  order  to  test  the  reliability  of  the  data  reduction 


scheme  and  to  ensure  the  consistency  of  the  measured  data 
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(in  particular,  the  high  absolute  value  of  the  0  reaction 

probability,  which  at  the  maximum  case  corresponds  to 

more  than  one  reacting  O-atom  for  every  four  such  incident 

particles) ,  the  reaction  probability  of  O  was  obtained 

for  different  source  conditions  that  represent  different 

0-to-02  flux  ratios  and  total  incident  pressures. 

Table  IV-2  shows  the  measured  e  for  different  incident 

o 

beam  and  target  conditions. 

Table  IV-2 


Reaction  Probabilities  of  O  at  Different 
Source  and  Target  Conditions 


Ps (torr) 

Ts(°K) 

vs 

tt(°k) 

e  / 
o  o2 

eo 

0.10 

2220 

0.49 

1355 

20.1 

0.32 

0.20 

2205 

0.31 

1300 

16.4 

0.26 

0.24 

2200 

0.27 

1390 

14.6 

0.24 

0.26 

2110 

0.14 

1380 

15.7 

0.25 

_2 

An  average  value  of  e  =  1.6  x  10  was  vised  in  calculating 

°2 

the  e  -to-e  ratio  in  column  5  and  the  target  temperatures 
o  o2 

are  all  within  a  range  of  +  50°K  of  the  most  reactive  target 

temperature  (i.e. ,  1350°K) .  No  attempt  is  made  to  account 

for  the  difference  in  e  in  view  of  the  same  small  target 

°2 

temperature  range  in  question.  It  can  be  seen  that  the 
measured  0  reaction  probability  is  reasonably  repeatable 
(within  a  maximum  fluctuation  level  of  25%)  by  beams  of 
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different  compositions  (which  vary  from  an  O-to-O^  flux 
ratio  of  0.14  to  0.49)  and  flux  levels  (which  cover  a 
2  1/2-fold  range  in  magnitude) .  The  apparent  higher 
reaction  probability  at  p  =  0.10  torr  (e  =  0.32  as 

b  U 

compared  to  an  average  tQ  of  0.25  for  the  rest  of  the  data) 

may  lead  to  the  suggestion  that  s0  is  pressure  dependent. 

However,  a  closer  examination  reveals  that  the  0  flux 

levels  in  all  experimental  runs  are  comparable  (e.g. ,  the 

0  intensity  at  pg  =  0.10  torr  is  within  25%  variation  cf 

all  the  remaining  0  intensities) ,  thereby  eliminating  .  ich 

a  possibility.  For  establishing  the  pressure  dependence 

of  a  reaction  probability,  a  variation  of  at  least  one 

order  of  magnitude  of  the  pressure  level  is  desirable. 

Furthermore,  in  view  of  the  low  pressure  level  in  the 

experiment  (which  means  a  low  surface  coverage,  as  there 

are  approximately  ICO  surface  carbon  atoms  for  every 

incident  oxygen  particle) ,  a  first  order  pressure  dependence 

of  the  carbon  atoms  removal  rate  is  likely.  The  variations 

of  eQ  may  be  attributed  to  the  experimental  accuracy, 

which  shows  a  25%  fluctuation  level  as  compared  to  15%  for 

e  .  Nevertheless,  the  results  clearly  establish  the  fact 
°2 

that  O  is  much  more  reactive  (by  more  than  ten  times)  than 

C>2  at  the  corresponding  target  temperature. 

To  the  best  of  the  author's  knowledge,  the  only  other 

available  data  for  high  temperature  oxidation  of  graphite 

3  10 

by  dissociated  oxygen  is  that  of  Rosner  and  Allendorf.  ' 


Using  a  microwave  discharge  tube  for  generating  dissociated 
oxygen  and  calibrating  the  O-atom  concentrations  in  the 
vicinity  of  the  graphite  by  a  N0o  light  titration  technique, 
Rosner  and  Allendorf  measured  the  oxidation  rates  of  both 
pyrolytic  and  isotropic  graphite  in  a  gas  fl^w  system. 

The  graphite,  in  the  form  of  a  circular  filament  and 
heated  resistively,  was  placed  in  the  gas  stream  and  the 
reaction  rate  was  calculated  from  microscopic  determina¬ 
tions  of  the  reduction  in  filament  diameter  after  known 
exposure  times.  A  typical  total  test  section  pressure  of 
1  torr  was  maintained  and  oxygen  diffusional  limitations 
were  apparently  ruled  out  on  the  basis  of  the  absence  of 
effects  of  flow  rate  and  carrier  gas  identi' 
observed  reaction  rates.  No  reaction  product  "  dtored. 

The  reaction  probabilities  for  isotrop.  ;e 

measured  by  Rosner  and  Allendorf  together  with  present 
results  are  shown  in  Fig.  33.  Apart  from  the  absolute 
magnitude  and  the  location  of  the  maximum  reactivity,  the 
essential  features  of  the  two  findings  are  surprisingly 
similar  in  view  of  the  completely  different  experimental 
techniques  and  oxygen  pressure  levels  (which  vary  by  almost 
five  orders  of  magnitude).  The  common  features  are: 

(1)  0  is  much  more  reactive  than  C>2  (by  approximately  one 
order  of  magnitude)  at  corresponding  target  temperatures; 

(2)  neither  0  nor  C2  exhibits  simple  Arrhenius  behavior 
with  the  reaction  probabilities  of  both  0  and  C>2  peak  at 
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a  common  target  temperature  within  the  temperature  range; 

(3)  the  surface  temperature  dependence  of  the  oxidation 
probabilities  by  0  is  much  less  pronounced  than  that  by  C>2. 

The  difference  in  reaction  probability  magnitude  may  be 
attributed  to  the  different  kinds  of  graphite  used  in  the 
two  experiments  —  the  more  porous  (and  hence  more  surface 
area  and  grai;.  boundaries)  graphite  used  by  Rosner  and 

3 

Allendorf,  with  a  density  of  1.73  gm/cm  as  compared  to 

3 

the  present  graphite  density  of  1.7S  gm/cm  ,  may  account 

for  its  higher  reactivity.  The  same  explanation  may  be 

valid  for  the  difference  in  maximum  reactivity  location, 

since  such  peak  behavior  is  presumably  due  to  the  competitive 

processes  of  thermal  annealing  and  oxideti  ‘-induced 

generation  of  active  sites  and  the  egui'  •; ~i  point  for 

the  two  processes  definitely  depends  on  •  *  .-  .xnd  of 

surface  structures.  The  reaction  probability  of  Rosner 

and  Allendorf  was  also  found  to  be  higher  than  predicted 

by  the  Nagle  and  Strickland-Constable  model  at  the  correspond- 

5  8 

ing  pressure  level.  Nevertheless,  Fig.  33  clearly 
establishes  the  consistency  of  some  dominant  features  of  the 
C-0  and  C-02  reactions  that  may  lead  to  an  understanding  of 
the  reaction  mechanisms . 

While  detailed  reaction  mechanism  propositions  are 
certainly  premature  based  on  tn<?  prece  ng  data,  the  large 
absolute  value  of  the  reaction  probability  and  its  weak 
target  temperature  dependence  strongly  suggest  that  0- 
oxidation  takes  place  along  a  different  path  than  does 
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C^-oxidation.  Since  CO  is  the  dominant  product  for  both 
O  and  O2  oxidation,  it  can  be  suggested  that  while  the 
reaction  rate  of  O2  is  limited  by  the  surface  migration 
and  dissociative  adsorption  processes,  the  O-oxidation 
occurs  via  a  Rideal-type  mechanism  in  which  direct  O-atom 
strikes  from  the  gas  phase  lead  to  formation  and  desorption 
of  product  molecules,  and  no  migration  along  the  surface 
to  reaction  sites  is  necessary.  To  test  such  a  view,  the 
so-called  "hysteresis"  behavior  of  the  reaction  probability 
in  regard  to  thermal  annealing  of  the  target  has  been  in¬ 
vestigated,  as  will  be  described  in  the  following  section. 
Rosner  and  Allendorf  also  suggested  that  the  surface 
structure  may  have  a  minor  effect  on  the  0  reaction 
probability  —  a  fact  that  they  proved  by  showing  the 
relatively  constant  0  reaction  probability  for  both 
pyrolytic  aud  isotropic  graphite.  However,  their  postula¬ 
tion  of  the  C-0  reaction  mechanism  in  which  a  gaseous 
O-atom  impacts  upon  an  adsorbed  O-atom  leading  to  the 
formation  and  desorption  of  a  C02  molecule  may  be  in  error, 
since  no  C02  is  observed  as  the  reaction  product  in  this 
case. 

2.  Hysteresis  Behavior 

To  further  investigate  the  influence  of  the  surface 
conditions  in  regard  to  the  reaction  rates,  a  series  of 
experiments  has  been  carried  out  in  which  the  thermal  anneal 
ing  of  the  surface  is  known  to  influence  the  O2  reaction 
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probability.  The  results  are  presented  in  Fig.  34,  in 
which  the  ordinate  is  not  the  reaction  probability  but 
a  normalized  CO  signal  flux  from  the  target  for  0  and  0^ 
separately  instead.  This  is  done  simply  to  present  a 
clearer  picture  and  to  give  a  common  datum  point  for  com¬ 
paring  the  0  and  0 2  results . 

The  sequence  of  data  taking  is  important  in  inter¬ 
preting  the  respective  curves  in  Fig.  34.  Firstly,  the 
curve  was  obtained  using  pure  0 2  at  a  beam  temperature  of 
1875°K.  Starting  at  a  low  surface  temperature  of  1000°K, 

Tt  was  increased  at  100 JK  intervals,  after  each  data  point, 
to  a  maximum  temperature  of  1700°K.  T^,  was  then  decreased 

to  1100 °K  at  150 °K  intervals.  The  waiting  periods  between 
points  for  the  increasing  target  temperature  sequence  are 
similar  to  the  ones  used  for  obtaining  steady  stc.te  02 
reaction  probability,  as  discussed  in  Section  IV-D;  namely, 

30  minutes  for  1000°K  <  T„,  <  1300°K  and  1  hour  for  T^,  >  1300°K. 
On  decreasing  from  1700°  to  1100 °K,  a  waiting  period  of 
also  1  hour  was  used  between  points.  A  similar  experimental 
procedure  was  carried  out  for  the  case  of  a  dissociated 
oxygen  beam.  The  data  reduction  for  the  0  results  is 
standard  except  when  dealing  with  the  decreasing  target 
temperature  regime  —  the  CO  signal  due  to  0 2  that  it.  sub¬ 
tracted  out  from  the  combined  signal  for  a  dissociated 
beam  is  based  on  the  previously  measured  02  signal  for  a 
decreasing  target  temperature  regime  rather  than  the  steady 
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state  C>2  signal.  This  ensured  that  the  CO  signal  due  to 
02  in  a  dissociated  beam  is  accounted  for  under  identical 
conditions  for  a  pure  02  beam. 

It  is  interesting  to  see  in  Fig.  34  that  on  decreasing 

the  target  temperature  from  1700°K,  the  CO  signals 

(equivalently,  the  reaction  probabilities)  due  to  both  O 

and  02  do  not  return  to  their  previous  values  at  the 

corresponding  target  temper atures ;  instead,  both  show  a 

decrease  in  magnitude.  This  "hysteresis-like"  behavior, 

upon  cycling  target  temperatures ,  has  been  witnessed  pre- 

29  28  27 

viously  by  Shih,  Olander  and  Duval,  all  using 
molecular  oxygen;  however,  this  is  the  first  time  that  0 
has  been  shown  to  exhibit  a  similar  behavior.  Such  behavior 
can  be  explained  along  the  two-site  model  of  Nagle  and 
Strickland-Constable.  At  high  surface  temperatures ,  thermal 
annealing  of  the  active  surface  sites  takes  place.  When  the 
target  temperature  is  lowered,  the  annealed  active  sites 
require  a  finite  time  for  them  to  be  regenerated  by  oxygen 
bombardment.  If  the  reaction  probability  is  measured  at 
a  time  when  the  surface  is  still  undergoing  adjustment, 
then  the  measured  f  v.’ill  invariably  be  smaller  than  the 
steady  state  value.  While  one  hour  seems  to  be  sufficient 
as  the  adjustment  time  for  the  increasing  temperature 
sequence,  it  appears  to  be  inadequate  for  the  decreasing 
te.nperaturc  direction,  as  shown  by  the  decrease  in  c, 
possibly  due  to  the  concomitant  reduced  oxidation  (therefore 
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site  regenerating)  rates.  The  readjustment  time,  of 
course,  depends  on  the  type  of  graphite  and  the  degree 
of  surface  annealing  which  is  highly  temperature 
dependent,  as  shown  by  the  decrease  in  e  at  high  surface 
temperatures.  For  the  present  isotropic  graphite,  a 
readjustment  time  of  several  hours  is  required  to  reach 
steady  state  at  TT  =  1300°K  after  previously  annealing 
the  surface  at  =  1700°K,  as  shown  by  Shih. 

While  an  apparent  hysteresis  exists  for  both  0  and 
C>2 ,  a  difference  in  the  degree  or  extent  is  definitely 
noticeable  for  the  two.  For  example,  at  a  target  tempera¬ 


ture  of  1300°K,  e  comes  back  to  within  80-85%  of  the 

o 

stationary  state  value  before  annealing  as  compared  to 
only  approximately  20%  for  eQ  .  The  significance  of  this 
result  is  either  or  both  of  two  possibilities:  (1)  the 
surface  structure  has,  a  relatively  minor  effect  on  the  O 
reaction  probability;  (2)  the  regeneration  rate  of  active 
sites  for  0  is  much  faster  than  that  for  C^.  Most  likely 
both  of  these  possibilities  contribute  in  view  of  the  much 
less  pronounced  dependence  on  surface  temperature  as  well 
as  the  higher  absolute  value  for  eQ. 

3.  Reaction  Mechanism 

Though  an  understanding  of  the  graphite  oxidation 


mechanism  is  far  from  compete,  the  results  in  this  experi¬ 
ment  nonetheless  provide  enough  information  for  the  postula 
tion,  at  least  qualitatively,  of  a  plausible  reaction  model 
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The  i-'XOco3se~  for  --02  an<*  c_c'  reactions  are  summarized 
in  tbo  following  mechanisms  respectively: 


0(g)  — »  0(a) 


CO  (a)  — »  CO  (g) 


(4.8) 


where  (g)  denotes  the  gaseous  state 

(a)  denotes  the  adsorbed  state 
* 

and  C  denotes  an  active  carbon  site. 

For  a  C-02  reaction/  the  first  step  is  02  adsorption, 
which  is  followed  by  surface  migration  to  an  active  site. 

02  then  dissociates  to  form  either  two  adsorbed  CO  or 
one  adsorbed  CO  and  one  adsorbed  0.  The  adsorbed  particles 
are  all  subsequently  desorbed.  This  mechanism  is 

28 

essentially  similar  to  that  proposed  by  Olander  et  al 
with  the  exception  that  the  adsorbed  0-atom  formed  by  02 
dissociation  either  desorbs  or  reacts  to  form  another  CO 
without  further  surface  migration.  The  surface  migration 
of  the  adsorbed  0-atom  to  another  reactive  site  is  'he 
mechanism  proposed  by  Olander  et  al.  This  step  is  unlikely 
in  view  of  the  high  reaction  probability  of  0  and  its 
relative  independence  on  both  target  temperature  and 
structure.  Since  02  tends  to  dissociate  and  react  around 


active  sites  (a  fact  shown  by  electron  microscopy 
studies  '  ) ,  the  most  probable  process  will  be  the 

formation  and  desorption  of  two  CO  molecules  at  the  same 
point,  as  represented  by  the  upper  branch  mechanism  in 
(4.7).  The  lower  branch  mechanism  of  (4.7)  is  possible 
but  probably  accounts  for  only  a  small  fraction  of  the 
adsorbed  0 2,  again  due  to  the  high  reaction  probability 
of  an  adsorbed  O-atom.  Unfortunate  whether  such  a 
mechanism  takes  place  or  not  cannot  be  verified  experi¬ 
mentally,  since  the  desorbed  0  signal,  if  any,  will  be 
buried  in  the  huge  0  (or  m/q  =  16)  signal  from  the  crack¬ 
ing  pattern  of  the  surface  scattered  C>2  molecules.  C02 
formation  has  not  been  considered,  since  it  was  found  to 
be  experimentally  insignificant,  thereby  leading  to  the 
present  postulation  of  02  adsorption,  migration  and 

dissociative  reaction  mechanisms.  The  role  of  the  active 
* 

carbon  sites,  C  ,  can  be  explained  along  the  two-site 
model  of  Nagle  and  Strickland-Cons table.  In  this  model, 
all  surface  sites  can  be  considered  as  possible  reactive 
sites  with  the  only  distinction  of  an  associated  degree 
of  reactivity  —  in  particular,  they  are  classified  into 
two  types:  the  more  reactive  A-site  and  the  less  reactive 
B-site.  The  balance  of  A  and  B  sites  is  considered  in 
Section  IV-D  and  will  not  be  repeated  here. 
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A  Rideal-type  mechanism  in  which  a  reaction  readily 
takes  place  after  the  adsorption  of  an  O-atom  is  proposed 
for  the  C-0  reaction,  as  shown  in  (4.8).  The  absence  of 
a  dissociation  process  probably  accounts  for  the  high  0 
reaction  probability  and  its  relatively  independent  nature 
on  surface  structure.  The  proposed  mechanism  is  in  dis¬ 
agreement  with  that  suggested  by  Rosner  and  Allendorf 
in  which  the  high  reaction  probability  of  0  is  attributed 
to  the  direct  impingement  of  0(g)  on  0(a)  that  leads  to 
the  formation  and  desorption  of  CO 2»  The  high  0  reaction 
probability  can  also  be  attributed,  in  part,  to  the  in¬ 
creased  sticking  probability  of  0  than  02,  as  shown  by 
its  lower  activation  energy  (at  surface  temperatures 
below  approximately  1300 °K,  the  0  and  C>2  data  can  be 
represented  by  straight  lines  with  corresponding  local 
activation  energies  of  about  7  kcal/mole  and  36  kcal/mole 

respectively).  As  to  the  question  of  what  constitutes  an 
* 

active  site,  C  ,  the  same  classification  of  surface  sites 
for  C>2  can  be  applied  here  with  a  modification  only  in  the 
reactivity,  namely 

A-sites:  e  =  1,  e  S  1 
o  o2 

B-sites:  e  <  1  -  0.2-0. 3,  e  <<  1. 
o  o2 

In  general,  ail  surface  atoms  are  fairly  reactive  for  0 
as  show,*  by  the  high  reaction  probability.  This  could 
account  for  the  absence  of  a  surface  migration  process  in 
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the  proposed  C-0  reaction  mechanisms  in  (4.8).  The 

suggestion  of  two  types  of  surface  sites  for  O-oxidation 

is  reasonable  in  view  of  the  observed  local  maximum  in 

reaction  probability  and  the  hysteresis  behavior  upon 

surface  annealing  (both  of  which  are  fairly  well  explained 

by  the  two-site  model) .  This  classification  of  surface 

sites  is  further  supported  by  studies  of  the  surface 

structures  after  exposing  single-crystal  graphite  to 

45 

streams  of  dissociated  oxygen.  Using  electron  microscopy, 
it  has  been  shown  that  while  02-oxidation  forms  distinctive 
hexagonal  pits  at  crystal  defect  points,  O-oxidation 
generally  creates  a  background  of  conical  pits  over  the 
entire  basal  plane  with  a  higher  concentration  of  pits 
along  grain  boundaries.  The  implication  of  this  result 
is  that  while  O-oxidation  dees  not  necessarily  initiate 
at  crystal  defect  points,  it  tends  to  have  a  higher  reaction 
probability  at  such  points. 
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CHAPTER  V 

CONCLUSIONS  AND  SUGGESTED  FUTURE  WORK 

A.  Conclusions 

CC  has  been  observed  as  the  primary  reaction  product 

for  bovii  C-0  and  C-02  reactions  (C  being  isotropic 

graphite  in  this  case)  within  the  target  temperature 

range  of  1000°  to  1700°K  and  a  beam  pressure  of  about 
-7 

2  x  10  torr.  No  C02  was  observed  by  either  a  phase 
sensitive  detection  technique  or  monitoring  of  the  back¬ 
ground  CO 2  peak.  Both  CO  and  C02  have  been  observed  by 
previous  investigators  as  products  of  graphite  oxidation; 
however,  as  the  target  temperature  is  increased  to  above 

1200°K,  the  ratio  of  CO  to  C02  is  drastically  increased, 

2  8  29 

especially  for  experiments  in  the  low  pressure  range.  ' 
C02  detection  from  previous  experiments  can  be  partially 
explained  by  the  presence  of  gas-gas  collisions,  where 
CO  in  the  gas  phase  is  oxidized  to  C02  by  the  background 
02 .  Such  a  mechanism  is  clearly  ruled  out  in  the  present 
experiment,  where  the  background  pressure  is  approximately 
10-8  to  10-9  torr. 

Both  O  and  02  reaction  probabilities  exhibit  non- 

Arrhenius  behavior  and  show  a  local  maximum  in  the  target 

temperature  range.  A  qualitative  explanation  of  such  a 

behavior  can  be  presented  along  the  lines  of  the  two-site 

37 

model  of  Nagle  and  Strickland-Cons table , 


though  their 
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suggested  parameters  for  the  various  rate  constants  do 
not  seem  to  fit  available  data  in  the  low  pressure  range. 

It  is  the  competition  between  generation  of  active  sites 
by  oxygen  impingement  and  their  deactivation  by  thermal 
annealing  that  results  in  the  observed  maximum  in  the 
reaction  probability  with  varying  target  temperatures. 
Physically,  the  more  active  surface  sites  are  grain 
boundaries  and  lattice  defect  points,  such  as  vacancies, 
edge  and  screw  dislocations. 

The  reaction  probability  of  0  is  more  than  an  order 
of  magnitude  higher  than  that  of  0 2,  reaching  a  maximum 
of  almost  0.3  at  a  target  temperature  of  about  1350°K, 

It  is  also  weakly  target  temperature  dependent.  The  large 
absolute  value  of  the  reaction  probability  and  its  weak 
temperature  dependence  strongly  suggest  that  O-oxidation 
takes  place  along  a  different  path  than  does  02-oxidation. 
Since  CO  is  the  dominant  product  for  both  O  and  02  oxida¬ 
tion,  it  is  proposed  that  while  the  reaction  rate  of  C 2 
is  limited  by  the  adsorption,  surface  migration  and  dis¬ 
sociation  processes,  the  O-oxidation  occurs  via  a  Rideal- 
type  mechanism  in  which  the  surface  adsorption  of  an  O-atom 
from  the  gas  phase  leads  directly  to  the  formation  and 
desorption  of  a  product  molecule,  and  no  migration  along 
the  surface  to  a  defect  point  is  necessary.  This  view  is 
further  supported  by  a  surface  annealing  experiment  in 
which  it  was  found  that  while  annealing  at  high  temperature 
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(1700°K)  apparently  reduces  the  C>2  reaction  probability 

subsequently  measured  by  virtue  of  a  reduction  of  the 

defect  points  in  the  graphite,  it  simultaneously  has  a 

minimal  effect  on  the  0  reaction  rate.  Recent  observation 

64 

by  McCarroll  and  McKee  using  electron  microscopy 
techniques  indicates  that  while  0 2  attacks  graphite  at 
lattice  defect  points,  0  produces  a  general  background 
of  conical  pits  in  the  surface  —  a  result  that  also  lends 
credence  to  the  suggested  reaction  mechanisms. 

While  not  surprising,  it  is  disappointing  to  find 
*-hat  the  reaction  probability  of  02  does  not  vary  as 
functions  of  beam  temperature  (300°K  to  1875°K)  and  angle 
of  incidence  (45°  to  60°)  within  the  accuracy  of  the 
experiment.  This  behavior  can  be  attributed  partly  to  the 
rough  and  uncharacteristic  surface  of  the  graphite  being 
used  in  this  experiment.  If  a  molecule  made  several 
collisions  with  a  roughened  surface  before  reacting  or 
being  re-emitted,  then  the  oxygen  would  be  preheated  to  the 
target  temperature  regardless  of  the  incident  beam  tempera¬ 
ture  and  a  beam  heating  experiment  would  be  inconclusive. 
Similar  arguments  hold  for  the  insensitivity  of  the  reaction 
probability  as  regard  to  angle  of  incidence.  The  roughness 
of  the  surface  is  also  confirmed  by  the  measured  re-emitted 
oxygen  flux  distribution,  which  is  a  diffuse  pattern  regard¬ 
less  of  beam  and  target  temperatures  as  well  as  angle  of 


incidence. 
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The  flux  distribution  of  CO  off  the  target  is  sur¬ 
prisingly  lobular  in  the  normal  direction  and  can  be 
fitted  quite  closely  by  a  cosine  square  distribution 
independent  of  target  and  beam  temperatures  and  incident 

angle.  Lobular  desorption  patterns  have  been  observed 

96  9 8 

for  H2  and  D2  desorbed  from  metal  surfaces  '  and 

96 

although  Van  Willigen  suggested  that,  the  activated- 
adsorption  model  of  Leonard-Jones  may  be  used  as  the 
basis  for  developing  an  approximate  model  of  the  desorption 
process,  actual  agreement  between  the  observation  and 
theory  was  fair.  It  is  questionable  whether  Van  Willigen1 s 
model  can  be  applied  in  our  case  and  clearly  additional 
studies  must  be  made  before  an  understanding  of  such  a 
desorption  phenomenon  can  be  expected. 

B.  Suggested  Future  Work 

Though  a  great  deal  of  informat  .on  has  been  learned 

about  graphite  oxidation  recently,  the  reaction  mechanisms 

are  far  from  fully  understood.  There  are  still  apparently 

contradictory  results  (and  hence  different  conclusions) 

obtained  by  various  investigators  using  different  experi- 

13 

mental  techniques.  Some  of  the  discrepancies  may  no 
doubt  be  attributed  to  the  complexities  of  the  oxidation 
processes,  but  much  may  be  due  to  the  vastly  different 
crystalline  structures  of  the  carbon  samples  used  by 
different  investigators.  The  present  isotropic  graphite 
target  may  not  be  the  ideal  material  for  investigating  the 


importance  (or  irrelevance)  of  some  experimental  parameters 
for  the  oxidation  reaction  due  to  its  macroscopic  roughness 
as  well  as  its  uncharacteristic  crystalline  structure. 

For  example,  whether  the  gas  energy  and  momentum  have  any 
effect  on  the  reaction  probability  is  still  open  to  ques¬ 
tions.  One  way  to  clarify  such  problems  is  to  use  single¬ 
crystal  (limited  by  available  size,  however)  or  pyrolytic 
graphite  for  better  surface  characterization.  This 
approach  may  also  serve  as  a  more  consistent  basis  for 
comparison  of  results  among  various  work.  One  must  be 
cautioned  that  continuous  monitoring  of  the  surface 
structure  may  be  necessary,  since  the  surface  condition 
of  a  sample  is  liable  to  change  as  the  experiment  is 
carried  on. 

The  primarily  steady  state  measurements  of  the  present 

experiment  fail  to  obtain  the  rate  constants  of  such 

processes  as  th  anal  annealing  and  oxidation-induced 

surface  reactivation.  Transient  studies  under  various 

29 

surface  conditions,  such  as  those  performed  by  Shih, 

may  be  helpful  in  estimating  some  of  the  rate  constants. 

In  addition,  the  incident  beam  chopping  scheme  of  Olander 
2  8 

et  at  is  particularly  promising  in  correlating  results 
with  proposed  mechanisms.  These  approaches,  however,  are 
plagued  by  the  ever-present  problem  of  low  signal-to-noise 
ratio,  especially  in  the  present  case  of  a  relatively  low 
incident  beam  flux  for  dissociated  oxygen. 


APPENDIX  A 


Due  to  both  the  finite  size  of  target  and  detector 
and  the  relatively  short  distance  between  them,  the 
target  can  no  longer  be  considered  as  a  point  source 
for  the  measurements  of  the  spatial  distributions  of 
the  re-emitted  particles.  The  flux  measurement  as  a 
function  of  the  detector  angle  (measured  from  the  target 
normal)  requires  a  correction  that  arises  from  a  non- 
uniform  flux  across  the  chopper  collimator,  resulting 
from  purely  geometrical  effects  as  well  as  non-uniform 
incident  beam  flux  on  the  target.  Since  the  incident 
beam  intensity  is,  for  all  practical  purposes,  uniform 
across  the  whole  target  (which  subtends  a  source  angle 
of  less  than  6°),  as  can  be  seen  in  Fig.  24,  no  correc¬ 
tion  is  necessary  for  non-uniform  incident  f?ux  in  this 
case.  Corrections  due  to  geometrical  effects  are  necessary 
when  the  detector  view  includes  either  one  or  both  target 
edges,  which  results  in  deficiencies  in  the  chopper 
collimator  flux  when  compared  with  that  of  an  infinite 
target.  The  corrections  for  the  exact  present  target- 
-cor  geometry  with  the  identical  chopper  have  been 
dered  in  Ref.  29  and  what  follows  will  be  a  summary 
of  the  steps  to  arrive  at  the  theoretical  distribution 
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curves ,  as  shown  in  the  figures  for  the  angular  distribu¬ 
tions  of  the  scattered  or  reaction  product  particles  from 
the  target. 

Fig.  35  shows  a  schematic  diagram  of  the  target- 
detector  geometry,  which  is  taken  directly  from  Ref.  29. 

The  notation  is  quite  self-explanatory  except  perhaps  the 
following:  Xg  is  the  point  in  the  chopper  collimator  where 
the  detector  first  sees  a  target  edge;  AG^  is  the  lateral 
angle  measured  from  the  target  edge,  over  which  there  is 
no  flux  contribution  to  the  chopper  collimator  and  results 
in  a  maximum  flux  deficiency  at  the  edge,  i.e.  X  =  0;  0O 
is  the  angle  subtended  by  the  detector  at  any  point  in  the 
collimator  and  is  independent  of  X  for  d1<<£^.  All  these 
parameters  have  expressions  given  in  Ref.  66  on  which  the 
analysis  is  originally  based. 

When  the  detector  is  rotated  about  the  target  normal, 
it  begins  to  see  target  edges  which  results  in  a  deficiency 
in  the  detected  flux  when  compared  with  that  of  an  infinite 
target.  The  deficiency  terras  are  given  by 
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where  F/F0  is  the  chopper  collimator  flux  per  unit  chopper 
ccllimatcr  width  divided  by  the  infinite  target  value. 

The  effective  shutter  function,  A(Q^) ,  that  accounts  for 
the  flux  deficiency  due  to  the  target-detector  geometry 
for  the  present  half  sine  wave  shutter  function  of  the 
tuning  fork  has  been  evaluated  in  Ref.  29;  this  is  done  by 
weighting  the  half  sine  wave  function  (which  corresponds 
to  an  effective  shutter  function  for  an  infinite  target) 
by  Eq.  (A. 1)  and  integrating  over  the  entire  shutter  per¬ 
iod.  Such  i  procedure  will  not  be  repeated  here  and  the 
reader  is  referred  to  Ref.  29. 

The  theoretical  curves  for  the  detected  signal  for 
given  flux  distributions  can  now  be  calculated.  If 
F„  (dr)  is  the  angular  flux  distribution  per  unit  target 
area  for  an  infinite  target,  then  the  detected  flux  by  a 
fixed  size  detector  (as  in  the  case  of  the  mass  spectro¬ 
meter  ionizer)  is  %  f0(6i)/cos  0r  since  the  area  of  the 
target  seen  by  the  detector  is  'v  1/cos  @r.  For  the  case 
of  a  finite  target,  the  actual  detected  flux,  D(0r), 
would  be  less  due  to  target  edge  effects  and  the  reduction 
j.n  signal  level  is  proportional  to  the  ratio  of  the  areas 
of  the  effective  shutter  function,  A(©r)  and  the  true 
shutter  function,  A0,  which  is  simply  the  area  of  the  half 
sine  wave  in  this  case,  resulting  in 


The  pumping  speeds  of  the  diffusion  pumps  as  installed 
in  the  vacuum  system  were  measured.  This  was  necessary  be¬ 


cause  the  manufacturers'  supp.ied  values  are  for  those  of 
an  unbaffled  system. 

Gas  at  atmospheric  pressure  was  trapped  in  a  gradu¬ 
ated  burette  immersed  in  an  open  glass  beaker  containing 

D.C.  704  oil  (vapor  pressure  of  D.C.  704  at  room  tempera- 

_  8 

ture  is  approximately  10  torr) .  The  gas  in  the  burette 
was  leaked  into  the  vacuum  chamber  through  a  small  adjust¬ 
able  leak  valve  with  the  gas  si'oply  pressure  at  almost 
constant  pressure  (atmospheric)  since  the  rise  in  oil 
column  in  the  b  :_ette  due  to  the  gas  evacuation  was 
always  kept  to  be  less  than  5  cm  as  compared  to  8.6  meters 
of  oil  for  one  atmosphere.  From  a  measurement  of  both 
the  time  for  a  giver,  volume  of  gas  to  leak  into  the  cham¬ 
ber  and  the  resulting  rise  in  chamber  pressure,  one  can 


obtain  the  pumping  speed,  v»  in  liters/sec  from 


P„  Av 
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where  p„  =  gas  supply  pressure  in  torr 

o 

Av  =  volume  of  gas  evacuated  in  liters 

t  =  time  of  evacuation  in  sec 

Ap  =  change  of  chamber  pressure  in  torr 
c 

Different  gases  were  used  for  the  calibration  of  the  test 
chamber  16"  diffusion  pump  and  only  O2  was  used  for  the 
source  chamber  6"  booster  pump.  The  results  of  the  cali¬ 
brations  are  shown  in  Table  B-l. 

Table  B-l 

Pumping  Speeds  of  the  Diffusion  Pumps 


Pump 

Gas 

V  (liters/sec) 

V/V 

0, 

(m  /m) 
0  2 _ 

16"  dia. 

1730 

z 

1.00 

1.00 

(Test  Chamber) 

CO 

1955 

1.13 

1.07 

C°2 

1510 

0.87 

C.  85 

6"  dia. 

(Source  Chamber) 

°2 

380 

Under  free-molecule-f low  conditions  in  the  vacuum 
chamber  as  in  the  present  case,  the  pumping  speeds  for  dif¬ 
ferent  gases  vary  as  the  inverse  of  the  square  root  of  the 
gas  molecular  weight.  This  result  is  confirmed  by  comparing 
the  last  two  columns  in  the  table. 
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APPENDIX  C 


Ion  Gauge  and  Mass  Spectrometer  Calibrations 

The  two  ion  gauges  (GE  hot  filament  miniature  ion¬ 
ization  gauges)  that  were  used  to  measure  test  and  source 
chamber  pressures  respectively  were  calibrated  by  means 
of  a  Knudsen  flow  method  and  the  results  checked  independ¬ 
ently  by  a  McLeod  gauge.  A  small  chamber,  approximately 
6"  x  4"  diameter,  was  mounted  on  the  side  of  the  test  cham¬ 
ber  from  which  it  was  separated  by  a  thin  stainless  steel 
plate  with  a  3/4”  diameter  orifice.  An  adjustable  leak 
system  with  a  measurable  flow  rate,  similar  to  that  used 
for  the  pumping  speed  measurements,  was  connected  to  the 
small  chamber  with  its  pressure  monitored  by  the  ion  gauge 
to  be  calibrated.  The  pressure  in  the  small  chamber,  under 
Knudsen  flow  conditions,  is  given  by 

p  =  —  y/  2iTmkT  (C.l) 

where  A  is  the  orifice  area  and  F  is  the  total  flux  out  of 
the  orifice  which  in  steady  state  conditions  is  equal  to 
the  flow  rate  measured  by  the  gas  supply  system.  The  mea¬ 
sured  gauge  sensitivities,  Si  =  i  )  ,  for  different 

gases  are  shown  in  Table  C-l. 
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Table  C-l 


Ion  Gauge  Sensitivities  for  Different  Gases 


Gas 

Si  (Torr) 

Vs*, 

N2 

1.83 

£ 

1.00 

°2 

1.54 

0.82 

CO 

1.87 

0.99 

C°2 

2.06 

1.09 

The  gauge  sensitivity  for  N2  was  further  checked  by 

a  liquid  nitrogen  trapped  McLeod  gauge  which  showed  an 

agreement  to  within  10%  of  the  measured  chamber  pressure 

-6  -4 

for  a  pressure  range  of  3  x  10  to  10  torr.  The  cali¬ 
brated  ion  gauge  was  used  for  subsequent  calibration  of 
the  mass  spectrometer. 

The  mass  spectrometer  sensitivities  for  different 
gases  were  calibrated  by  leaking  the  gases  into  the  test 
chamber  through  the  beam  source  with  the  test  chamber  pres¬ 
sure  monitored  by  the  calibrated  ion  gauge.  The  ionizer 
electron  emission  current,  electron  accelerating  voltage 
and  the  electron  multiplier  voltage  were  set  at  1  ma,  70 
volts  and  2750  volts  respectively,  whereas  the  ion  draw- 
out  and  deflection  plates  1  and  2  were  set  nominally  at 
5.6,  3.2  and  5.3  units  respectively,  which  correspond  to 


! 
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a  maximum  output  signal  for  N2  (^  =  28) .  These  mass  spec¬ 
trometer  settings  were  kept  constant  throughout  the  exper¬ 
iment  to  provide  a  consistent  basis  for  comparing  the  var¬ 
ious  signal  levp  -  and  they  were  maximized  on  only  the  N2 
peak  since  28  is  about  the  medium  mass  number  for  all  the 
gases  of  interest  in  this  experiment.  Table  C-2  shows  the 
mass  spectrometer  sensitivities  of  the  different  gases  in 

terms  of  S1  as  defined  in  Eq.  (2.20).  The  units  for  S1 
ms  ms 

are  pulses/sec-amp-torr. 

Table  C-2 

Mass  Spectrometer  Sensitivities  for  Different  Gases 


Gas 

m/q 

si  Q015) 

ms 

s1  /(s1 

ms7  ms 

N2 

28 

2.05 

1.00 

°2 

32 

1.72 

0.84 

CO 

28 

2.26 

1.10 

C°2 

44 

1.56 

0.76 

Ionization  of  the  gaseous  molecules  by  electron  bombard¬ 
ment  also  results  in  a  dissociation  of  the  molecule  into 
different  components.  The  resulting  "cracking  pattern"  of 
each  gas  of  interest  in  this  experiment  is  tabulated  in 
Table  C-3;  the  cracking  pattern  is  shown  by  the  ratio  of 
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the  component  signal  level,  S^,  to  that  of  the  parent  gas 
signal  level,  S. 


Table  C-3 

Mass  Spectrometer  Cracking  Patterns  for  Different  Gases 


Parent  Gas 

Sc/s 

Vs 

Vs 

Ws 

°2 

0.087 

CO 

0.032 

0.016 

C°2 

0.030 

0.091 

0.076 

N2 

0.067 
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Fig.  1.  Degree  of  Dissociation  of  Oxygen  at  Thermal  Equilibrium 
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Fig.  3.  Timing  Diagram  for  Chopper  Mode  Operation  of  DSC  4 
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Fig.  11.  Source  Chamber 
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Interior  View  of  Source  Chamber 
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Fig.  15a.  Upper  Trace:  Reference  Signal 
(Chopper  Driver  Vcltage) 

Lower  Trace:  Measured  Shutter 
Function 

Time  Scale:  500  psec/cm 


Fig.  15b.  Measured  Shutter  Function 
by  Light-Photocell  Scheme 


Time  Scale:  200  psec/cm 
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Fig.  18.  Timing  Diagram  for  the  Phase  Adjusting  Scheme 
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Fig.  19.  Cross-Section  Diagram  of  Beam  Source 
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Fig.  21.  Schematic  Diagram  of  Gas  Inlet  System 
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Fig.  27.  Angular  Distribution  of  Scattoied  Room  Temperature 
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Fig.  30.  Angular  Distribution  of  Dasorbed  CO 
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Fig.  33.  Reaction  Probabilities  of  0  and  02  vs .  Target  Temperature 
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